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The absorption spectrum of furan vapor between 2200 and 1550A shows three band systems 
which overlap slightly. The second of these with a 0—0 band at 52,230 cm™ has sharp bands 
which could be analyzed to reveal three prominent totally symmetrical vibration frequencies 
of the excited molecule of 1395 cm=, 848 cm and 1068 cm™. Other characteristic frequency 
differences in the more diffuse bands are noted. Since the furan was found to be comparatively 
stable to photochemical decomposition, measurements of the intensity of absorption were made. 





HE ultraviolet absorption spectrum of furan 
at wave-lengths above 2100A has been 
studied by a number of investigators who have 
used both solution and the vapor phase.'! Some 
of these have found evidences of narrow shallow 
bands while others have observed only a broad 
stepout with no evidence of structure. All agree, 
however, that the absorption in this longer 
wave-length region is of very low intensity and 
thus one may safely conclude that any electronic 
transition here is of the forbidden type. The 
measurements also show that the absorption is 
increasing rapidly to a band of high intensity 
beyond the short wave-length limit of observa- 
tion and hence indicate that the lowest excited 
electronic states to which transitions are allowed 
are to be studied in the Schumann region. 
The present paper describes an investigation of 
the absorption spectrum of furan vapor at wave- 
lengths 2200-1550A which shows three distinct 





* W. Hartley and J. Dobbie, J. Chem. Soc. 73, 598 (1898) ; 
J. E. Purvis, J. Chem. Soc. 97, 1648 (1910); S. Menczel, 
Zeits. f. physik. Chemie 125, 161 (1927); I. Kaswagi, Bull. 
soc. chim. Japan 1, 145 (1926); E. C. Hughes and J. R. 
Johnson, J. Am. Chem. Soc. 53, 737 (1931); O. V. Fialkov- 
skaya, Acta Physico-chimica (U.S.S.R.) 9, 215 (1938); 
V. Henri and D.C. Foster (private communication). 


band systems of high intensity, in one of which 
a vibrational analysis was possible. The results 
of this study are to be correlated with an analysis 
of the vibrational frequencies of the normal furan 
molecule from infra-red data which will soon be 
published. 


EXPERIMENTAL 


Preparation and purification of furan 


Furan was prepared by heating the highest 
grade of furoic acid obtainable from Eastman 
Kodak Company and passing the vapor through 
a tower of sodium hydroxide pellets in an 
apparatus previously filled with nitrogen. A layer 
of glass wool placed over the furoic acid together 
with a careful regulation of the temperature 
prevented any appreciable quantity of the acid 
from subliming. The product was dried over lime 
and fractionally distilled in a one-meter Crismer 
column in a nitrogen atmosphere; the furan so 
obtained boiled at 31.3° at 760 mm and had a 
refractive index up”°=1.4214. The measure- 
ments of absorption spectra were made immedi- 
ately following distillation. The data were later 
duplicated with some furan kindly supplied by 
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Fic. 1. Absorption spectrum of furan. 
solution in hexane. 


Professor G. B. Kistiakowsky which had been 
carefully purified for heats of hydrogenation 
measurements.” 


Spectral measurements 


The absorption spectrum below 2200A was 
measured with a Hilger fluorite prism vacuum 
spectrograph and a hydrogen discharge as light 
source. Photographs were first made of a con- 
tinuously flowing stream of furan vapor at 
different pressures in a 60-cm absorption tube 
whose windows were protected by a layer of 
flowing nitrogen, using a technique which has 
been described.* Next, furan vapor was irradi- 
ated in the absorption tube for periods up to 
thirty minutes and the spectrum photographed 
at intervals of ten, twenty and thirty minutes in 
order to determine the extent of photochemical 
decomposition. Finally, since this substance was 
found to be quite stable to the action of the 
ultraviolet light, some measurements were made 
on the nonflowing vapor at different measured 
pressures. 

The Ilford Q plates which were used were so 
satisfactory as to permit of densitometer records, 

2M. A. Dolliver, T. L.Gresham, G. B. Kistiakowsky, E. A. 
11938) and W. E. Vaughan, J. Am. Chem. Soc. 60, 440 


’E. P. Carr and H. Stiicklen, J. Chem. Phys. 4, 760 
(1936). ; 
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the first of which were made on the Moll micro- 
photometer at the Massachusetts Institute of 
Technology used through the courtesy of Pro- 
fessor G. R. Harrison. Those made from the 
photographs of the flowing vapor were used to 
determine the position and relative intensities of 
the bands. The results were checked by measure- 
ments of the positions of the sharp bands with a 
Hilger comparator. In order to get a measure of 
absolute intensities, exposures of furan at dif- 
ferent pressures were compared photometrically 
with definitely timed exposures of the hydrogen 
lamp through the evacuated tube and extinction 
coefficients calculated for the wave-lengths where 
there was equal blackening of the two spectra. 


RESULTS 


Since measurements of the furan vapor in the 
quartz region have been recently made and will 
be published by Henri and Foster‘ only pre- 
liminary experiments in this region were made. 
No evidence of narrow bands was obtained 
either in vapor or in hexane solution. A set of 
15 narrow bands at wave-lengths 2500 to 2650A 
were found in the vapor by both Purvis and 


TABLE I. 








WAVE 
NUMBER 
IN CM! INTENS. 


48,890 
50,105 
51,200 


52,230 S. 0-0 band 
52,685 A V4 
53,078 S. 848 V2 
53,298 S. V3 
53,625 ‘ 1395 "4 
53,930 .W. 1700 1696 | 2 
54,136 , 1906 1916 vetvs 
54,455 , 2225 2243 vritv2 
54,690 . 2460 2463 vitvs 
55,010 S. 2780 2790 | 2 
55,195 2965 7 
55,317 3087 3091 vit2v2 
55,435 3205 3255 | 2vni+v 
55,540 3310 3311 vitvetys 
55,875 3645 3638 | 2v1+v2 
56,055 3825 3858 | 2n1+03 
56,290 4060 4103 2vnitvetr: 
56,410 4180 4185 | 3 

50,535 4305 4360 yeti 
56,660 4430 4486 | 2v1+2»2 
56,968 4738 4739) | 2vitvetys 
57,325 5095 5033 | 3vnit+ve 


ASSIGNMENT 
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4 Private communication, D. C. Foster. 
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Fialkovskaya! but the present measurements 
showed complete transmission up to 2300A at 
the highest vapor pressures measured® (400 mm 
in a 60-cm tube). The solution spectrum which 
shows a broad step-out from about 32,000 cm 
to 42,000 cm~ and a steep rise to a maxmium 
beyond 45,000 is reproduced in Fig. 1. 

In the Schumann region, three distinct band 
systems were noted, the first two of which may 
be seen in Fig. 2 which is a microphotometer 
record of a photograph of the region from 
48,000 to 59,000 cm-!. The bands on this are 
lettered for identification, and in Table I are 
given the corresponding frequencies with an 
estimated possible error of +10 cm which may 
be somewhat greater for the broader bands. The 
dispersion of the spectrograph is 200 cm per 
mm at 56,000 cm—. 

The first region of absorption consists of the 
three broad bands a’, b’ and c’ with separations 
of 1215 cm and 1095 cm. These resemble in 
appearance the bands of cyclohexadiene, 1,3 ® and 
cyclopentadiene’ which are to be observed at 
much lower frequencies in the quartz ultraviolet. 
Similar broad diffuse bands at slightly lower 
frequencies than those of the furan are to be 
observed in the spectra of all of the aliphatic 
hydrocarbons with conjugated double bonds 
studied in this laboratory’ and would seem to be 
characteristic of this grouping. There are also 
three step-outs or very shallow bands to be 
observed on the long wave-length side of these 
maxima which are not shown in the record. 
These step-outs were found at 45,960 cm-, 
46,760 cm—! and 47,480 cm~' with separations of 
800 cm! and 720 cm-, respectively. 

The second band system comprising bands a 
to x is characterized by sharply defined maxima. 
The densitometer record gives a somewhat mis- 
leading picture of the relative intensity of a as 
compared with the other strong bands because 
there is total absorption, and records of photo- 
graphs taken at lower pressures show that band a 
is much more intense than any of the other 
bands of this series. Thus it is assumed that a 
represents the vibrationless or 0—O transition. 





* Measurements of the vapor were made by Miss Elinor 
Sackter, Mount Holyoke College. 
°\V. Henri and L. Pickett, J. Chem. Phys. 7, 439 (1939). 
’E. P. Carr and H. Stiicklen, unpublished work. 


The next strongest band is e with a wave number 
difference of 1395 cm (v,;) which is interpreted 
as the totally symmetrical ring vibration in the 
excited molecule corresponding to 1499 cm™, the 
wave number of the very intense band of furan 
vapor in the infra-red. Bands c and d are assigned 
as fundamentals of wave number 848 cm (v2) 
and 1068 cm (v3), respectively. The band 
labeled 6 is so weak as to be almost questionable 
but it is tentatively assigned as a fundamental 
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Fic. 2. Microphotometer tracing of furan spectrum from 
48,000 to 59,000 cm showing second band system and a 
part of the first and third, together with exposure of 
hydrogen lamp for comparison. 








of 465 cm~'. It is not impossible to interpret 
this either as a transition from a higher vibra- 
tional level of the normal molecule (i.e., vo’’ — v;’) 
or as a part of the first band system of broader 
bands which overlaps this one. Thirteen of the 
remaining seventeen bands seem to be transitions 
to combinations or overtones of these three 
vibrational states, 71, v2 and v3, from the lowest 
vibrational level of the ground state. The 
strongest of these, somewhat surprisingly, corre- 
sponds to the combination frequency, v:+ 2, 
while that which has been identified as 2v,;+ v2 
is also abnormally strong. The overtones of all 
bands show a rapid and regular falling off of 
intensity so that no overtone higher than the 
third can be distinguished with certainty. Thus 
all of the prominent bands may be represented 
by the series 1/A = 52230+1395m+848n+ 1068p 
(where m=0,1,2,3, n=0,1,2, p=0,1,2). The 
assignment of the other bands, all of which are 
very weak, is less certain. Band k with a differ- 
ence of 2965 cm~ from a@ may measure a sym- 
metrical C—H stretching frequency of the ex- 
cited molecule (v;) and the other four bands be 
combination frequencies of vs and v;. They may 
represent two quanta of vibrational energy of 
any of the other possible nontotally symmetrical 
furan vibrations alone or in combination with 7. 
Table I gives the assignment considered the 
most probable in each case with calculated 








TABLE II. 








WAVE NUMBER OF LIGHT ABSORBED 


Beyond 45,600 cm 

46,250-55,550; beyond 59,120 
peak at 55,875 

Peak at 52,230; beyond 62,700 


LoG € 


2.2 mm 2.76 
.60 4 3.33 


.27 3.78 


PRESSURE 




















frequencies without anharmonicity corrections 
for comparison with the observed values. 

The best test of the probability of the assign- 
ment is in the comparison of vibrations with 
those of the normal molecule. If furan has a C2, 
symmetry and it is difficult to think of it other- 
wise, ® group theory methods predict? four totally 
symmetrical internal vibrations of the ring. 
Reitz” has tentatively selected these four fre- 
quencies from a study of depolarization factors 
of Raman lines of liquid furan as 1137 cm“, 
986 cm, 1483 cm and 724 cm," the last of 
which has however a depolarization factor of 
0.63. Two totally symmetrical C—H stretching 
frequencies are to be expected of a furan molecule 
of C2, symmetry; at least two frequencies have 
been reported between 3100-3200 cm. 

Since the moderately high intensity of the 
Schumann band system under consideration indi- 
cates that it represents an allowed transition, 
selection rules would permit totally symmetrical 
vibrational ‘changes only and experience has 
shown that only the ring vibrations are likely to 
appear with any intensity. The finding of three 
prominent vibration frequencies in the excited 
molecule which correspond to three of the four 
totally symmetrical ring frequencies of the 
normal molecule with a decrease in each case of 
about 100 cm~! seems strong confirmatory evi- 
dence for the present assignment. That the 
possible excited vibrational level of 465 cm is 
the counterpart of the normal 724 cm level is 
much more open to question not only because 


8G. B. Bonino who has made careful studies of the 
Raman spectrum of furan has suggested that the molecule 
may exist in a dimeric form with C2, symmetry to explain 
the anomalous Raman results. Atti x° Cong. Inter. Chim 
2, 141 (1938). 

®This number is found by applying the formula 
n=1/NZpRXpXp’, where N is the number of symmetry 
operations (4), Xz is the character of each and X,’ is the 
character as applied to each of the original coordinates of 
the 5 atoms of the ring. 

10 Reitz, Zeits. f. physik. Chemie B38, 381 (1938). 

11 Infra-red measurements of the vapor give correspond- 
ing frequencies which are slightly larger. 
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of the difference of these numbers but also 
because of the uncertainty that this is a real 
frequency of the excited molecule as previously 
discussed and because of the further question 
raised by the high depolarization factor of 0.63 
as to the validity of the assignment in the 
normal molecule. It may be that a frequency so 
near in magnitude to one of the other three that 
it is not resolved, is the missing fourth totally 
symmetrical ring vibration. 

All of the bands in this system have been 
considered to be transitions from the lowest vibra- 
tional level of the normal molecule. The reason 
for the absence of others is thought to be due to 
two factors. First, the energy of the totally 
symmetrical vibrational levels of the normal 
molecule is sufficiently high so that the number 
of molecules existing in them would be compara- 
tively low at room temperatures. The Boltzmann 
factor for a 724 level is 0.031, for a 986 level is 
0.009. Secondly, the overlapping band system 
on the long wave-length side which has been 
previously discussed would be likely to mask such 
weak bands. 

The third band system begins at the point 
where the hydrogen spectrum is discontinuous 
and this makes the bands difficult of measure- 
ment so that only strong bands may be deter- 
mined with certainty. There is a band of high 
intensity at 59,580 cm! and a second band of 
nearly equal intensity at 60,810 cm followed by 
an absorption region which includes a maximum 
about 62,800 cm~! and seems to have other less 
well defined peaks. There is also a shallow broad 
band at 58,765 cm= or 815 cm! to the long 
wave-length side of the first strong band. This 
may represent a transition from a higher vibra- 
tional state of the normal molecule although it 
does not correspond exactly to any of the vibra- 
tions which have been identified. The separation 
of the two strongest bands, 1230 cm, seems to 
be the vibration frequency corresponding to the 
1483 of the normal molecule and 1395 of the 
first excited state. 

The results of the intensity studies are em- 
bodied in the curve in Fig. 1, while Table II 
gives the pressures, exposure times, log ¢ values 
and absorption limits of a series of measurements 
in an absorption cell of 8.5 cm length. 

While these measurements are somewhat crude 
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and do not take into consideration the photo- 


chemical decomposition during the exposure 


(which preliminary experiments indicated is 
small) or certain minor optical errors, they are 
presented here as giving a semi-quantitative 
measure of intensity in the Schumann region 
where almost no intensity data exists because of 
the difficulties of measurement. It is of interest 
to note that the solution curve is slightly higher 
than the vapor curve in the region of over- 
lapping data. 

Mulliken” has predicted that there should be 
an absorption maximum for furan in the neigh- 
borhood of 58,500 cm-! corresponding to an 
N-V, or rather four N--V transitions of very 
similar frequency, which would have a total f 
value of 0.3. It seems possible that the third 
band system may represent this transition, 
although the observed f value, uncertain because 
of lack of knowledge of the high frequency limit 
of the band, would seem to be below this pre- 
dicted value. 

The second band system has the characteristics 
attributed to Rydberg series bands in that the 
bands are sharp with an intense 0—0 band and 


2 R. Mulliken, J. Chem. Phys. 7, 121 (1939). 


short vibrational progressions starting from a 
ground vibrationless state. This may be a 
transition of the Rydberg type superimposed on 
the other band system. Studies of the unsatu- 
rated hydrocarbons'* have shown that Rydberg 
series bands are to be expected in this part of 
the spectrum often superimposed on other bands. 

Examination of the photographs of the furan 
which had been exposed to the short wave- 
length radiation for periods up to thirty minutes 
showed that the furan was not destroyed in this 
time. Unlike the unsaturated hydrocarbons 
studied.in this laboratory no absorbing layer was 
formed on the windows of the absorption cell. 
However, there was evidence of the development 
of a substance with an absorption band from 
about 57,000 to 60,700 cm, a doublet with the 
stronger peak at 55,860, and a peak at 55,500 
during the exposure. 

I wish to acknowledge with appreciation the 
grant of a fellowship from the Lalor Foundation 
during which these measurements were made, 
and to thank Doctors E. P. Carr and G. B. 
Kistiakowsky for helpful suggestions. 


%E,. P. Carr and H. Stiicklen, J. Chem. Phys. 7, 631 
(1939). 
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Bond Strength and Potential Energy Relationships in Acetylene 


HyMAN HENKIN AND MILTON BurRTON 
Department of Chemistry, New York University, New York, New York 


(Received January 15, 1940) 


The spectroscopic value for the heat of formation of free CH and Herzberg’s value for the 
heat of sublimation of graphite lead to a value of ~137.9 kcal. for the C=C bond strength in 
acetylene. This value is employed in a study of the potential energy levels of that compound 
and is shown to be the only one consistent with the most recent studies of its absorption 
spectrum. The fact that the known energy levels of free CH fit into the scheme is further evi- 
dence for the 125.1 kcal. value for the heat of sublimation of graphite. 


renee widely divergent values for the C=C 
bond strength appear in the literature.! The 
value of 220.5 kcal. given by Norrish? is based 
upon an older value for the heat of combustion 
*Cf. G. K. Rollefson and M. Burton, Photochemistry 
and the Mechanism of Chemical Reactions (Prentice-Hall, 


Inc., New York, 1939), pp. 349 and 422. 
*R. G. W. Norrish, Trans. Faraday Soc. 30, 103 (1934). 


of acetylene and the value of 150 kcal. for the 
heat of sublimation of carbon, taking into 
account the reorganization from tetravalent to 
divalent carbon as is probable in free CH. That 
of ~200 kcal. given in an early work by Pauling® 


3L. Pauling, J. Am. Chem. Soc. 54, 3570 (1932). Cf. 
L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 (1933), 



























derives from a similar calculation, but takes no 
account of this reorganization energy. A later 
calculation** yielding a value of ~123 kcal. for 
the C=C bond strength, while using the more 
recent considerations of Herzberg for the heat 
of sublimation of graphite still employs a value 
for the heat of formation of free CH derived 
from that of methane instead of from spectro- 
scopic data. 

On the other hand the value of ~187 kcal. 
given by Price is based simply on the observation 
of predissociation broadening in the ultraviolet 
absorption spectrum of acetylene.* The broaden- 
ing places an upper limit on the strength of the 
bond concerned. The nature of the spectrum 
makes it improbable that a C—H bond is in- 
volved. Thus, the 187 kcal. value becomes an 
upper limit for the strength of the C=C bond. 
On the assumption that the true value lies close 
to this maximum,® Price was able to suggest a 
consistent system of potential energy curves 
which seemed to fit the known data on acetylene. 

We have re-examined the available data re- 
lating to the chemistry and spectroscopy of 
acetylene and have concluded that the C=C 
bond strength is ~137.9 kcal. In this discussion 
we retain the definition of bond strength pre- 
viously used:* the strength of a bond in a 
particular compound is the energy difference 
(per mole) between the lowest energy state of 
that compound and the summed lowest energy 
state of the products formed by rupture of that 
bond ; thus, the calculation is relative to two CH 
radicals in the II state.” It should be emphasized 
that the extremely low value we have obtained 
derives essentially from the fact that we have 
accepted Herzberg’s 125.1-kcal. value® for the 
heat of sublimation of graphite.® 


who give a value of 198.5 kcal. for all acetylenic hydro- 
carbons. 

s¢ 1. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, New York, 1939), p. 122. 

4W. C. Price, Phys. Rev. 47, 444 (1935). 

5In general such correspondence is not to be expected. 
Cf. M. Burton, J. Chem. Phys. 7, 1072 (1939). 

6M. Burton, J. Chem. Phys. 6, 817 (1938). 

7R.S. Mulliken, Rev. Mod. Phys. 4, 1 (1932). 

8G. Herzberg, Chem. Rev. 20, 145 (1937). Herzberg 
actually gives 124.1 kcal. See, however, reference 6. 

® Recent further investigations by Professor H. L. Johns- 
ton on the sublimation of graphite, as well as other con- 
siderations, lead him to the conclusion that the 125.1-kcal. 
value is correct. This work will be the subject of a shortly 
forthcoming paper. The authors wish to thank him for the 
opportunity to quote this result prior to publication. 
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In the estimation of the C=C bond strength 
the reliability of the following data has been 
assumed.!° 


CeHo+ (5/2)02>2CO2.+H,0(1) ; 


AH=-—311.1 kcal. (1)" 
C(gr) +O2=COz; 
AH = —94.030+0.011 kcal. (2)” 
H.+302=H.O(1) ’ 
AH = — 68.37 kcal. (3)"! 
C(gr,0°K) = C(gr) ; 
AH =0.2 kcal. (4)§ 
H2(g,0°K) = He ; 
, AH=1.9 kcal. (5)§ 
2H = (n)H2(g,0°K) ; 
AH = —102.48+.02 kcal. (6)'8 
CH=C+H; 
AH =80.0 kcal. (7)'4 
C=C(gr,0°K) ; 
AH = —125.1 kcal. (8)8 
In addition, we assume 
C.H2= C:2H2(g,0°K) ; 
AH = —2.2 kcal. (9)! 
From these data the relationship 
C2H2(g,0°K) =2CH: 
AH=137.9 kcal. (10) 


follows directly. This value may properly be 
compared with the C—C bond strength in 
ethane given variously as 72 kcal.? and 97.6 


10 Unless otherwise indicated all chemically stable com- 
pounds are referred to their standard states. C(gr) repre- 
sents graphite at 18°. Radicals and atoms are taken in the 
states corresponding to 0°K. 

1F, R. Bichowsky and F. D. Rossini, The Thermo- 
chemistry of the Chemical Substances (Reinhold Publishing 
Corp., New York, 1936). 

2F. D. Rossini, J. Chem. Phys. 6, 569 (1938); F. D. 
Rossini and R. S. Jessup, Nat. Bur. Stand. J. Research 
21, 491 (1938). 

13H. Beutler, Zeits. f. physik. Chemie B29, 315 (1935). 

4G. Herzberg, Molecular Spectra and Molecular Struc- 
ture. I. Diatomic Molecules (Prentice-Hall, Inc., New York, 
1939), p. 485. G. W. King, J. Chem. Phys. 6, 378 (1938), 
uses the value 78.9 kcal. Here AH is from spectroscopic 
data for Do. ; 

16 This is an approximation obtained by comparison with 
the cases of CH, and C2He (see reference 6). The error 
probably does not exceed +0.2 kcal. The same result is 
obtained by approximation fiom available thermal data. 
(Int. Crit. Tab.) 
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Fic. 1. Potential energy relationships in acetylene. The 
upper curves are from the data of Price (ref. 4). The lower 

rt of the normal curve is obtained from data on the 
undamental frequency and the moment of inertia. For the 
K curve, r, is given by Kistiakowsky (reference 17). 


kcal.!® The striking fact is that, if we accept the 
125.1-kcal. value for the heat of sublimation of 
graphite, we must conclude that the strength of 
the C=C bond is less than two times the strength 
of the normal C—C bond. 


SPECTRAL DATA AND POTENTIAL CURVES 


It will now be shown that the value of 137.9 
kcal. accords with data derived from a study of 
absorption bands in the near ultraviolet (at 
2400-1900A), first described by Kistiakowsky.” 
These bands have been further investigated by 
Jonesco® who was able to represent the band 
heads by the formula” 


v= v9 +1034.60’ — 9.60”, 1 


where v=0, 1, 2, 3---. From this equation he de- 
rived a value of ~204 kcal. as the convergence 
limit for the band heads. 

If this convergence limit is taken as the 
beginning of the region where the split into one 
normal (7II) and one excited (2A)CH radical 
occurs, then it may very simply be correlated 





®Trenner, Morikawa and Taylor, J. Chem. Phys. 5, 
203 (1937). 

’G. B. Kistiakowsky, Phys. Rev. 37, 276 (1931). 

8M, A. Jonesco, Comptes rendus 199, 710 (1934). 

In Jonesco’s article, the coefficient, 9.6, was omitted, 
Perhaps because of a typographical error. We have deter- 
mined it directly from his data. 
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to the C=C bond strength of normal acetylene. 
Mulliken’ has given 2.85 ev (65.5 kcal.) as the 
energy difference *A—*II for the CH radical. If 
we use the 137.9 kcal. value for the C=C bond 
strength, it follows that the dissociation in the 
first excited state should require 203.4 kcal., in 
very good agreement with Jonesco’s value of 
~ 204 kcal. 

The acetylene spectrum has been more thor- 
oughly investigated in the same region by Woo 
and his co-workers.”° They found that their data 
for what they designated as C bands could be 
represented by the relationship 


v=42,227.51+529.56v;' —2.4502;”, 2 


where v;’=0, 2, 4---. The convergence limit de- 
rived from this equation is 201.6 kcal. Since v; 
represents transverse vibration of the C=C bond 


t t 
V5: H-—-C=C-H. 
t v 


and the C bands in Woo’s classification are those 
with v,' =2, 


t ft 
vm: H-C=C-H, 
4 V 
the convergence limit corresponds to somewhat 
distended radicals with vibration frequencies 
slightly higher than in the normal state. This 
additional energy amounts to ~3 kcal. Thus 
using the data of Woo and his co-workers we find 
that the energy required to dissociate C2H:2 into 
CH(?II) and CH(?A) in their lowest states is 
~199 kcal. Considering the extrapolations in- 
volved and the assumptions made, the agreement 
with our ~ 203 kcal. is good. 

From the above data we may draw a more 
nearly correct picture of the C=C positions of 
the acetylene potential energy hypersurfaces. 
The assumption made by Price that the be- 
ginning of the predissociation in the spectrum 
corresponds to the bond strength has recently 
been shown! to be invalid in several cases where 
it might have been expected to hold; this 
assumption is at best an approximation, not 
holding well in a relatively large number of 
cases. However, Price’s value does fix an upper 
limit for the bond strength. 


20 Woo, Lui, Chu and Chih, J. Chem. Phys. 6, 240 (1938). 
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The potential energy curves which we have 
drawn for the various states of acetylene are 
shown in Fig. 1.2! The normal curve flattens out 
at ~137.9 kcal., not at 187 kcal. as previously 
suggested by Price.‘ He postulated that the K 
curve also approached 187 kcal. asymptotically 
and that the interaction of that curve with the 
upper excited states produces predissociation, 
observed both at 1520A and 1350A. He thought 
it likely that this curve might represent the 
upper state of the Kistiakowsky bands which are 
found in the near ultraviolet. From the work of 
Jonesco and of Woo and his co-workers, it now 
appears, however, that this level dissociates, not 
at 187 kcal., but at ~203 kcal., so that we 
cannot explain the observed predissociations as 
due to interactions of this state with the higher 
electronic levels. The Kistiakowsky bands re- 
quire a minimum in the K curve at slightly less 
than 120 kcal., as shown on the diagram; this 
value is fixed by the work of Jonesco and of Woo 
for the 0—O transition from this state to the 
ground state. 

To explain the predissociations at 1520A and 
1350A, we postulate a repulsive curve as shown 
on the graph, coming down to the asymptotic 
value for the ground state. The interaction of this 
curve with the =3 levels of Price’s Class I and 
II bands would account for the predissociation 
at 1350A and 1520A, respectively. There has 
been some question as to whether another pre- 
dissociation region occurs at ~1900A. Herzberg” 
has observed a continuum in this region at 
higher pressures. It may be that this is an 
example of an induced predissociation, which 
might be expected under such conditions. On the 
basis of these potential energy curves it would be 


21 Cf. G. B. B. M. Sutherland, Phys. Rev. 43, 883 (1933) ; 
Herzberg, Patat and Spinks, Zeits. f. Physik 92, 87 (1934). 
2G. Herzberg, Trans. Faraday Soc. 27, 378 (1931). 
According to Price (reference 4), Dieke was unable to 
detect predissociation at 1900A. 
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explained as an interaction of the K curve with 
the repulsive curve. This predissociation is 
possible, but not required, according to our 
diagram; that it has not been observed by all 
workers may be due to different pressures used 
by various investigators or perhaps even to the 
occasional presence of paramagnetic impurities, 
e.g., oxygen. 

Price has shown that the K state would fit 
into the Rydberg series for his Class I system of 
electronic bands, those involving a '2—>'2 transi- 
tion. He doubts that this is anything more than 
a coincidence, inasmuch as the intensity of the 
bands he observed is at least 100,000 times as 
great as that of the Kistiakowsky bands, and 
the latter contain P, Q, and R branches whereas 
Price’s are double-branched. 

It is immaterial for the purpose of our dis- 
cussion just what the quantum designation of the 
K level is. The excellent work by Woo and his 
collaborators is unfortunately incomplete; it is 
thus impossible to check Jonesco’s 'Z designa- 
tion?* which, in view of the work of Woo, of 
Gopfert,?4 and of Kistiakowsky," seems improb- 
able. However, it is clear from the weakness of 
the Kistiakowsky bands that the transition 
1'¥+hv—K seems to be forbidden by some 
selection rule for it is unlikely that application 
of the Franck-Condon principle could by itself 
account for such an effect. 

In conclusion, we may emphasize that the 
choice of Herzberg’s 125.1-kcal. value for the 
heat of sublimation of graphite has led to the 
137.9-kcal. value for the C=C bond strength. 
This is the only value so far suggested for that 
bond strength which is consistent with the known 
energy levels of free CH. Thus the system of 
potential energy relationships here demonstrated 
lends further support to the 125.1-kcal. value for 
the heat of sublimation of graphite. 


23M. A. Jonesco, Comptes rendus 200, 817 (1935). 
4H. Gopfert, Zeits. f. wiss. Phot. 34, 156 (1935). 
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The Exchange of Oxygen Between NO and NO,* 


EpGaR LEIFERT 
Department of Chemistry, Columbia University, New York, New York 


(Received February 19, 1940) 


The reaction N“O+ N50, = N“%O+ NO, has been stud- 
ied in the gas phase at —35°C and at a total pressure of 
2.5 cm of Hg. Although the reaction rate constant could 
not be determined due to the speed of reaction a lower 
limit for its value is 10° cc moles second. The equi- 
librium constant for the exchange at 298.1°K calculated 


from spectroscopic data is 0.967 as compared with 0.96 
+0.02 which was observed. By combining the data ob- 
tained in this investigation with the equilibrium data for 
N20; which is assumed to be the intermediate in the 
exchange reaction an approximate value for the energy of 
activation for the formation of N2O; is determined. 





HE existence of N2O3 as a stable compound 

in the gas phase has been shown by 

measurements of the equilibrium constant for the 
reaction 


NO+NO2=N,.0; (1) 


by Verhoek and Daniels! and by a study of the 
absorption spectra of mixtures of NO and NO; by 
Melvin and Wulf.? Although its structure has 
never been determined, a study of the electronic 
states for NO and for NO» molecules makes it 
seem highly reasonable that the bond which is 
formed in their association is between oxygen of 
NO: and nitrogen of NO. This is most probable 
since the main resonating structure of NO is the 
one having an odd electron on the nitrogen atom, 
and that for NO» having an odd electron on 
either of the two oxygens. Proof for this structure 
can be obtained from a study of the oxygen 
exchange reaction between NO and NOs as 
traced by isotopic nitrogen. 

In this research an attempt was made to study 
the kinetics of the exchange of nitrogen between 
NO and NOs in the gas phase. Owing to the speed 
of this reaction, only a lower limit for the rate 
constant could be established. However, the 
exchange itself established some interesting 
points about N2O3, which will be reported here. 
The reaction studied was 


N¥“O+N*O2=N4O+N"Oo. (2) 


The equilibrium constant for this reaction has 





* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

t Fellow of the Charles A. Coffin Foundation, 1939-1940. 
ass and Daniels, J. Am. Chem. Soc. 53, 1250 


(1938) H. Melvin and O. R. Wulf, J. Chem. Phys. 3, 755 


been calculated from spectroscopic data for the 
ordinary compounds by the method of Urey and 
Greiff.* The molecular constants used in the 
calculation are given in Table I. The molecular 
constants for the rarer isotopic molecules have 
been calculated from the spectral data of the 
more abundant molecules. The equilibrium con- 
stant for reaction (2) is given by 


K= fy f ws o2)/fox"o2 favo). 


The ratio of the partition functions is given by 
the following formulae :* 


/ M, ; I; 
foxres/fo=(—) (-) 
M2 Is 
he 1 —e—hews/kT 
Xexp | -—tr-o» |(—_), 
2kT 1 —e-hew/kT 


/ M;, 3 A,B,C, 2 
Jno.) /fos 169) = (—) (——") 
M, AoB2Ce 





hc 1 —emhewai/kT 
Xexp | ——— i- ||, os : 
P| Te ” Lo ] 6] —g-hew,i/kT 
At 298.1°K Kyaic=0.967. 
TABLE I. 
wo(CM~!) Io/]h* 
N¥O4 1891.3 1.0369 
N40 1857.4 


@1 w2 


w3 a(half-angle) A2/A1 B2/B; C2/C1 


N¥O,° 1373 641 1615 60° 1.0044 1.0486 1.0000 
N®O, 1356 634 1581 60° 


*H. C. Urey and L. G. Greiff, J. Am. Chem. Soc. 57, 
321 (1935). 

* Subscript 1 refers to ordinary molecule and 2 to rarer 
isotopic molecule. « 

* Nielsen and Gordy, Phys. Rev. 56, 781 (1939). 

5 Sutherland and Penney, Nature 136, 146 (1935). 
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EXPERIMENTAL 


NO was prepared by the method of Johnston 
and Giauque.* NO, was prepared by the 
catalytic oxidation of ammonia (2.4 percent N), 
prepared by Thode and Urey,’ over a single layer 
of 80-mesh platinum gauze. The products of 
oxidation were caught in a trap surrounded by 
a Dewar flask containing liquid nitrogen. They 
were then distilled through one meter of a P2O; 
drying tube, oxidized with excess oxygen and 
then frozen to dry ice temperature. The excess of 
oxygen was then pumped off with a mercury 
diffusion pump. The NO» was kept at tempera- 
ture of dry ice at all times to prevent its attacking 
the ‘“‘Apiezon”’ stopcock grease which is attacked 
if allowed to remain in contact with NO, for any 
period of time. Since the reaction discussed is so 
very rapid, the ‘“‘Apiezon”’ grease used to lubri- 
cate the stopcocks of the reaction vessel and 
filling chamber was not attacked even after many 
runs. 

In the attempted kinetic runs measured 
amounts of NO were admitted from a filling 
chamber into a reaction vessel which already 
contained a known amount of NO: (heavy in 
N). The reaction vessel was immersed in a 
Dewar flask containing bath of chloroform at 
— 35°C. The total pressure in the reaction vessel 
was 2.5 cm. More than half of the pressure was 
always due to NO. After about 15 seconds the 
reaction vessel was opened to an evacuated 
system which contained a trap and sample tube 
in series. The trap, which was kept at —118°C 
with solid-liquid ethyl bromide, separated quan- 
titatively all N.O3 and the excess NO was frozen 
into the sample tube which was immersed in a 
Dewar flask containing liquid nitrogen. The NO 
was then analyzed for its isotopic ratio with a 
mass spectrometer of the Bleakney*® type. Since 
the reaction was found to be very rapid, it was 
necessary to design the reaction vessel so as to 
eliminate the error which might arise due to the 
rate of diffusion of gases. The rate of diffusion of 
gases in an ordinary bulb proved to be slower 
than the rate of reaction even at the low pressures 
in the experiments. The reaction vessel was a 


“_ and Giauque, J. Am. Chem. Soc. 51, 3194 
(1929). 

7 Thode and Urey, J. Chem. Phys. 7, 34 (1939). 

8 Bleakney, Phys. Rev. 40, 496 (1932). 
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125-ml bulb; through a ring seal in it a one- 
millimeter capillary tube extended about two- 
thirds of the way into the bulb. Attached to this 
tube was a three-way stopcock, one opening of 
which was to the bulb, another to the filling 
system, and the third to the analysis system. The 
jet effect produced by allowing NO from a higher 
pressure in the filling system into the reaction 
vessel, and the fact that NO is lighter than NO, 
already in the bulb caused the gases to mix 
rapidly so that diffusion was completed within 
15 seconds. 

The equilibrium constant for the reaction was 
determined by mixing larger amounts of NO and 
NQ.. Since the original isotopic ratios were 
known for both NO and NOs, a mass spectrome- 
ter analysis for the isotopic ratio of NO at 
equilibrium would give us the equilibrium 
constant. 

Experiments were also performed on solid 
N2O3. Solid NO; containing a high isotopic 
content of nitrogen was mixed with NO at 
—118°C using a solid-liquid ethyl bromide 
mixture as the cooling bath. The excess of NO 
was analyzed for its isotopic ratio after remaining 
in contact with N2O; for several minutes. 


RESULTS AND DISCUSSION 


The mere fact that we have obtained exchange 
proves that in its structure N2O; has an oxygen 
bridge between both nitrogens. This follows 
immediately from the assumption that the ex- 
change which was observed occurs through the 
formation of N2O3. Although the kinetics for this 
reaction could not be measured by the method 
used since all samples were at equilibrium within 
15 seconds, a lower limit for the rate constant ‘‘k” 
defined by 


d(N150) /dt = K ,k(N“O) (N"°Ox) 
—k(N¥9O)(N"O3), 


where K, is the equilibrium constant for reaction 
(2), is approximately 10° cc moles second at 
— 35°C and at a total pressure of 2.5 cm of Hg. 
Since the rate constant for this reaction is s0 
high even at so low a temperature, it must mean 
that the energy of activation for the formation of 
N.Os3 is very small and is most probably less than 
one kilocalorie. Since AH for reaction (1) is 
10,300 calories,! the energy of activation for the 
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dissociation of NO; into NO and NOs is ap- 
proximately 10.5 kilocalories. These results com- 
pare favorably with those for the association of 
NO: to form N2O,4*® where a similar situation 
exists. We should normally expect these results 
since it is probable that very little energy of 
activation is required to form a bond between 
two molecules, both of which have odd electrons. 

The equilibrium constant determined experi- 
mentally is 0.96+0.02 as compared with the 
value 0.967 at 25°C calculated from spectroscopic 


® Richards and Reid, J. Chem. Phys. 1, 114 (1933). 


data. This reaction might, therefore, also be used 
for concentrating isotopic nitrogen, although it 
seems less feasible from an experimental point of 
view than the method used by Thode and 
Urey.’ 

From the exchange which was found between 
NO and solid N2O3 at —118°C, it appears that 
solid N2O; is unstable and can exist only in 
equilibrium with its components. 

The author wishes to express his appreciation 
to Professor Harold C. Urey for his helpful 
advice throughout the course of this investigation. 
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The Molecular Structure of Isobutane 


J. Y. BEAcH AND JOHN WALTER 
Frick Chemical Laboratory, Princeton, New Jersey 


(Received January 15, 1940) 


The two published electron diffraction investigations of isobutane led to slightly conflicting 
results. We have repeated the calculation of the theoretical scattering curves using Z—F 
instead of Z for the atom form factor and including the temperature factor. The final results 
are: C—C distance, 1.54+0.02A; C—C—C angle 111° 30’+2°. 


WO electron diffraction investigations of 
isobutane have been published. Both led to 
a C—C bond distance of 1.54+0.02A. There was 
a slight disagreement, however, on the value of 
the C—C—C angle. Beach and Stevenson! re- 
ported 113° 30’+2°. Pauling and Brockway? 
reported 111° 30’+2°. In this paper we shall 
attempt to resolve this discrepancy. The posi- 
itons and intensities of the electron diffraction 
maxima and minima obtained in the two inves- 
tigations do not differ significantly. In this 
re-examination of the problem we shall use the 
average of the two sets of data. The averages of 
the observed values of s (sops=(4m sin 0/2)/d) 
and J (visually estimated intensities) are given 
in Table I. 
In both investigations the usual visual method® 
was used to interpret the photographs. In this 





1J. Y. Beach and D. P. Stevenson, J: Am. Chem. Soc. 
60, 475 (1938). 


- 2 li é “ ® . . . 
1223 (1937). © and L. O. Brockway, J. Am. Chem. Soc. 59, 


*L. Pauli 0. eae 
867 (1934), and L. O. Brockway, J. Chem. Phys. 2, 


procedure the visual estimates of the positions 
and intensities of the maxima and minima are 
compared with the positions of the maxima and 
minima on the simplified theoretical curve, 
I=CDY id ZZ; sin sri;/sr:i;. In this expression 
the atom form factors (Z;—F;)/s? have been 
replaced by Z;. This involves two simplifications. 
First, the neglect of 1/s? changes the curve from 














TABLE I.* 
Max. Min. I Ck Sops_ | Scatc* | Scatc/Sops 
1 10 16 3.00 2.77 | (0.923) 
2 4.36 4.10} (€.940) 
2 28 80 5.48 5.75} 1.049 
3 6.99 6.94} 0.993 
3 11 42 7.88 8.04} 1.020 
4 8.63 9.11} (1.056) 
4 (shelf) 8 33 9.31 9.60} 1.031 
5 11.38 | 11.13} (.978 
5 19 85 12.88 | 12.94) 1.005 
6 12 44 17.69 |17.79| 1.006 
7 6 17 21.12 | 20.83} 0.986 
Average 1.009 
C-C 1.55 























* Values in parentheses not included in the average. 
« Curve D, Fig, 1, 
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a rapidly falling one to an approximately hori- 
zontal one possessing a series of maxima and 
minima. This curve represents the visual ap- 
pearance of the photographs much more closely 
than the complete curve. Second, it is assumed 
that Z — Fis proportional to Z for different atoms, 
with the same proportionality factor, at a given 
angle. It is easy to check the accuracy of the 
first simplification. We can be sure that Z—F 
is proportional to Z by considering a molecule 
in which all the atoms are of the same kind. This 
is true of the halogen molecules. These molecules 
have been treated by the visual method* and 
interatomic distances are obtained which agree 
with the spectroscopic values. This is good proof 
that the first simplification introduces no serious 
error. The second simplification is not com- 
pletely justified,* although it introduces only a 
slight error for atoms in the same row in the 
periodic system. In order to eliminate a possible 
source of error in isobutane we have used Z—F 
in the scattering function. The Pauling and 
Sherman! F values were used in the calculations. 

In both of the investigations of isobutane the 
effect of the thermal motion of the molecule on 
the scattering curve was not considered. We 
have now included the temperature factor in the 
calculations. To include the temperature factor 
it is necessary to know the mean of the square of 





J. Y. BEACH AND J. 





WALTER 





the deviation of each interatomic distance in the 
molecule from the equilibrium value. To know 
this accurately requires a complete normal coor- 
dinate treatment of the molecule. The mean 
square displacement of each interatomic distance 
can then be calculated for each normal vibration. 
The total mean square displacement for an inter- 
atomic distance is the sum of the mean square 
displacements for all the normal vibrations.® 
Expressions have not been obtained for all the 
normal vibrations of isobutane. We have split 
the problem into two parts. We consider first the 
vibrations of a pyramidal XY; molecule. As 
normal coordinates we have used the symmetry 
coordinates of Howard and Wilson.’ Their sym- 
metry coordinates correspond quite closely to 
either stretching vibrations or bending vibra- 
tions. This is a fairly good approximation. We 
use this treatment to obtain the two C—C mean 
square displacements. The two shortest C—H 
and the shortest H—H mean square displace- 
ments we obtain from a similar treatment of the 
X YZ; molecule. 

This accounts for all the terms in the scattering 
formula that do not change with the orientation 
of the methyl groups. The scattering formula 
including Z—F and temperature factor but not 
including rotation terms, in the case of a tetra- 
hedral C—C—C angle, is: 


I=3(Zo— Fe)*(sin 1.54s/1.54s) exp (—0.0014s?) +3(Z¢— Fe)*(sin 2.52s/2.52s) exp (—0.0020s?) 
+10(Z¢— Fe)(Zu— Fu) (sin 1.09s/1.09s) exp (—0.0028s?) +12(Z¢— Fe)(Zu — Fu) 
X (sin 2.16s/2.16s) exp (—0.0034s?) +9(Zy — Fu)*(sin 1.78s/1.78s) exp (—0.0074s’). 


This curve is shown in Fig. 1, curve B. It is only 
slightly different from the curve calculated with- 
out Z—F and temperature factor, curve A. The 
introduction of Z—F and temperature factor 
does not have an important effect on the 
C—C-—C angle. As before! we find that the 
assumed orientation of the methyl groups has 
a considerable effect on this angle. Previously! 
the methyl groups were supposed to be sta- 
tionary in an orientation which allowed the 
H—H distances to be as large as_ possible. 
Pauling and Brockway omitted the rotation 
4M. H. Pirenne, J. Chem. Phys. 7, 144 (1939). 


5L. Pauling and J. Sherman, Zeits. f. Krist. A81, 1 
(1932). 








terms. To determine the most reasonable orien- 
tation of the methyl groups we have constructed 
a model of isobutane from a set of atom models.’ 
From an examination of the H—H distances we 
found that the most stable orientation giving the 
molecule the symmetry C3 (not C3y) is the one 
used before.! However, it is evident that the 
methyl groups can rotate to all positions more 
or less as cog wheels without having the H—H 


®R, W. James, Physik Zeits. 33, 737 (1932). 
asdiy” Howard and E. B. Wilson, J. Chem. Phys. 2, 630 
8 Fischer-Hirschfelder models, obtained from the Fischer 
Scientific Company. The radii in bond directions are 
obtained from observed interatomic distances. The radii in 
other directions are kinetic theory radii. 
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distances become seriously smaller. For any 
orientation of a given methyl group the other 
two can be oriented in such a way that the H—H 
distances remain reasonably large. We have 
therefore assumed, in calculating the long C—H 
terms, that the methyl groups rotate freely. 
Each methyl group was assumed to have an 
equal chance of occupying: one of the twelve 
orientations differing by 30° (or m 30°). The long 
C—H terms might have been accounted for by 
considering stationary methyl @roups and using 
a large ((61:;7)4 large) temperature factor. Free 














Fic. 1. Theoretical intensity curves for isobutane. Curve 
A:C—C—C angle, 109° 28’; Z-F, temperature factor and 
long C—H terms not included. Curve B: C—C—C angle, 
109° 28’; Z-F and temperature factor included. Curve C: 
C—C—C angle, 109° 28’; Z-F, temperature factor and long 
C—H terms included. Curve D: C—C—C angle, 111° 30’; 
Z-F, temperature factor and long C—H terms included. 
Curve E: C—C—C angle, 113° 30’; Z-F, temperature 
factor and long C—H terms included. The vertical lines 
indicate the observed positions of the maxima and minima. 
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Fic. 2. Modified radial distribution function for iso- 
butane. Vertical lines indicate the final values for the 
interatomic distances and their importance in the scatter- 
ing function. 


rotation is probably a much better representation 
of the motion of the molecule. Only the long 
H —H terms were not included in the calculation. 
The short H—H terms were included. 

Curves C, D and E were calculated using 
Z— F, temperature factor and long C—H terms. 
The C—C—C angles are 109° 28’, 111° 30’ and 
113° 20’, respectively. The fourth maximum on 
the photographs is reproduced most closely in 
the 111° 30’ curve. In curve C it is too low and 
in curve E it is too high. Quantitative comparison 
of curve D and the photographs is made in 
Table I. The radial distribution function cal- 
culated by the method of Schomaker® is shown 
in Fig. 2. The C—C bond distance is 1.55A in 
Table I and 1.54A in Fig. 2. The long C—C 
distance is 2.56A in Fig. 2 giving an angle of 
$47". 

The final results are: C—C distance, 1.54 


+0.02A; C—C—C angle 111° 30’+2°. 


®V. F. Schomaker, Baltimore meeting of the American 
Chemical Society, April, 1939. 
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Errors in the Determination of the Dissociation Constant of a Weak Acid by the 
Extrapolation Method 


Mary L., KILpatRick 
Department of Chemistry and Chemical Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 


(Received February 12, 1940) 


In the extrapolation method for determining the dissociation constant of a weak acid the 
value of the limiting equivalent conductance is not known from independent measurements, 
and both the limiting equivalent conductance Ao and the dissociation constant K are obtained 
from conductimetric measurements on solutions of the acid. While clearly convenient, the 
reliability of the method has been questioned. In the present papgr equations are derived 
which express the effect on K and on Ao of errors in the experimental quantities. For water as 
solvent, tables are given showing the numerical values of the errors. The possibility of detecting 


constant errors is discussed. 





ITH the realization that the relationship 

between the structure of an acid and its 
dissociation constant requires a knowledge of the 
thermodynamic quantities, free energy, heat con- 
tent, and entropy, the importance of the accurate 
determination of dissociation constants over a 
range of temperature becomes evident. The 
conductimetric determination of the thermo- 
dynamic dissociation constant of a weak acid 
HX commonly requires a knowledge of the 
equivalent conductance of the acid at infinite 
dilution, obtained by measurement of the con- 
ductances of the salts MX and MY and that of 
the strong acid HY. It has been proposed by 
Ives,! Fuoss and Kraus,” Fuoss,* and Shedlovsky ,* 
that the dissociation constant may be evaluated 
from measurements on the weak acid alone. 
This procedure, known as the extrapolation 
method, has however been criticized by Shed- 
lovsky and by Belcher.5 After comparing the 
values of the dissociation constant computed by 
the extrapolation method for a number of weak 
acids in water with the values obtained in the 
customary manner using the conductance of the 
alkali salt, Belcher concluded that the extrapola- 
tion method is reliable if the electrolyte is not 
too weak (K210-*), but that if the electrolyte 
is much weaker (K=10~-5) the method may give 
Ao and K values that are seriously in error. 


1D. J. G. Ives, J. Chem. Soc. 731 (1933). 

2 R. M. Fuoss and C. A. Kraus, J. Am. Chem. Soc. 55, 
476 (1933). 

3R. M. Fuoss, J. Am. Chem. Soc. 57, 488 (1935). 

* T, Shedlovsky, J. Frank. Inst. 225, 739 (1938). 

5 D. Belcher, J. Am. Chem. Soc. 60, 2744 (1938). 
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In solvents in which all electrolytes are weak 
the experimental determination of Ao is not 
possible, and the extrapolation method is per- 
force applied. Although Ao is experimentally 
determinable in water, the experimental determi- 
nation presents a laborious task if the conduc- 
tances of HY and MY are not already available. 
It therefore appears worth while to examine the 
problem of the accuracy with which the con- 
ductance of the weak acid HX must be measured 
in order that the extrapolation method may 
yield a reliable value of the dissociation constant. 

The extrapolation method is based upon the 
assumption of obedience to the equation of 
Onsager for the equivalent conductance of a 
completely dissociated solute,* and to the equa- 
tion of Debye and Hiickel for the activity coeffi- 
cient of an ion. Let A be the equivalent con- 
ductance of the weak acid HX present at total 
concentration C. The sum of the ion conduc- 
tances at this acid concentration, A., is given by 
the equation 


Ae=Ao—(LAo+M)(CA/A-)! 

; =Ao—(LAo+ M)Ci, (1) 
where L and M are constants for a given solvent 
and temperature, and C;=CA/A, represents the 
concentration of each ionic species. The dissocia- 
tion constant K is given by the equation 


Céf* A?Cf? 
Cal, idd~t) 





*In the development given by Shedlovsky, reference 4, 
an extended form of the Onsager equation is employed. 
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where f is the mean ionic activity coefficient, 
that of the undissociated molecules being taken 
as unity. Upon setting 


logio f= —AC#= —A(CA/A.)}, 


A being a constant for a given solvent and tem- 
perature, we have 


A2C- 10-24¢8 
K= 





or 


A+(LAo+M)(CA/A,)! 





1 FA2C-10-24(CAsAe)} 
| } @ 


.—-— 
A. 


K 
[A+(LAp+M)(CA/A,)*]=y 
is a linear function of 
[A2C- 10-24 (CA/d0)8/A =x, 


the slope of the line being —1/K and the inter- 
cept Ao. The solution of (1) and (2) for Ap and K 
is effected by a series of approximations from an 
initial approximate value of Ao.* The reader is 
referred to the papers cited for details of solution. 

Two points (Ai, C; and As, C2) furnish two 
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equations of the form of (1) and two of the form 
of (2); from these values of the four unknowns 
Ae, Aez, Ao, K may be obtained. For n experi- 
mental points there are n+2 unknowns, Ae, «+: 
Aen, Ao, K and 2n equations, m of the form of (1), 
n of the form of (2). The ” equations of the form 
of (2) may be written briefly as 


1 
yi=Ap—-—*1, 

K 

1 
Yn=Ao——Xn. 


For 1 points of equal weight the theory of least 
squares gives as the best value of K and of Ao 


n>x*?— (=x)? 
K =———__—_ (3) 


’ 
rxzy—n=xy 


LTxrlxy —Ix*Zy 
Ao = . (4) 
(Sx)? —nrx? 





Eqs. (3) and (4), together with the ” equations 
of the form of (1), furnish +2 equations. 
From these one may compute the effect upon K 
and Apo of error in the experimental quantities, 
concentration and specific conductance. 


THE PARTIAL DERIVATIVES 


Upon differentiating (1), there results for the kth point 











1—L(CyAp/Aex)? Aex(Ao—Acx)fdAx dC, 
dc. =2Acr lig cnneeanconiel hats | (5) 
3Aez— Ao 3Nea—Ao LA, C, 
Upon differentiating 
Ap2C 10724 (Ce Ak! Nex) 3 A,2C, 
x-= = 
, rm 1024 ((Ao—Acy)/(LA0+M)) Ay, 
and introducing the expression above for dAc, there results 
dx), = upd Ay +v,dC,+ wd Ao, (6) 





*It may be asked what is the effect upon the linearity of (2) when x and y are computed from an incorrect value of 
Ao. The effect is small, as the following figures show. Brockman and Kilpatrick (J. Am. Chem. Soc. 56, 1483 (1934)), 
found Ao for benzoic acid in water at 25° to be 382.10; the values of ordinate and abscissa in (2) computed therefrom are 
designated as Yoand Xo, respectively. Letting the values computed from Ao=1000 be y’ and x’, it is found that the maxi- 
mum deviation of y’—0.4 from yo, and of 2.563x’—0.00019 from xo, is less than 0.2 percent. Since x’ and y’ are obtain- 
able (within the accuracy stated) by linear transformation of x» and yo, respectively, y’ is an approximately linear 


function of x’. 
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where, letting 


- 2.303(2A) 
Jovemmmnmen, 
LAg+M 
han Ag— Ack , 1 
“es 2+————_(1 — jAex) 
Axl 3! er — Ao ll 


xf Ao— Ack : 
(1 — jAex) 


2 =—| 1+———. 
CiL 3Aezn— Ao 
1 (—— 


3Ae—Ao\ LAo+M 











Lion) —23] 


w= 


and similarly from 
ve=Anct (LAg +M)(CyAz/Ae)*= Ayt+Ao— Acx 


there results 
dy; =rydAp+sidC,+tidAo, 


1 Ao— Aex 
'.~= 1+—| Nex (—) 
Ak 3Aex— Ag 


1 Ao — Nex 
5, = | Nex (——~) 
C;. 3Ae.— Ao 


QNer(LAee+M) 
(3Aex—Ao)(LAo+M) 


where 








th=1 
From Eq. (3) it follows that 


n (OK 0K 
dK=> |—dn-+—dyu}, 
k=11 OX, OV: 


where 

0K 1 
—— =———— {2[nx,—Ix]+K[ny.—Xy]} 
Ox, Uxty—n=Zxy 


0K K{[nx,.— =x] 
OY; = rxLTy—nN=xy 


and from (4) it follows that 

n ONo OAo 

dAy= > —dx;, +—dy,, 

k=1 OX, OY: 

where 
OAo 1 

— = ———_—§_ | [2xyt yy, 2x —2x,.2y]+2Aclnx,— =x }} 
Ox, (2x)?—nz=x? 


OAg x,2x—Dx? 


Oy, (2x)?— ndx? 





(7) 


(8) 


(9) 


When the expressions for dx, and dy, from (6) and (7) are introduced into (9), there results 


1{ n 
seit’ D> mdAyt+ndC;y , 
k=1 


(10) 





(7) 


(8) 


(9) 


(10) 
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where 
n OAg OAg OAo OAo OAg OAg 
j=1-|5 Wy— +x ; M,=—Uyz+—?*r; 1 = Ve +—S. 
k=1 OX, OY OX; OV: Ox; OV: 


We are now in a position to express dK in terms of the dC;,’s, dA,’s. Upon introducing into (8) the 
expressions for dx; and dy, given in (6) and (7), with dAo replaced by its value from (10), there results 


n 0K W,| 
dK=> | nda tod. —| > mad naCy 


k=1 OX, j=1 


n OK ty. { 

+> | ndtutsdCe+-| ss mas +nacy| | (11a) 
k=1 OY;, j=1 

Since the specific conductance x, is the directly measured quantity, rather than the equivalent 

conductance A,;, dA; will be replaced by its value in terms of d& and dC,. 


A,= 





1000%, di, dC, 
’ =a | 
Kk Cy 


Ci 
Eq. (11a) then becomes 


n OK dk, dC, Wy { 2 dk; dC; 
dK =>" —| uss —-—) +dCi+—| ym ms(—- )+naci} | 
j=1 K 


k=1 OX; Ki k Kj j 





n OK dk; dC; te { 
+¥ —| rats t/t sdcit“ S ait (11b) 
k=1 OY, K Ll li=i 


Ki k 
from which it follows that 


——— mt Mgt tig ty 2 Oy hy} FO 


+] 


OK -| | OK OK m,n» OK }] QO; 
Ax = 














OK, =o OX: OV: Ll i=1 Ox; OY; Kk (12) 
'“@K 1 0K OK nn. » OK OAK 1 
—=—| forts. > w—+—| ~0.|=—[P.- 011. 
OC; C. OX: OY: l j=1 Ox; OY; Cy 
Upon introducing into (10) dA,;’s value in terms of dz, and dC;, there results 
1 n di; dC; 
ary] 5 mA + tm Cx—mds)— |, (13) 
LLk=1 Kk k 
from which it follows that ° 
OAgo myAyp 1 OAg [npCy. my,A, 71 
ens —.| i E (14) 
OR: lL k& OC; l l C. 


EVALUATION OF ERRORS FOR A WEAK AcipD _ is plotted against C=(CA/A,)', a linear relation- 

IN WATER ship obtains up to a certain concentration, after 
which log K. drops off. Thus for benzoic acid in 
water at 25° 





Before proceeding further it is necessary to 
discuss the range of applicability of the equations 


derived. It is known that when for a weak acid logiy Ke=logio K+1.011C;3 
of the type of HX 
A2C for solutions of ion concentration less than 
log —————-= log K. 2.9 10-4N. ® Since the slope 1.011 is close to the 


A.(A.—A) theoretical value of 2A for this solvent and 
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temperature, the conclusion is drawn that in the 
case of benzoic acid in water, Eq. (2) is valid for 
solutions sufficiently dilute so that C; does not 
exceed 2.9 10-4N. When the data for the acids 
listed by Belcher? in Table IV of his paper are 
similarly examined, it is seen that for carbonic 
(K =4.31X10-") the limiting ion concentration 
is not reached experimentally, i.e., (C;),>1.2 
X10-*; for propionic (K=1.343X10-5), (Ci): 
=3.7X10-‘; for u-butyric (K=1.508X10-5), 
2.5 10-4; for acetic (K = 1.753 X 10-5), 2.0K 10-4; 
for benzoic (K=6.312X10-5), 2.910-*; for 
o-chlorobenzoic (K=1.197X10-%), 1.5X10-%; 
and for monochloroacetic (K=1.396X10-%), 
2.7X10-*. For weak acids in water the range 
C;=4X10-> to C;=2X10~ is selected for ex- 
amination, that is to say, the assumption is made 
that over this range Eqs. (1) and (2) and the 
partial derivatives derived therefrom are valid. 
At these low concentrations of univalent ions it 
does not appear necessary to use Shedlovsky’s 
extended equation for A, in place of Eq. (1); 
according to Shedlovsky,‘ Onsager’s equation is 
valid up to a concentration of 0.001 or 0.002N for 
uni-univalent electrolytes in water. 

In Table I are given values of the partial 
derivatives for an acid for which in water at 
25°C, K=1X10-* and Ay=400. The constants 
L, M, and A have been taken as 0.2274, 59.79, 
and 0.5065, respectively. The concentrations 
have been chosen so as to give equal increments 
of x on the graph of Eq. (2), and the range of ion 
concentration is from 4X10-> to 2X10-‘N. 
In all instances concentration is expressed in 
molarity, specific conductance in ohm! cm-. 
The partial derivatives listed are computed from 


TABLE I. Values of the partial derivatives. K = 1X 107; 
Ao= 400; n=11; C;=4X10~ to 2X10. 











OK 0K .0Ao0 OAo 
aad ‘aap — me | 3G 2 
Cix106| C104 Kx<:108 | Kaz * 104} CGE XIOHKGE X10 | CoE X10" 
0.4000 | 0.5577| 1.596 | —2.423 | +1.890 | +3.426 | —2.672 
0.4860 | 0.7184| 1.939 | —1.985 | +1501 | +2'874 | —2.173 
0.5720 | 0.8935| 2.282 | —1.568 | +1.154 | +2348 | —1.728 
0.6720 | 1.115 | 2.680 | —1.112 | +0.797 | +1.772 | —1:270 
0.7860 | 1.391 | 3.133 | —0.628 | +0.438 | +1.160 | —0.809 
0.9130 | 1.728 | 3.639 | —0.131 | +0.092 | +0:531 | —0.372 
1.064 | 2.169 | 4.240 | +0.406 | —0.269 | —0.149 | +0.099 
1.239 | 2.735 | 4.935 | +0.964| —0.625 | —o.855 | +0.554 
1.442 | 3.464 | 5.742 | +1.536 | —0.973 | —1:580 | +1.001 
1.693 | 4.474 | 6.739 | +2.149 | —1.330 | —27358 | +1.459 
2:000 | 5.870 | 7.957 | +2.785 | —1.685 | —3.165 | +1.915 


























7D. Belcher, reference?5. References to the original 
papers are given by Belcher. 


formulas (12) and (14). Since the points lie on 
the line of Eq. (2), 


OK 1 0K OAo 1 OAo 


=—— and —=——. 


OX}, K d0yx Ox, K dy; 


The probable error in K resulting from random 
errors in the experimental quantities is 


AK= (x mae) of acs) ). 
ORE 


If the relative errors in the x’s and C’s are less 
than e, 


AK « 11 
—<( E (a) +(o--) ) <6.6¢«, 
K K =1 OK: 


and similarly 
AAo 





Ao 


Turning to the consideration of constant 
errors, we have for the case of a constant relative 
error in specific conductance 


.- 0K Ak 0K Ak 
[AK ]}= > —Ak.=— D &— = — 7X 107 — 
k=1 OR}, K OR, K 
AK | Ak 
[ ne wn 





K 


i.e., an error of 1 percent in all specific con- 
ductances causes an error of —0.007 percent in 
K. Error in the cell constant is therefore of little 
consequence in the determination of K by the 
extrapolation method. This point has been 
stressed by Ives.! For the case of a constant 
relative error in concentration 


— ma AC J 
AK |\c=>. —AC,=— 
Cc ac, k C "Ce 
AC 
= +0.990 X is 
[AK ]c AC 





ie., an error of 1 percent in all concentrations 
causes an error of practically 1 percent in K. 
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TABLE II. Relative errors in K, from measurements at two 
concentrations. 


TABLE III. Relative errors in Ao, from measurements at two 
concentrations. 








1.28]|1 1,98 
Kane} |KL ‘aca 





i aaX| 1 [c of | 
K K [ase KL~*ac: 


C1 X104| C2 X 108 

















ral a2] a oe [ose] t[c:dt 
K C1 1041C2103|Ao lL" dei J |Ao| ” dx2 | |Ao| ~'8C1 | |AoL "aCe 

















Cir =4X10, k1=1.596 X10; Ciz=2X10, k2=7.957 X1075. 


Ci =4X10-5, ks: =1.596 X10; Ci2z=2X10-4, k2=7.957 X 10-5, 





1X10-3 | 0.4158 | 0.2387 | —9.729 | +9.717 +9.365 —8.358 
1X10~4 | 0.5577 | 0.5870 | —5.622 | +5.615 +4.385 —3.397 
5 X10-5 | 0.7153 | 0.9740 | —7.245 | +7.236 +5.034 —4.051 
1X10-5 | 1.977 | 4.070° | —22.32 | +22.30 +12.44 —11.49 
1X10-6 | 16.17 | 38.90 —168.5 | +168.3 +85.48 —84.82 
1X1077 | 158.1 | 387.2 —692.4 | +691.6 +347.6 —347.8 


0-3 | 0.4158 | 0.2387 | +1.354 —0.354 — 1.303 +0.305 


1X1 

1X 10-4 | 0.5577 | 0.5870} +2.238 — 1.237 —1.745 +0.748 
5 X10-5 | 0.7153 | 0.9740} +3.216 —2.214 —2.235 + 1.239 
1X10-5 | 1.977 | 4.070 +10.92 —9.90 —6.085 +5.104 
1X10-6 | 16.17 | 38.90 +84.30 —83.20 —42.78 +41.94 
1X1077 | 158.1 | 387.2 +347.3 —345.9 —174.3 +174.0 





Ci =2X104, ki =7.957 X10-5; Ciz=1X10-3, K2=3.952 X10-4, 


Ci; =2 X10, ki =7.957 X10; Ciz=1X1073, k2=3.952 X1074, 





















































1107 | 2.039 | 1.093 | —16.86 | +16.82 | 416.51 | —15.45 1X10~ | 2.039 | 1.093 | 41.311 | —0.312 | —1.283 | +0.287 
510-8 | 2.077 | 1.186 | —10.09 | +10.06 | +9.707 | —8.691 510-2/ 2.077 | 1.186 | +1.359 | —0.360 | —1.307 | 40.311 
1X10-3 | 2.387 | 1.929 | —5.357 | +5.343 | +4.608 | —3.628 1X10-8 | 2.387 | 1.929 | +1.743 | -0.744 | —1.500 | +0.505 
Similarly, In contrast to neutral salt, a small amount of 


Ak 
[Ado ]k= +4.004 x 10?—, 


K 
[Av Ak 
and chan ae a 
Ao i 


AC 
[AAo]e= —3.996X 10, 


. [AA AC 
and tAede .00—; 
Ao 


i.e., 1 percent errors in all x’s and in all C’s cause 
errors of +1 percent and —1 percent, respect- 
ively, in Ao. For the case of a constant absolute 
error in specific conductance, such as would be 
caused by the presence of a small amount of 
neutral salt, 


OK 
[AK ]x avs = Ak >> —= —26.53Ak, 
OK. 
[AK ]x abs ‘ 
————- = — 2.65 X 10° Ak; 


for Ak=4X10-8 ohm cm, the error in K is 
—1.1 percent. The corresponding error in Apo is 


aA 
[AAo]k ave = Az © — = +4.604X107 Ai, 
Ki 


[AAo |x abs 
————— = $1.15 X 105 Az; 





Ao 


for Ag =4X 10-8 ohm-! cm-!, the error in Ao is 
+0.5 percent. 


basic impurity causes a diminution in specific 
conductance and an increase in the calculated 
value of the dissociation constant. For example, 
if it is assumed that the water used as solvent 
contains NaHCO; and H2COs, and has a specific 
conductance of 210-7 ohm! cm- and a 
hydrogen-ion concentration of 4X10-7N, the 
Ai’s for the solutions listed in Table I range from 
—5.5X10-* for the most dilute to —3.910-% 
ohm! cm~ for the most concentrated solution of 
acid. The resultant errors in K and Ao are 


AK 1 0K 
—=— )) —Ak,,= +1.55 percent, 
x -« OR; 


—=— } —Ak,. = —0.64 percent. 
Ao Ao OR} 





COMPARISON OF ERRORS FOR ACIDS OF 
DIFFERENT STRENGTHS 


It can be shown that small constant errors in 
the experimental quantities have little effect 
upon the linearity of (2). That being the case, it 
appears justifiable to approximate the errors in 
K and Ao by regarding K and A» as determined by 
two points, which correspond to the extremes of 
the range of ion concentration. For n=2 the 
computation is considerably shorter than in the 
previous instance where ” was taken as 11. The 
values of the partial derivatives obtained from 
(12) and (14) for n=2 are listed in Tables II and 
III. The solvent is water, the temperature 25°, 
and Ao is 400. 








312 MARY L. 





KILPATRICK 











4020 


4100 4018 


Log Ke 


4092 


1010 



































O6 OBuixlO* Je 149 oY 


Fic. 1. Curve A, down 0.00200. 
Curve - - -, down 0.00300. 


For n=2 the errors in K and Ao due to a 
constant absolute error in specific conductance 
are, for an acid of dissociation constant 110-4, 


[AK x abs 
——— = — 2.82 X 10° Ak, 
[AAo |x abs 
and —— a ae, 
0 


as compared with —2.65X105Ak and +1.15 
X10°Ak, respectively, for »=11. The errors re- 
sulting from a constant relative error in specific 
conductance or in concentration are the same for 
n=2 as for n=11. The effects of constant errors 
in the experimental quantities are in general so 
little different for n=2 and n=11 that the 
approximation appears warranted. 

It has already been mentioned that for 
o-chlorobenzoic and monochloroacetic acids, 
which in water have dissociation constants of the 
order of 1X10-*, log K. is a linear function of C;? 
up to ion concentrations somewhat above 
1X10-*N. The lowest ion concentrations studied 
in these investigations were 5X10-‘N for o- 
chlorobenzoic® and 1X10-4N for monochloroac- 
etic acid.® It seemed worth while, therefore, to 
include in Tables II and III errors computed by 
use of the ion concentrations 2X10-* and 
1X10-*N in the case of the stronger acids. 

8B. Saxton and H. F. Meier, J. Am. Chem. Soc. 56, 
1918 (1934). 


®T. Shedlovsky, A. S. Brown, and D. A. MacInnes, 
Trans. Electrochem. Soc. 66, 165 (1934). 


is... 1 1 1 1 
08 Ver x /0e = /2 14 06 


Fic. 2. Curve B, up 0.00500. 
Curve - - -, up 0.00500. 


i i i 
08 7 x/0e le i4 


Fic. 3. Curve C, up 0.00900. 
Curve - - -, up 0.00350. 


The effect upon K of a constant relative error 
in specific conductance is seen to be minor. For 
acids of dissociation constant between 1X10-% 
and 1X10-* a 1 percent error in k causes in K an 
error not exceeding 0.02 percent, and for the 
weakest acid listed (K=1X10-") a 1 percent 
error in k causes in K an error of 0.8 percent. The 
error in K caused by a constant error of 1 percent 
in concentration drops from 1 percent for the 
strongest to —0.2 percent for the weakest acid. 
The effect upon K of a constant absolute error in 
specific conductance increases markedly with 
diminution in K for acids weaker than 5X10~-°. 
As regards Ao, constant relative errors of 1 
percent in x and C cause errors of approximately 
1 percent in Ao, except in the case of the weakest 
acid (K=1X10-’) where the effects are 1.4 
percent and 0.3 percent, respectively. The error 
caused in Ap by a constant absolute error in 
specific conductance increases with decrease in K. 

The most striking conc usion to be drawn from 
Table II concerns the random errors. For acids of 
dissociation constant between 2 or 3X10-* and 
510-5, the individual errors in K are of 
approximately the same magnitude for a given 
percentage error in % or C. For acids in this 
range, therefore, the extrapolation method should 
yield about equally reliable values of the dis- 
sociation constant provided the experimental 
accuracy remains the same. The high values of 
the individual errors given in Table II for acids of 
dissociation constant 1X10-* and 110-7 show 
that the determination of K by the extrapolation 
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method is not practicable in the case of acids 
which in water are very weak. 


DETECTION OF CONSTANT ERRORS 


As already mentioned, small constant errors in 
the experimental quantities have little effect 
upon the linearity of (2). They have a greater 
effect upon the linearity of the relationship 
between log K. and C;}. In the accompanying 
figures log K, is plotted against C; for an acid of 
dissociation constant 1 10-‘, for the cases where 
(Fig. 1, curves A and A’) there is neutral salt 
present in the water, where (Fig. 2, curves B and 
B’) the water contains basic impurity, and where 
(Fig. 3, curves C and C’) the concentrations 
employed in the computation are 1 percent too 
high. 

Data for the figures are taken from Table I; in 
A and A’, 4.0 10-* ohm—! cm is added to each 
specific conductance; in B and B’, from each 
specific conductance is subtracted an amount 
varying from 5.5 X10-* ohm cm-! for the most 
dilute solution to 3.9 10-8 for the most concen- 
trated, the impurity in the water used as solvent 
being regarded as NaHCO; and H2COs; of such 
concentration as to give the water a specific 
conductance of 2X10-? ohm-! cm™ and a 
hydrogen-ion concentration of 4X10-7N; in C 
and C’ the recorded concentrations are increased 
by 1 percent. The curves designated by the 
unprimed letters are constructed by use of 
Ao=400 and A,.=400—(400L+ M)(1000%/A.)! 
= 400.00 —150.75C;3. Those designated by the 
primed letters are constructed by use of Ao 
obtained by the extrapolation method, i.e., from 
the partial derivatives of Table I we have 


for y o- AAo= +1.84, Ao = 401.84, 

A.=401.84 — 151.17(1000%/A.)'; 
for B’, AAg=—2.56, Ao=397.44, 

A. = 397.44 —150.17(1000%/A.)!; 
for C’, AAgp=—4.00, Ao=396.00, 

A. = 396.00 — 149.84(1000k/A.)?. 


These cases are selected because they cause in K 
as given by the extrapolation method errors of 
approximately the same magnitude; from the 
partial derivatives of Table I we have 


for A’, AK=—1.06X10-*; 
for B’, AK=+1.55X10-; 
for C’, AK=+0.99X10~. 


The broken lines in the figures are drawn with 
intercept 4.0000 and slope 1.013, ie., they 
represent the case of noerrors. Continuous curves 
have not been drawn through the points of 
A’, B’, C’ since they would obscure the lines 
drawn to show departure from linearity. It was 
not considered worth while to include a figure for 
the case of a constant relative error in specific 
conductance, since an error of 1 percent in all x’s 
causes in K as given by the extrapolation method 
an error of less than 0.01 percent for an acid of 
this strength. 

The curves A, B, C show large deviation from 
linearity, and if lines are drawn through the 
points of higher ion concentration, their slopes 
differ from the theoretical. Deviation from 
linearity is less pronounced for the curves desig- 
nated by the primed letters. A’ and B’ neverthe- 
less show definite departure from linearity, A’ 
being concave upward toward lower ion con- 
centrations and concave downward toward 
higher, and B’ being concave downward. The 
curve C’ shows only slight departure from 
linearity; it is slightly concave upward. The 
slope of the line drawn through the points of 
higher ion concentration is greater than the 
theoretical for the curve A’, and less than the 
theoretical for the curve B’. From curves A’ and 
B’ it appears that, if consistent data for the acid 
are at hand, it should be possible to detect the 
presence of neutral salt or basic impurity in the 
water, even though the plot of Eq. (2) exhibits a 
linear relationship. The detection of a constant 
relative error -in concentration appears more 
difficult; the detection of a constant relative 
error in specific conductance is comparatively 
unimportant when the extrapolation method is 
employed to determine the dissociation constant 
of a weak acid in water. 
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A Note on the Efficiency of the Thermal Diffusion Process for Separating Isotopes 


J. W. WESTHAVER* AND A, KEITH BREWER 
Bureau of Agricultural Chemistry and Engineering, U. S. Department of Agriculture, Washington, D. C. 


(Received February 1, 1940) 


The mechanism proposed by Gillespie for the thermal separation of gases is correlated with 
the attendant conduction of heat; this is done by describing gaseous heat conduction as a heat 
exchange between two oppositely bombarding Maxwellian groups which differ in temperature 
at each isothermal surface by a finite AT. The resulting formulas enable a computation of AT 
the bombardment ‘‘equipartition distance ” and the transport yield. An ideal system of cells 
each having temperature limits of 300°K and 500°K is found to require a minimum of 880 kw 
hours of conducted heat per gram CH, extracted to a concentration of 20 percent. 





HE transport of the components of a 

gaseous mixture by action of a tempera- 
ture gradient is receiving special attention as a 
method for segregating isotopes. Until the recent 
paper by Gillespie! no adequate attempt had 
been made to explain this action on the basis of 
elementary considerations.? The success of Gil- 
lespie’s theory, particularly its agreement with 
observed equilibria separations, appears sufficient 
to warrant an application of its underlying 
mechanism to a study of other factors of the 
transport process. 

The mechanism assumes that the difference dT 
in the average temperatures of contiguous thin 
gas layers normal to the temperature gradient 
can be regarded as the difference in the effective 
temperatures of the molecules which bombard 
the opposite sides of the plane separating the 
two layers. If these layers are of finite thickness 
the bombardment on the two sides of the plane 
results in a net flow of molecules in the direction 
of increasing temperature. To preserve the con- 
dition of uniform pressure this net flow of mole- 
cules must be exactly compensated in number by 
an opposite superimposed mass return flow. This 
action occurs in any stagnant gas when subjected 
to a temperature gradient. In the case of homo- 
geneous mixtures of gases of unequal molecular 
weight, the numbers of molecules of a particular 
component flowing.in the two directions are not 
equal ; for this reason a transport occurs. 

The following discussion relates particularly to 
the case of a vertical temperature gradient 
through a cell of gas in which the top wall is 


* Now at the U. S. Patent Office. 
1 Gillespie, J. Chem. Phys. 7, 530 (1939). 
* Furry, Jones and Onsager, Phys. Rev. 55, 1083 (1939). 


heated and the bottom wall is cooled; the con- 
fined gas is necessarily stagnant. The molecules 
in each element of volume can be classified into 
two groups one having components of velocity 
toward and the other away from the heated wall. 
In virtue of the fact that heat is conducted across 
the plane there must exist a finite difference in 
the heat energies carried by the two groups. 
Assuming that each group is composed of the 
same kind of molecules it follows that the 
groups differ in temperature at the plane by a 
finite amount, designated as AT. 

It can be shown that the heat transport due to 
the net bombardment and corresponding mass 
return as above defined is proportional to the 
square of AT and hence can be neglected. In 
consequence, AT is found by dividing the rate 
of heat conduction kdT/dx across the plane by 
Ycipi(M;/27RT)' where c; is the specific heat at 
constant pressure and ;(M;/2rRT)? is Lang- 
muir’s*® formula for molecular bombardment of a 
surface, expressed in g cm~ sec.—! when ; is the 
partial pressure in dynes cm~*, M; the molecular 
weight when O= 16, and R=8.315 X10’ ergs per 
mole—°K. The summation is over-all compo- 
nents of the mixture. As an example of its magni- 
tude, AT is found to have the value 7.1 X 10-* °K 
in methane at 40 cm Hg where T=400°K, 
c=0.60 cal. g- deg, R=1.0X10-* cal. cm™ 
sec.—! deg.—! and dT’ /dx=200°K cm-. 

If dx is the finite thickness of the thin layers 
defined above, then dT is related to AT by 


kdT/dx 


T= . (1) 
Lcip:(M;/2xRT)} 





’ Langmuir, Phys. Rev. 2, 329 (1914). 
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Also if dx is the mean distance along the gradient 
over which the two bombarding groups can 
move before they partition their energy, the 
medial planes of the contiguous layers may be 
regarded as the emission sources for the groups, 
and, since these sources have the average tem- 
perature of their respective layers, AT becomes 
identical with dT. The equipartition distance dx 
is then given by the expression 


k 
dx= . (2) 
Yeip:(M;/22RT)' 





When the same constants are used as in the 
previous example, dx has the value 3.6 X 10-5 cm. 
This is about $ the mean free path for collision 
in methane under the conditions stated; dx and 
dT are therefore finite quantities of reasonable 
value. 

After an expression for AT has been obtained, 
Langmuir’s formula can be used to find the net 
bombardment for any component g whose re- 
spective group temperatures at the plane are 
T+AT/2 and T—AT/2. By subtracting from the 
mole net bombardment of g the mole mass 
return of g the initial rate of transport of com- 
ponent g across the plane of temperature T is 
found. This computation gives for the initial 
transport of light component 1 toward the hot 
wall in g cm-? sec.“ 


dQ M.zpi/M3 
it the A ws 
dt 2Trc:p:M;3 Zp: 


where dQ/dt=kdT/dx is the heat passing be- 
tween the boundaries in cal. cm sec.~! as 
conventionally computed. For a two-component 
system : 

dQ PipeM13(M2} — M;}) 


"dt 2T(pi+-b2) Mak(cipiMy-+-cop>M!) * 





(4) 


which is dependent on relative but not’ total 
pressures. When the molecular weights and 
specific heats are nearly the same, Rp becomes a 
maximum when p= po. 

Since the rate of transport across the mid- 
section decreases with time in the manner 
R,=Roe~*, the corresponding change in partial 
pressure across the cell for either component is 
Ap=Ap.(1—e-**) where Ap, is the equilibrium 
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change. Maximum rate of separation occurs 
when R,Ap is maximum. This is at the time 
when e~**=3; the maximum yield is therefore 
R,o/2dQ/dt and the small change in partial 
pressure is Ap,/2. 

The separation factor is conveniently defined 
as the ratio of the ratios of partial pressures at 
the hot and cold walls; for maximum rate of 
separation it may be written 


Ap.(pit pe) 
4+-— —---— 
2pipe 


where Ap, can be computed from Eq. (7) given 
by Gillespie. Since Ap, is small, his equation can 
be rewritten 


pips(M2!— M}') 
* 2(p1Mi+p2M2!) 


j= 


(5) 





In 72/7}. (6) 


For such isotopes as heavy carbon and 
nitrogen, pz will be small and J, nearly equal 
to M2, so that the denominator of (6) is nearly 
2(pitp2)Mi' and (5) becomes 


M2}—M;}} 
Sy= 1+———— In T?,/T3. (7) 
4M;} 


By assuming a mixture 0.01C%H,+0.99C"H, 
at 400°K average temperature, where ¢,=Ce2 
=0.60, the maximum yield 3Ro(dQ/dt)— is found 
from Eq. (4) to be 3.110-7 gram CH, toward 
the hot wall or 3.3 10-7 gram CH, toward the 
cold wall per calorie conducted heat. The corre- 
sponding energy requirement is 3.5 kw hours 
per gram C!*H, transported across the cell. For 
wall temperatures of 300°K and 500°K the cell 
separation factor S; from Eq. (7) is 1.0040. This 
haif-equilibrium condition is reached in a time 
which varies as the square of the wall separation. 
In the present case ¢; is about 0.15 second at 40 
cm Hg and wall spacing of 1.0 cm. 

As an example of accumulative yield, assume 
that N cells of equal height, heated at the top 
and cooled at the bottom to the same tempera- 
tures, are connected in series so that the hot gas 
layer of one cell and the cold gas layer of the 
following cell undergo a continuous interchange 
of gas and thereby acquire the same composition. 
If the hot gas layer of the first cell mixes with a 
constant source of 0.01C%H,+0.99C"H,, a 































a eee 


heavy mixture may be continuously bled from 
the cold layer of the final cell. Since Eq. (4) shows 
that the rate of transport for unit of transverse 
area varies as pif, the transverse area of the rth 
cell must be a fraction pip2/(pipe), of the first 
cell area; this is necessary to insure a total rate 
which does not vary along the series of cells. 
If 20 percent C!*H, is desired an over-all separa- 
tion factor of 25 is required. For this case 
(1.0040)* = 25 and 810 cells are required. 

The total heat conducted across the gas is 
S Pipe(pip2); dr times that of the first cell. 
The integration is made by noting that 


Pipo( Pipe), * = (pet piSy") (pit p2S)’) 


where (/1+/2),=1, the cell separation factor S, 
being the same for all cells. The result is [27pipe 
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— (prSy-’ — p2S,’) (In Si) Jo” which equals 250 
for the present example. The maximum yield of 
the system is therefore 1/250 of the first cell, 
namely 3.310-7/250=1.3X10-® gram CH, 
per calorie for 20 percent concentration. The 
corresponding energy requirement is 880 kw hr. 
per gram C!*H,. This neglects all heat transfer 
other than by conduction through the gas. 

These calculations have been carried through 
primarily to show the applicability of the 
mechanism suggested by Gillespie for the thermal 
separation of gases. However, the theoretical 
yield obtained suggests the feasibility of a system 
in which the operating principle is that of a 
series of cells. 

The writers wish to thank Dr. Arthur Bramley 
for his interest and constructive criticisms. 
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The Specific Heat of Tantalum at Low Temperatures and the Effect of Small Amounts 
of Dissolved Hydrogen* 


K. K. KELLEyt 
Metallurgical Fundamentals Section, Metallurgical Division, Bureau of Mines, United States Department of the 
Interior, Berkeley, California 


(Received January 8, 1940) 


Specific heat determinations of pure tantalum were made 
in the temperature range 53° to 298°K and the entropy 
was found to be Sos.1=9.9+0.1. Tantalum containing 
hydrogen in amounts ranging from Ta+0.0284 H to 
Ta+0.0958 H was shown to behave anomalously, the 
specific heat having a high maximum in the region 220° 
to 265°K depending on the composition. The height of the 
maximum and the difference in measureable heat content 
between hydrogenated and normal tantalum were shown 
to be approximately proportional to the amount of dis- 
solved hydrogen. The temperature of the maximum in- 


N pursuance of one of the programs of study 
of the Pacific Experiment Station of the 
United States Bureau of Mines, low-temperature 
specific-heat measurements of a sample of elec- 
trolytic tantalum were made. This material, 
labeled ‘‘sample F’’ below, was found to have 
* Published by permission of the Director, Bureau of 
Mines, U. S.. Department of the Interior. (Not subject to 
copyright.) 


+t Chemist, Metallurgical Division, Bureau of Mines, 
U. S. Department of the Interior. 





creases as the hydrogen is increased. After dehydrogena- 
tion, the same specific heat was obtained as for the original 
unhydrogenated material. The density of the dehydro- 
genated product also was normal, whereas the hydro- 
genated materials had lower densities. It is postulated that 
differences in entropy between hydrogenated and pure 
tantalum may be attributed to differences in the number 
of possible ways of arranging the hydrogen particles among 
the positions in the tantalum lattice available at the 
higher and lower temperatures. 


an unexpected, marked “hump” in its specific- 
heat curve centered near 275°K. The sample was 
known to contain about 6.4 atomic percent of 
hydrogen, and moreover it is probable that it 
did not consist of a single crystalline species. The 
desirability of obtaining a more adequate entropy 
figure for tantalum and the interest attached to 
the “hump” itself were the incentives for study- 
ing the effects of relatively small amounts of 
hydrogen on the low-temperature specific heats 
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of tantalum with which this paper is concerned. 
Six samples were investigated, the electrolytic 
material already mentioned, three containing 
different amounts of hydrogen, one of wrought 
tantalum wire containing virtually no hydrogen, 
and one made by pumping the hydrogen from 
material previously hydrogenated to 8.74 atomic 
percent. 

The scope of this investigation has been re- 
stricted intentionally to the range 0 to 12 atomic 
percent hydrogen, as Hagg' has shown by x-ray 
studies that it is near the latter composition that 
the structure begins to change from a _ body- 
centered cubic lattice to hexagonal close packing, 
the change being complete near 33 atomic 
percent. Hagg has stated that in the range below 
12 atomic percent hydrogen the introduction of 
hydrogen merely causes an increase in the 
lattice parameter without change of structure. 
This swelling also has been noted by Pietsch and 
Lehl? and Chaudron, Portevin, and Moreau.® 


MATERIALS 


The tantalum was loaned for this work by the 
Fansteel Metallurgical Corporation through the 
courtesy of Mr. J. H. Harper and Mr. F. L. 
Hunter. 

Sample A consisted of pure wrought tantalum 
in the form of 0.006’ diameter wire. The wire 
was sheared into lengths averaging about 0.5 
cm, washed consecutively with HCl solution, 
NaOH solution, distilled water, and alcohol, and 
then dried near 100°C. 

Sample B was made from sample A by heating 
to 710°C in vacuum (10-' mm Hg at room 
temperature to 2X10-* mm Hg at 710°C) and 
then admitting hydrogen to give the atomic 
ratio 0.0676 : 1. The hydrogen was purified by 


TABLE I. Density measurements. 











SAMPLE DENSITY (d4°%5°) 
A,noH 16.557 
E, 0.0284 H 16.552 
B, 0.0676 H 16.436 
C, 0.0958 H 16.330 
D,H pumped off C 16.614 











’ Hagg, Zeits. f. physik. Chemie B11, 433 (1931). 
* Pietsch and Lehl, Kolloid Zeits. 68, 226 (1934). 


* Chaudron, Portevin, and Moreau, Comptes rendus 207, 
235 (1938). 


SPECIFIC HEAT OF TANTALUM 
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Fic. 1. Specific heats: A, D—pure Ta; E—Ta+0.0284 H; 


B—Ta+0.0676 H; C—Ta+0.0958 H. 


passing through activated charcoal immersed in 
liquid nitrogen and was measured by filling a 
calibrated flask sealed to the system through a 
stopcock. About half the hydrogen was absorbed 
in a few seconds at 710°, and the remainder was 
taken up as the tantalum was cooled slowly to 
300°. 

Sample C was made from sample B by evacu- 
ating to 10-> mm Hg at room temperature, heat- 
ing to 700°C, and admitting additional hydrogen 
to raise the atomic ratio to 0.0958 : 1. 

Sample D was prepared by pumping the 
hydrogen from sample C at 720°C. The pumping 
was continued until the pressure had fallen to 
and remained static at 2X10-* mm Hg at 720° 
and, of course, during the cooling to room tem- 
perature. 

Sample E was made by rehydrogenating 
sample D at 700°C to the atomic ratio 0.0284 : 1. 

Sample F was the electrolytic powder men- 
tioned in the introduction. According to the 
makers, this material contained 99.85 percent 
Ta, 0.11 percent C, 0.005 percent Fe, 0.02 percent 
residual impurities including oxide, and about 
75 cc hydrogen, measured at atmospheric 
pressure and room temperature, per cc of metal. 
The last figure corresponds to 0.037 weight 
percent or 6.4 atomic percent hydrogen. 


DENsITY MEASUREMENTS 


The densities of all samples except the elec- 
trolytic powder were determined with a pyc- 












TABLE II. Specific heats. 
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SAMPLE A D E B Cc F 

T °K Cp Cp Cp Cp Cp Cp 
53 2.81 2.71 
55 2.95 2.84 
60 3.26 3.14 
65 3.54 3.44 
70 3.80 3.72 
75 4.03 3.96 
80 4.22 4,21 4.19 4.16 4.13 4.16 
85 440 | 4.38 4.37 4.34 4.32 4.33 
90 4.55 | 4.53 4.53 4.49 4.48 4.48 
95 4.68 | 4.66 4.66 4.63 4.62 4.62 
100 | 4.78 | 4.77 4.77 4.76 4.74 4.75 
105 4.87 | 4.86 4.88 4.88 4.86 4.87 
110 4.97 | 4.96 4.98 4.98 4.97 4,98 
115 5.06 | 5.05 5.08 5.08 5.07 5.08 
120 §.14 | 5.13 5.17 5.17 S17 5.18 
125 5.21 5.20 5.27 5.26 .| 5.26 5.27 
130 Bsr i S27 5.36 5.35 5.36 5.36 
135 5.33 | 5.33 5.44 5.44 5.44 5.44 
140 5.39 | 5.38 5.53 5.52 S00 5.52 
145 5.44 | 5.42 5.61 5.59 5.60 5.59 
150 5.48 | 5.46 5.69 5.66 5.68 5.66 
155 §.52 5.50 5.76 5.74 5.76 5.74 
160 So0 | 5.53 5.84 5.81 5.84 5.82 
165 5.58 | 5.57 5.91 5.89 5.92 5.90 
170 5.62 | 5.61 5.99 5.97 6.00 5.98 
175 5.65 5.65 6.07 6.05 6.08 6.06 
180 5.68 5.67 6.15 6.14 6.16 6.13 
185 5.70 | 5.70 6.24 6.23 6.25 6.21 
190 $.73 | 5.72 6.33 6.32 6.35 6.30 
195 5.76 | 5.74 6.43 6.42 6.46 6.39 
200 5.78 | 5.76 6.54 6.52 6.57 6.48 
205 5.80 | 5.78 6.65 6.64 6.70 6.57 
210 5.82 | 5.80 6.78 6.76 6.84 6.67 
215 5.84 | 5.82 6.91 6.91 6.99 6.76 
220 5.86 | 5.84 7.04 7.08 7.16 6.86 
225 5.88 5.86 6.74 7.27 7.36 6.96 
230 5.89 | 5.88 6.35 7.48 7.58 7.04 
235 5.91 5.90 6.03 7.73 7.81 7.13 
240 5.92 | 5.91 5.97 8.00 8.08 7.22 
245 5.94 | 5.92 5.98 8.29 8.38 7.30 
250 5.95 | 5.94 5.98 8.58 8.72 7.39 
255 5.96 | 5.95 5.99 8.76 9.10 7.48 
260 5.98 | 5.96 5.99 8.01 9.53 7.56 
265 §.99 | 5.97 6.00 7.23 9.80 7.64 
270 6.01 5.98 6.00 6.48 8.98 7.72 
275 6.02 6.00 6.01 6.21 8.07 7.78 
280 6.03 6.02 6.02 6.17 7.16 7.70 
285 6.03 6.03 6.02 6.14 6.33 7.59 
290 6.04 6.04 6.04 6.12 6.28 7.48 
295 6.04 | 6.05 6.06 6.10 6.25 7.37 
300 7.26 
310 7.12 
320 7.10 

















nometer, using distilled water as the fluid. 
Corrections were made for buoyancy and the 
temperature coefficient of the volume of the 
pycnometer flask. The measurements were made 
at 25+1°C. The results are contained in Table I. 

These measurements show that as hydrogen is 
introduced in these rather low concentrations the 
' tantalum swells, the rate of swelling increasing 
with the hydrogen content. If Hagg’s' lattice 








parameter for tantalum containing 12 atomic 
percent hydrogen (hydrogen-tantalum ratio, 
0.136 : 1), the composition at which the structure 
begins to change, is correct, then the density of 
such material is 15.82 and the swelling must 
become much more pronounced as this composi- 
tion is approached. 

Sample D, which was made by pumping the 
hydrogen from sample C, has about the density 
given by the best direct measurements in the 
literature* (ranging from 16.55 to 16.64) and by 
Hagg’s! lattice parameter, 16.64. Sample A, 
which also contains virtually no hydrogen, has 
a lower density which may be attributed to the 
effects of cold working. These effects apparently 
are entirely removed by hydrogenation followed 
by dehydrogenation at 700°C. Consequently, the 
density of sample D is more significant than that 
of sample A for comparison with the other results 
in Table I. 


SPECIFIC-HEAT MEASUREMENTS 


The apparatus and methods are very similar 
to those described in such detail elsewhere® that 
repetition is unnecessary here. The temperature 
scale employed is that of the laboratory of Pro- 
fessor Giauque.? The over-all error in the 
measurements is considered to be 0.5 percent or 
less. 

The results, expressed in defined calories (1 
calorie = 4.1833 int. joules) are given in Table II 
and shown graphically in Fig. 1. The figures in 
Table II were read from smooth curves while the 
plotted values are actual measured results. This 
procedure was followed to save space. All results 
are for 180.88 grams of tantalum plus contained 
hydrogen when present. 

A number of items of interest are immediately 
evident on examination of Fig. 1. First, sample 
A, which is the original unhydrogenated material, 


* Landolt-Bérnstein, Physikalisch-chemische Tabellen (Ju- 
lius Springer, Berlin, 1923), Vol. 1, p. 290; ibid. second 
suppl. (1931), Vol. 1, p. 212; ibid. third suppl. (1935), 
Vol. 1, p. 285. 

5 For example, see Maier, Trans. Am. Inst. Min. Met. 
Eng., Inst. of Metals Div. 122, 121 (1936). 

6a. Gibson and Giauque, J. Am. Chem. Soc. 45, 93 
(1923); b. Giauque and Wiebe, ibid. 50, 101 (1928); c. 
(1930)" ibid. 50, 1875 (1928); d. Anderson, ibid. 52, 2296 

7a. Giauque, Buffington, and Schulze, J. Am. Chem. 
Soc. 49, 2343 (1927); b. Giauque, Johnston, and Kelley, 
ibid. 49, 2367 (1927). 
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and sample D, which was prepared by pumping 
the hydrogen from sample C at 720°C, show the 
same specific heat within limits too small to be 
evident in Fig. 1 (average deviation observed 
<0.3 percent). From curves B, C, and E it is 
seen that the specific-heat maxima are of a higher 
order of magnitude than the quantity of hydrogen 
present. The ratio of the height of the maximum 
above the normal curve to the gram atoms of 
hydrogen present ranges from 39.7 to 41.3 and 
averages 40.7 for the three cases. Consequently, 
the virtual coincidence of the specific heats of A 
and D may be taken as evidence that the hy- 
drogen was nearly completely removed in pre- 
paring sample D, as something of the order of 
10-* gram-atom of hydrogen per gram-atom of 
Ta should have resulted in a detectable effect. 
This removal of hydrogen is also indicated by the 
density measurements; and, in addition, it may 
be stated that rough pressure observations made 
with both increasing and decreasing tempera- 
tures in the range 400° to 700°C on samples B 
and C indicate that the hydrogenation process is 
approximately reversible. 

Curves B, C, and E also show that both the 
temperature and the height of the maximum are 
increased by increasing the hydrogen content. 
Moreover, except for second-order differences, 
too small to be shown in Fig. 1, these three 
curves coincide below the maximum in £, and 
Band C coincide below the maximum in B. This 
indicates that the same quantity of change, 
whatever it may be, has taken place when a 
given temperature is reached independently of 
the amount of hydrogen present provided the 
temperature is below that of the lowest specific 
heat maximum. If the effect produced is at- 
tributed solely to the hydrogen, then the quan- 
tity changed from the low temperature to the 
high-temperature condition is independent of the 
total amount present, or of the amount un- 
changed, until the maximum in the specific-heat 
curve is reached. 

At the lowest temperatures studied all the 
samples are shown in Fig. 1 to have the same 
specific heat. Actually this is not the case and 
small, but regular differences are evident on 
examination of Table II, showing that increasing 
hydrogen content is accompanied by slightly 
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decreasing specific heats at temperatures below 
100°K. 

Curve F, for the electrolytic powder, cannot 
be considered a member of this group, although 
it very nearly coincides with the others at inter- 
mediate temperatures. It is believed that this 
sample did not consist of a single crystalline 
species, which is not at all unusual for electrolytic 
deposits of metals, and that there is present the 
resultant of the effects of hydrogen on the two 
(or more) crystalline species, which broadens the 
maximum and lowers its height. It was not con- 
sidered advisable to attempt to stabilize this 
sample or to remove the hydrogen therefrom, as 
it was a special material of considerable value 
and its owners stated that removal of hydrogen 
would result in a material very active toward 
oxidation by air or moisture. 

It was found possible to supercool the hydro- 
genated samples in the regions of the maxima. 
On rapidly cooling sample B from room tem- 
perature to 230°K the following results were 
obtained: 6.42 at 231.5°, 6.70 at 236.5°, and 7.01 
at 241.3°. Similar results were obtained with 
sample F; namely, 6.34 at 222.7°, 6.78 at 232.3°, 
7.17 at 242.3°, 7.39 at 252.5°, 7.60 at 262.6°, 
7.70 at 270.5°, 7.75 at 274.6°, 7.67 at 282.0°, and 
7.48 at 289.6°. Also the value 6.86 at 298.3° was 
obtained with sample F after rapidly cooling 
from 320°K to room temperature. However, the 
temperature extent of supercooling was limited, 
and it was found that quite reproducible results 
always were obtained after cooling to liquid 
nitrogen temperatures. In the case of sample B, 
two sets of determinations that agreed within 0.1 
percent were made in the range 80°-90°K after 
cooling to 78° under virtually identical conditions. 
For this same sample duplicate determinations 
also were obtained in the range 215 to 273° after 
warming from 85°K to 215°K at rates differing 
by at least a factor of 3. These agreed within 0.3 
percent. For sample F, duplicate determinations 
were made throughout the range 78° to 290°K 
after cooling from room temperature to 78°K 
under very different conditions. The first cooling 
was made in 2.5 hours, using helium for conduc- 
tion, while the second was made in vacuum (ca. 
10-*— 10-* mm Hg), the process taking 50 hours. 
Agreement within 0.2 percent was obtained 
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except near the maximum, 265° to 275°K, where 
the greatest deviation was 0.7 percent. 

The time required to reach equilibrium after 
energy input in the calorimeter for samples B, 
C, and E was 6 min. at the lower temperatures 
but increased to 15 min. at the temperatures of 
the maxima. Just above the temperature of. the 
maximum this time interval decreased rapidly in 
each instance, returning to 6 min. In the case of 
sample F, the sluggishness of energy distribution 
was more marked in the neighborhood of the 
maximum and did not decrease above this 
region. In fact, for sample F, the uncertainty in 
obtaining equilibrium above 280°K was so great 
that the inaccuracies in the specific-heat meas- 
urements were several times as large as those 
normally present, as is apparent on examination 
of Fig. 1. It is to be emphasized, however, that 
this last remark does not apply to any of the 
other samples. 

There is but little earlier work with which to 
compare the present measurements on pure 
tantalum. Good agreement is obtained with the 
value 6.024 calories per degree at 0°C that Jaeger 
and Veenstra® have reported as the extrapolation 
of their measurements in the range 300 to 1600°C. 
Simon and Ruhemann® have made six deter- 
minations in the range 71° to 79°K, employing a 
special apparatus that involved dropping samples 
from temperatures obtained with boiling oxygen 
and nitrogen to those obtained by boiling these 
liquids under reduced pressures. Their results are 
nearly 8 percent lower than those reported here. 
The reason for this difference is not apparent, as 
measurements made in this laboratory on 
chromium agree with the work of these inves- 
tigators while measurements on antimony differ, 
as in the present instance. 


TABLE III. Entropy of pure Ta. 








Ss3.1 (extrapolation) =1.40 
S298.1— 553.1 (graphical) = 8.54 


Sos. 1=9.9+0.1 








8 Jaeger and Veenstra, Proc. Acad. Sci. Amsterdam 37, 
61 (1934). 

®Simon and Ruhemann, Zeits. f. physik. Chemie 129, 
334 (1927). ; 
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TABLE IV. Entropy calculations for hydrogenated Ta. 








S298.1 —So 
MATERIAL (DEBYE) |S29s.1—So] 5S298.1 | CASE 1 | CASE 2 


Ta+0.0284 H | 10.17 10.18 | 0.24 | 0.26 0.27 
Ta+0.0676H | 10.43 10.45 | 0.50 | 0.49 0.52 
Ta+0.0958 H | 10.62 10.65 | 0.69 | 0.63 0.67 


























ENTROPY CALCULATIONS 


The entropy of tantalum was calculated from 
the results for samples A and D in the usual 
manner. Extrapolation below 53°K was made by 
means of the Debye function with 6=224, which 
adequately represents the measurements in the 
lower temperature range. The results are shown 
in Table III. The error has been taken as +0.1 
to allow safely for a possibly large error in the 
extrapolated portion. The value of So9s.1 is 0.5 
unit higher than that previously in use.” 

It cannot be stated that the entropies of the 
hydrogenated materials are zero at 0°K. How- 
ever, calculations of S2e9s.1.—So9 may be made in 
the same manner, account being taken, of course, 
of the small decrease in specific heat at the lower 
temperatures accompanying hydrogenation. In 
Table IV, values are given in columns 2 and 3 
based on two methods of extrapolation, the first 
involving Debye functions and the second the 
procedure of Kelley, Parks, and Huffman." The 
two methods yield substantially the same results, 
and those of the first method are adopted. The 
differences between So93.1—So for the hydro- 
genated samples and So9s.; for the pure metal are 
given in column 4 and designated 6S295.1. 


DIFFERENCES IN MEASURED HEAT CONTENTS 


Again using Debye extrapolations, the heat 
involved in raising the temperature of these sub- 
stances from 0° to 298.1°K, /o7C,dT, may be 
computed. The results are given in column 2 of 
Table V. Differences between the values for the 
hydrogenated samples and pure tantalum, desig- 
nated 6Ho3.1, are given in column 3, but it is 
not meant to imply that these substances have 
the same heat content at 0°K. Values in the last 
column give 1890, 2010, and 2000 calories, 


1 K. K. Kelley, Bureau of Mines Bull. 394, 1936, p. 43. 
( we Parks, and Huffman, J. Phys. Chem. 33, 1802 
1 ’ 
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respectively, when computed to the bases of 1 
gram-atom of hydrogen, the mean being 1970. 


DISCUSSION 


It has been accepted practice for some time to 
consider hydrogen dissolved in metals as dis- 
sociated, because in many instances it has been 
found that the amount dissolved is proportional 
to the square root of the pressure, at least over 
certain temperature and pressure ranges. From 
the recent statistical mechanical treatments of 
the solution process by Fowler and Smithells!” 
and Lacher," it appears very highly probable 
that the dissolved hydrogen is largely dissociated 
into protons and electrons, the latter adding to 
the metal electrons already present. Conse- 
quently, in considering the change in any 
property, after hydrogen of the amounts under 
consideration is absorbed by tantalum, three 
factors are involved—the presence of protons, 
the increased number of metal electrons, and the 
swelling of the lattice without crystal structure 
change. When, as in the present instance, only 
some of the resultant effects are known, it may 
seem absurd to attempt to reach any definite 
conclusions concerning the relative importance of 
these three factors. However, the contribution to 
the specific heat made by electrons is ordinarily 
taken to be very small in the temperature range 
under consideration. Also, in spite of differences 
(although not uniform) known to occur in lattice 
parameters, Eucken and Werth," for copper and 
nickel, and Maier and Anderson," for copper and 
aluminum, have found no difference in specific 
heats between cold-worked and annealed metal 
in this temperature range. Consequently, it 
appears probable that the anomalous behavior 
observed here may be attributed primarily to the 
presence of protons. 

If it is assumed that at temperatures far below 
the maximum in the specific-heat curve a single 
fixed stable position or set of fixed stable posi- 
tions is available to each proton while at tem- 
peratures above the maximum a greater number 
of positions is available, then the entropy dif- 


$7 (4 ae and Smithells, Proc. Roy. Soc. (London) A160, 
@Lacher, Proc. Roy. Soc. (London) A161, 525 (1937). 
*Eucken and Werth, Zeits. f. anorg. allgem. Chemie 

188, 152 (1930). 

** Maier and Anderson, J. Chem. Phys. 2, 513 (1934). 
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ference involved should be calculable from the 
number of possible ways of arranging the protons 
among the available positions in the two states. 
The expression for this entropy difference, deriv- 
able from probability considerations or by 
analogy with the ordinary entropy of mixing ex- 
pression for gases or ideal solutions, is 


asain (TF) 
-m()\(S)"] 


in which m is the number of gram-atoms of 
hydrogen, R is the gas constant, p is the number 
of positions available to each proton at the 
higher temperatures, and gq is the number of 
positions available to each proton at the lower 
temperatures. The number of positions available 
for either state is, of course, not determined by 
the data presented here. Two somewhat extreme 
cases will be considered, and the calculation will 
be restricted to whole numbers for p and g. Thus, 
if there is one position per tantalum atom avail- 
able at the higher temperatures and if it is 
assumed that at the lower temperatures the 
protons arrange themselves in some stable fashion 
that limits each proton to a single fixed position, 
the calculations, given as case 1 in Table IV, are 
in agreement with the measured results. In case 
2, 8 spaces per tantalum atom are considered 


TABLE V. Calculation of f,?C,dT. 











MATERIAL So CpdT 5H 298.1 
Ta 1362.7 
Ta+0.0284 H 1416.4 $3.7 
Ta+0.0676 H 1498.9 136.2 
Ta+0.0958 H 1553.9 192.2 








available at the higher temperatures; but it must 
then be assumed, to obtain agreement, that as 
the lower temperatures each proton has a choice 
of 3 positions. The agreement obtained cannot be 
considered as proof of this point of view, but it 
may be noteworthy that the trend in the meas- 
ured values with increasing hydrogen content as 
well as their magnitudes have been represented. 

For any two compositions considered, as men- 
tioned earlier, the extent of the transition reac- 
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Fic. 2. Plot of —log @ against 1/T: closed circle—Ta 
+0.0284 H; small open circle—Ta+0.0676 H; large open 
circle—Ta+0.0958 H. 


tion at a given temperature is the same regardless 
of the concentration of hydrogen, for all tem- 
peratures below that of the lower specific-heat 
maxima. This phenomenon may be considered 
somewhat analogous to heating different small 
amounts of pure liquid in a closed tube of such 
size that the liquid volume is always negligible. 
At a given temperature for all cases the vapor 
pressure would be constant as long as any 
appreciable amount of liquid were present, and 
the amount vaporized would be virtually con- 
stant. Such an analogy, however, does not por- 
tray the slowness with which the specific heats 
drop with the temperature above the maxima. 
The observation under consideration makes it 
reasonable to expect, however, that all data for 
temperatures below those of the maxima may be 
represented by a straight line in a log a vs. 1/T 
plot, where a is the quantity of protons changed 
to the higher temperature condition. 

At any temperature, a may be computed ap- 
proximately as the ratio of the observed 
Jo7AC,dT for a given sample to the average 
S?AC,dT=1970 calories computed for 1 
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gram atom of hydrogen present, the integrals 
being the areas between the hydrogenated 
specific heat curve and that of pure tantalum. 
The result is shown in Fig. 2. The slope of the 
line as drawn is too great compared with the 
1970 calories average excess heat-content change 
previously mentioned, but not unreasonably so 
considering the possible errors in the a’s, 
especially at the lower temperatures where AC, 
is a small fraction of the measured specific heats. 


SUMMARY 


Specific-heat measurements of pure tantalum 
have been made in the range 52° to 298.1°K, 
and the entropy of tantalum has been computed 
as Soo93.1:=9.9+0.1. 

Specific-heat measurements of tantalum con- 
taining, respectively, 0.0284, 0.0676, and 0.0958 
gram-atom of hydrogen have been made in the 
range 78° to 298°K. These substances have 
marked maxima in their specific-heat curves, the 
average excess specific heat at the maxima being 
40.7 calories per degree per gram-atom of dis- 
solved hydrogen. The excess specific heats at 
temperatures below all maxima concerned were 
shown to be independent of the amounts of 
hydrogen present, but both the temperatures and 
heights of the maxima are increased as the 
quantity of dissolved hydrogen is increased. 

On dehydrogenation the same specific heat 
was obtained as for unhydrogenated tantalum, 
and the density of the dehydrogenated product 
was normal. 

The entropy differences between the hydro- 
genated materials and pure tantalum perhaps 
may be attributed to differences in the number 
of positions in the tantalum lattice available to 
the hydrogen at the lower and higher tem- 
peratures. 
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Physics of Stressed Solids 
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The internal energy of a system is subdivided into a 
work or potential function and a thermal or kinetic func- 
tion, the former expressed in terms of the current electro- 
static theory of intercrystalline bonding, and these 
functions then examined for variations of temperature, 
hydrostatic pressure, unidirectional stress and combined 
hydrostatic and unidirectional pressure. From these con- 
siderations a theory is evolved which not only seems 
satisfactorily to explain and correlate phenomena of defor- 
mation, creep or plastic flow, cold working, elastic after- 
working, rupture, shear and certain other phenomena 
hitherto described as ‘“‘anomalous” effects but has been 
corroborated experimentally in some of its predictions, in 
particular for the effect of hydrostatic pressure on deforma- 
tion and compressive strength. The mechanism evolved 
consists of two processes—one an elastic deformation 
which is a function of the strain or potential energy of the 
system. Failure occurs here by “‘brittle’’ rupture wherein 


the maximum extension or maximum internal tension is 
the criterion. The other is a deformation by means of a 
two-phase transfer mechanism and is a function of the 
thermodynamic potential relations of the system. This 
latter type is also a function of time and therefore a func- 
tion of the rate of application of load. When both processes 
of this mechanism are operative failure occurs by shear; 
the criterion for this type of failure is given by a function 
of time, the strain or potential energy and the thermo- 
dynamic potential relations of the system. Expressions are 
derived for creep or plastic flow of polycrystalline sub- 
stances from the thermodynamic potential relations which 
not only satisfy the well-known phenomena of creep 
in metals but also express recent empirical creep data of 
some substances immersed in liquids in which they are 
somewhat soluble. An expression is also derived for the 
“brittle’’ potential type of rupture under combined thrust 
and hydrostatic pressure. 





HE thesis presented in this paper is that the 

phenomena of deformation, flow or creep, 

and rupture may be interpreted and correlated 

from a study of the internal energy stored up in 
the lattice as a result of deformation. 

The energy of deformation may be subdivided 
into energy of work and energy of heat, that is, 
into an energy of position and an energy of 
motion. If this energy of deformation exceeds a 
certain critical limit, determined by the physical 
characteristics of the material and the amount of 
externally applied compensation, a release of the 
strain energy takes place in such a manner that 
the system does the least work. If the energy 
relations are such that release is effected through 
the potential component, i.e., work function, the 
system exhibits potential or brittle rupture and 
the system relieves itself of its strain energy by 
doing work in the direction of greatest extension 
Which is the direction of least external con- 
straint. If the stress and energy relations are 
such that release may be effected through the 
kinetic component, i.e., the work done on the 
substance is largely dissipated by heat transfer, 
the specimen deforms by means of a “two- 
Phase” flow mechanism. This is the dominant 


effect in phenomena of creep and plastic de- 
formation. 

In general, deformation and failure are effected 
by a combination of these two mechanisms and 
observed as shear and as “gliding along shear 
planes.” 

Many of the ideas incorporated in this dis- 
cussion have been tentatively put forth by 
other writers but the picture as a whole with its 
correlations and the theory of rupture presented 
appears to be new and is therefore presented 
with the hope that it will prove of interest. 
This subject needs stimulus for further work on 
deformation phenomena under hydrostatic con- 
fining pressure because it is believed by the 
writer, perhaps also by others, that further 
advances in our knowledge of crystal lattice 
forces will come from such studies. 


EMPIRICAL HISTORY 


Volume compressibility, like density or heat 
capacity, is insensitive to variations in structure. 
By that we mean the volume compressibility of 
a single crystal is the same as that of a poly- 
crystalline aggregate of the same material. Shear 
strain, creep, and other related phenomena on 
the other hand, will vary with crystallographic 
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Fic. 1. Force vs, extension. 


orientation, the geometry of grain boundaries, 
and other structural factors. 

Under certain conditions the material may 
react like a ‘‘perfectly elastic body’’ whereas 
under other conditions it may react like a 
“viscous liquid.’’ In general our material will be 
found to behave in some intermediate manner 
for which several words have been coined such 
as ‘‘elastico-viscous,”’ ‘‘firmo-viscous,’’ and even 
“elastic flow.” 

Some other phenomena which must also be cor- 
related are: Strain hardening (work hardening) 
an effect for which yield stress increases with 
strain; and elastic afterworking. The stress is 
also found to vary with rate of deformation; 
for liquids the stress-flow rate slope is a constant 
related to the coefficient of viscosity. The type 
of deformation may also vary with the speed of 
loading and a ‘ductile’? or ‘malleable’ sub- 
stance may, under certain conditions of rapid 
loading, behave as a “‘brittle’’ substance. Again 
“aging,” which is presumably a slow transition 
toward a more stable state, is observed as a 
change of yield stress with time. A similar 
phenomenon may be observed in glass which 
although breaking readily along a fresh scratch 
does this with more difficulty after a lapse of 
time and eventually will no longer break cleanly 
along the scratch. 


THE ENERGY FUNCTIONS 


The current electrical theory of intercrystalline 
bonding forces is based essentially on an electro- 
static model. The electrostatic potential ¢ at a 


distance r will include the sum of charge poten- 
tials proportional to r~', of dipole potentials 
proportional to r~*, of quadrupoles proportional 
to r~*, of octupoles proportional to r~‘, and so on. 
There are also additional terms arising from 
interactions, induction and dispersion effects. 
For example the character of the van der Waals 
forces has been correlated with polarization 
forces (always attractive) produced by quadru- 
poles on molecules which are regarded as de- 
formable distributions of charge; the potential 
from this effect is proportional to r~*. These 
forces have also been correlated with the polariza- 
tion of one molecule by the time varying dipole 
moment of another; these interaction potentials 
have been computed as proportional to r~®, with 
also terms of higher order as r~* and r-". 

From quantum-mechanical considerations the 
repulsive potential is given by the encroachment 
energy from overlapping wave functions of the 
atoms and expressible, to fairly high pressures, 
by an exponential! term. The Thomas-Fermi 
atom model has been suggested and applied? to 
calculations of densities at very high pressures, 
i.e., exceeding 10° bars (1 bar=10° dynes cm~'). 

The force field is given by the negative gradient 
of the electrostatic potential. Our present pur- 
pose will be satisfied by noting that we can write 
the attractive force F on a particular element in 
the form 


F=> Cri-> C;r—, (1 
‘ 7 


where the C; and C; are, in general, complicated 
functions of space and direction. The derivative 
of F with respect to 7/7ro will denote the reciprocal 
of the linear compressibility, assuming linear 
symmetry. Volume cannot be arbitrarily sub- 
stituted for r* because of asymmetry with respect 
to the space coordinates. 

The coefficients C; and C; cannot, in general, 
be directly determined or computed; we may, 
however, expect to learn something of their 
respective importance from studies of deforma- 
tion under high hydrostatic pressure. Initial 
compressibilities have been computed for some 


1P. M. Morse, Phys. Rev. 34, 57-64 (1929). —_ 
2 J. C. Slater and H. M. Krutter, Phys. Rev. 47, 559-568 
(1935). H. Jensen, Zeits. f. Physik 111, 373-385 (1938). 
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ionic cubic crystals. Approximate correlations‘ 
of structure with certain physical properties such 
as viscosity have been obtained from an assumed 
two term expression. Although these simplified 
arbitrary expressions cannot be other than rough 
approximations they are able to extend the 
range of elastic theory® and also give us a very 
useful qualitative picture. Some explicit account 
may be taken of thermal energy pressure by 
adding a third term.* From such an assumed 
relation it is then possible to calculate the 
relative stability of certain types of lattices as 
functions of r. 

Schematic plots of these two functions F’7 and 
@ are given in Figs. 1 and 2 for an arbitrarily 
assumed two-termed expression wherein 7=3 and 
j=9. The horizontal scales are the same but the 
vertical scales differ. 

The internal energy may be divided into two 
portions, a work function and a heat function, 
and expressed as 


e=¥+Tn, (2) 


where ¢ denotes the internal energy, y the 
“maximum work” function, T the absolute tem- 
perature and 7 the entropy. Differentiation of 
this expression gives us 


(d€/dr) p= —2r+T7(07/0T),, (3) 


where 7 is a compressional stress (a tensional 
stress will be expressed as a negative pressure). 
The integral 


f rdr=(r) —¥( ©) = (1) (4) 


is the portion used herein as the potential 
function, and e=y~=¢ for T=0. 

The curve of Fig. 2 will also denote y as a 
function of r if the zero ordinate is replaced by 
the value ¥(«) where dy(«)/dT=—k, k denot- 
ing Boltzmann’s constant. 





*M. Born and J. E. Mayer, Zeits. f. Physik 75, 1-18 
(1932). A. May, Phys. Rev. 52, 339 (1937). 

‘References in R . Fowler, Statistical Mechanics 
(Cambridge University Press, 1936), Chapter 10. 

°F. D. Murnaghan, Am. J. Math. 59, 235-260 (1937), 
has used an arbitrary function similar in form to Ar*— Br~* 
for expressing pressure volume relations. 

* J. Bardeen, J. Chem. Phys. 6, 372-378 (1938). 

’ The term —F will hereafter in this article be spoken of 
as internal pressure and F as internal tension. The terms 
Pressure and tension without qualifying terms will refer to 
the external forces. 


Similarly the integral of Tdy from T=0 to 
T=T may be considered the thermal or kinetic 
function. A harmonic oscillator at temperature T 
has an average thermal energy of kT according 
to classical statistical mechanics. In summing up 
these energies for a system of oscillators, how- 
ever, the mutual coupling terms are neglected. 
Some calculations made for metals using Debye 
functions wherein the Debye temperature de- 
pends on the elastic constants have indicated the 
usefulness of this method but further develop- 
ment is needed. 

The mean internal energy curve will have a 
shape similar to the potential curve but with a 
minimum which lies at a higher level and is 
displaced in the direction of smaller r. Conse- 
quently the internal energy will at first decrease 
on application of pressure at constant tempera- 
ture by loss of thermal energy exceeding that 
gained from the work function y. Although in a 
large system the fluctuations in energy of the 
oscillators tend to cancel out, across any surface 
a certain fraction of them will have larger 
energies than the mean and, under favorable 
circumstances, can migrate from their positions. 

It might be anticipated that F,,, the maximum 
internal tension, of the F—r plot should bear 
some relation to the cohesive tensile strength 
and fm—?o, where rm is the separation distance 
for F=F,, and ro the distance for F=0, to the 
maximum elastic strain. This must however be 
of the nature of an upper limit because of 
(a) localized high stress regions resulting from 
inequalities of load distribution and from dis- 
tortions in the lattice structures caused by 
imperfections resulting from holes, fissures and 
inclusions of foreign elements, and (b) localized 
high energy levels resulting from thermal energy 
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Fic. 2. “Strain energy function” vs. extension. 
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gradients and fluctuations. The calculated tensile 
strengths, F,,, are all much higher than the 
observed values which are indicated roughly by 
So in the plot. In the following discussion F,, is 
used as the ideal elastic tensile strength and 
’¥m—?o as the ideal maximum safe elastic exten- 
sion. The actual values may be only about a 
tenth or less of these maxima. 

It should be noted that these considerations 
are made for the case where the “two-phase”’ 
transfer mechanism does not operate. Care must 
therefore be exercised in drawing comparisons 
between theoretical and empirical results and 
this will be shown more clearly later. 


TEMPERATURE 


If the temperature of the system is increased, 
at constant pressure thermal energy is absorbed 
with increase in the internal pressure. The 
system expands doing work against the external 
surroundings until a new positional equilibrium 
is attained for the new F. In this way we obtain 
a family of curves for different T. Ignoring for 
the time being the change in slope of these 
curves and considering only the large scale 
differences, we note that the effective zero 
ordinate of the F—r plot moves up and 7 
increases to some larger (70) 7. Similarly the mean 
minimum potential (¢o)7 moves to a larger 
(ro)r and higher up the right limb of the potential 
energy well which is therefore becoming shallower 
and shallower in the direction of increasing r. 
Eventually the right wall of this well would 
vanish and no further readjustments yield a 
stable configuration ; the system actually breaks 
down long before this point is reached. If in this 
breakdown a reorganization in potential and 
kinetic functions of the components can take 
place for a new minimum (¢o)7 the system melts, 
and if such a redistribution cannot be effected 
the system sublimes. The size and nature of the 
elemental units in this breakdown will depend 
on the variations in the interatomic bonds, the 
structure splitting up across bonds of least 
cohesive strength. 

In this connection it might be of interest to 
note that if the thermal energy be assumed as 
proportional to v~!, or to r~, the internal pressure 
contribution from this source will be propor- 


tional to r~*. Under increasing hydrostatic pres- 
sure the other internal pressure terms of higher 
power, here assumed as 7~°, eventually dominate 
over the r~* term and at an increasing rate. The 
thermal expansion should therefore decrease 
under hydrostatic pressure and approach zero at 
very high pressures. This conclusion, in con- 
junction with the third law, means that, for this 
model, the entropy must be zero at infinitely 
high pressure for all temperatures.* At moderate 
pressures the smaller powers of 1/r will play 
more effective roles and the initial rate of de- 
crease of dilatation should, as for compressi- 
bility, be less for simple ionic lattices since here 
the attractive r~ term is initially more effective 
than the r~* term. 


HYDROSTATIC PRESSURE 


The same type of reasoning applies to systems 
under hydrostatic pressure at constant tempera- 
ture but in this case the increase in potential is 
equal to the work done on the system. The 
system decreases in volume until an equilibrium 
between internal and external pressure has been 
achieved. The most stable pressure or phase 
configuration is the one best able to withstand 
the external pressure and therefore the one of 
higher density. A change of phase under pressure 
may take place either by reorganization of form 
(potential or work energy) or in part by change of 
form and in part by change of momenta (con- 
version to thermal energy). The pressure at 
which the latter type of phase change occurs 
will be a function of temperature; the former 
should be independent of temperature. 

Silicate glasses should in general behave some- 
what abnormally because these glasses are in a 
metastable energy state, i.e., the configuration 
stable for the high temperature kinetic and 
potential energy relationships has been frozen 
into the system. The internal forces should 
therefore be tensions but the configuration is 
presumably prevented from collapsing by the 
randomness of orientation inherited as a result 
of the previous thermal kinetic energy, that is 
to say, the glass is not only metastable from 
considerations of heterogeneous or phase equi- 


8 This was also suggested by G. N. Lewis, Zeits. f. 
physik. Chemie, Cohen Fest. Band 130, 532-538 (1927), 
from other considerations. 
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librium but is also metastable from considera- 
tions of homogeneous equilibrium. 

This is suggested as the reason why viscosity 
bears a relation to the previous thermal history 
and why anomalous heat capacity effects are 
observed in the annealing range. 

Resistance to motion and reorientation is a 
function of viscosity so we should expect 
anomalous effects in highly viscous substances 
quenched from high temperatures. In fact we 
might expect a possible increase of compressi- 
bility with pressure since the system is in a 
state which is up on the right limb of the energy 
curve plotted in Fig. 2 and should be under 
internal tension. 

Bridgman® observes anomalous effects in the 
compressibilities of quartz glass and of basalt 
glass but he suggests a different explanation for 
these behaviors. The results of Birch and Dow,” 
however, support the conclusion arrived at here 
because at higher temperatures, where rearrange- 


‘ments for a more stable configuration will take 


place more readily, they find the ‘‘abnormal”’ 
observed increase of compressibility with pres- 
sure becomes less pronounced and eventually 
disappears. 


UNIDIRECTIONAL PRESSURE 


If, instead of a hydrostatic pressure, a uni- 
directional pressure load 7 is applied to the 
system no external constraints will exist in the 
plane at right angles, resulting in an unbalanced 
energy distribution in the system. The spacing r 
along the line of thrust decreases by an amount 
r—ro and the strain energy increases, climbing 
up the left limb of the @ function plotted in 
Fig. 2. 

It is possible to set up hypothetical unsym- 
metric arrangements of attractive and repulsive 
forces such that on contraction of the longi- 
tudinal elements the effect laterally is an in- 
creased net attraction with consequent contrac- 
tion of these elements also. Such a condition 
yields a negative Poisson’s ratio, o, the ratio of 
lateral elastic extension to longitudinal elastic 
contraction under compression. A few examples 
of this phenomenon exist, but in general o is 

*P. W. Bridgman, Am. J. Sci. 237, 7-18 (1939). 


*F. Birch and R. B. Dow, Bull. Geol. Soc. Am. 47, 
1235-1256 (1936). 


positive and varies in magnitude between 0.25 
and 0.50. As a rule then, the symmetry will be 
such that, as the longitudinal elements are 
contracted, the effective lateral r for Eq. (1) 
decreases and thus the system extends in the 
plane at right angles to the compression. The 
effective left limb of the strain energy curve ¢ of 
Fig. 2 moves to the right a distance o(r,—70). 
When this extension reaches the value (7m—7o), 
the strain energy increasing to ¢m, the system 
becofnes unstable and, on further extension, 
ruptures by moving into the region of no stress. 
This “‘brittle’’ potential rupture, observationally 
called breaking along tension cracks, will occur 
across surfaces that tend to parallel the axis of 
compressive thrust. 

The ideal ‘‘potential” rupture condition, ac- 
cording to this presentation, is given by 


t2—to=(ro—%m)/o (5) 
and the elastic or potential strength by 
(F,)=Ar,-"—Br,. (6) 


Thus the elastic strength in compression 
should be 1/¢, or roughly three, times the elastic 
strength in tension. In general there is evidence 
of more or less ‘‘plasticity’’ in compressional and 
tensional test pieces of steel. Glass hard tool 
steel, in which a minimum of plasticity is ob- 
servable, has a compressive strength of about 30 
kilobars but even here there is some evidence of 
shear in the rupture. Direct comparisons for 
potential rupture cannot therefore be made but 
the tensile strength of this steel is not very 
different from the expected value. Such tests 
should be carried out at low temperatures. 

In considering the effect of hydrostatic pres- 
sure alone we had only one region available to 
the system, i.e., the entire region was subjected 
to the hydrostatic pressure ». Here, however, 
there are two regions available to our system— 
one subjected to compressive stress 7 and one to 
zero compressive stress. The system therefore, 
if free of internal and external constraint, will 
move into and occupy the region of lower energy. 
If then, before the lateral strain energy can 
attain the value ¢,, the longitudinal energy has 
increased to a value such that its thermodynamic 
potential becomes equal to that of the liquid for 




















the region corresponding to zero compressive 
stress, i.e., 


¢(T,2,p) =¢'(T,p), (7) 


the prime referring to the liquid phase, then the 
solid would melt to the liquid phase at (T,p) 
provided the stressed.region is permeable to, i.e., 
not be active on, the liquid. 

The thermodynamic potential is defined by the 


expression 
o 


c=y+er (8) 


and 0{/dr=r. 

A schematic plot of ¢ as a function of r at the 
stressed face is given in Fig. 3. ¢ will be observed 
to increase with increasing compression but to 
decrease with increasing tension until it reaches 
a minimum, fo, at fm. At the free face ({)rr 
=(Wo)rr and thus at this face ¢ increases with 
Yo for both tension and compression. It should 
perhaps be emphasized that in the direction at 
right angles to the stress y is always at its 
minimum value on the curve which however 
changes by virtue of the change in our expression 
for F. 

The use of the terms “liquid” and ‘‘melting”’ 
may seem confusing and should perhaps be am- 
plified. The process of the two-phase transfer 
type of mechanism suggested here involves the 
transition from an immobile solid phase to a 
mobile phase and thence back again to the solid. 
The mobile phase should more aptly be called 
the “fluid” phase and includes not only the 
condition ordinarily accepted as melting but also 
any other similar mechanism such as the one 
known as “migration of lattice points’ and 
diffusion. To make this point somewhat clearer 
the creep relations will be derived later for a 
solid in contact with a liquid in which it is 
slightly soluble and, for this case, the fluid phase 
is the solution. 

We have then two possible mechanisms to 
transport our system into the region of lower 
internal energy—one by snapping the cohesive 
bonds, determined by the strain energy poten- 
tial; and the other by transition to, and flow of, 
a fluid phase, solid—fluid—solid; this is deter- 
mined by the thermodynamic potential relations. 
The ‘‘fundamental strength’’ of our material is 
therefore not only a function of the strain energy 
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Fic. 3. Free energy for the stressed face vs. extension. 


potential but also of the thermodynamic po- 
tential. 


The internal energy of the system is merely’ 


a statistical mean of the lattice point energies 
and, on any surface, elements will be found with 
energy levels exceeding the mean value. The 
lattice structure is also more or less distorted 
because of impurities and other irregularities 
with resulting non-uniform stress distribution. 
Phase change is initiated at such regions of 
localized high energy values and results in still 
further localization of high stress areas. The 
process therefore tends to accelerate until the 
concentration of strain energy becomes too high 
to be borne by the remaining bonds and the 
specimen fails by shear, a combined “‘fluid’”’ and 
“brittle’’ release. 

An increase in compressive load thus acts 
similarly to an increase in temperature with 
respect to our ‘‘two-phase’’ mechanism. Phase 
change is initiated at, and proceeds from, loci of 
high energy levels at a finite rate, and not as an 
instantaneous disintegration of the lattice con- 
figuration. The time gradient of energy inter- 
change set up by change in the external condi- 
tions is also an important factor. The transfer of 
thermal energy will lag behind that of strain 
energy. This is evident from explosive phenomena 
wherein the rotational and vibrational energies 
can be observed to lag behind the translational 
energy. 
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The fact that this phase-change type of 
deformation is a function of time, and the elastic 
type is not, has important consequences. Failure 
as a result of a rapidly applied large compressive 
load is by brittle fracture across surfaces parallel 
to or making small angles with the axis of load; 
as the rate of application is decreased failure will 
tend to occur as shear along surfaces which are 
at larger angles to the axis of load; if the appli- 
cation is slow enough and the thermodynamic 
potential relations are favorable the specimen 
will deform by ‘‘gliding along 45° planes,’’ or by 
“flow.” 

If the load is not too great the shear may 
proceed in steps—‘‘melting” at high energy 
points, shear of the remaining cohesive bonds, 
slipping, recementation by solidification, the 
cycle repeating itself as conditions again become 
favorable. The most favorable surface for this 
process is the 45°-plane along which the resultant 
shear stress is a maximum and equal to one-half 
the compressive load. 

These statements need some qualification 
because the force and potential functions are not 
symmetric in space and therefore vary with 
orientation of the system. The curve of Fig. 1 is 
thus a function of a vector 7 and the lattice will 
therefore, if conditions are favorable, rupture or 
shear along surfaces across which the cohesive 
bond, Fn, is the least such as cleavage, parting, 
and twinning planes. 

““Melting’”’ can occur in specimens under ten- 
sion, viz. along shear planes, but the phenomenon 
should not be so evident here because it will take 
place only at the free surfaces (see Eqs. (11b), 
(12b), (14b)). 


DERIVATION OF EXPRESSIONS FOR ‘‘CREEP”’ OR 
“PLASTIC FLow” 


The creep relations will be derived first for a 
solid under-compressive load and immersed in a 
liquid in which it is somewhat soluble. 

Assume an initially ‘ideal solid” in which the 
thermodynamic potential is the same for all the 
faces and equal to the potential of the solid in 
solution (solute), i.e., assume that the solid is in 
equilibrium with the saturated solution at tem- 
perature T and hydrostatic pressure p. If now 
the solid is loaded by a longitudinal compressive 





force + the thermodynamic potentials at the 
stressed and free faces of the solid and of the 
solute will no longer be equal and the system no 
longer in equilibrium so long as the stress exists. 

In order to derive our thermodynamic equa- 
tions the system is first divided into hypothetical 
isolated parts, namely the regions at the stressed 
surface, at the free face, and of the solution bulk. 
We also assume any coexistence of phases neces- 
sary for our derivations. The physical interpreta- 
tion follows readily from this procedure and 
moreover we avoid any confusion that might be 
introduced from a compromise! between the 
thermodynamic and the physical picture. 

The following formulae were derived by ordi- 
nary thermodynamic methods applied to stressed 
systems” and have been somewhat simplified 
for purposes of clarity. 


We have 
(a In a2/ dr) sr = M/pRT (9a) 
and 
(0 In a2/0T) FF = Mr/pERT, (9b) 


where Eq. (9a) refers to the stressed face (sub- 
script SF) and (9b) to the free face (subscript 
FF). a2 denotes the activity of the solute, 7 the 
compressive stress, M the mole weight, p 
the density of the solid, R the gas constant 
(R=83.156 bar cm), and T the absolute tem- 
perature. E is Young’s modulus of elasticity in 
compression, i.e., E=dX/de where X is the 
stress (negative for pressure) and e the extension 
per unit length in the direction of the stress. For 
small stresses E may be approximated by a 
constant, generally of the order of 10° bars. 

On integrating (11a) we have, at the stressed 
face, 


mt —mo=(pRT/M) In (a2/a2,), (10a) 


where p is the mean value for the integration 
limits. 


1 P, W. Bridgman, Phys. Rev. 7, 215 (1916), derives an 
expression which may be correlated with (14a) and (14b). 
E. D. Williamson, Phys. Rev. 10, 275 (1917). H. C. 
Boydell, Ec. Geol. 21, 1 (1926). Boydell derives Poynting’s 
expression which applies to a solid under hydrostatic 
pressure ~: and liquid at p2 where p:>2. Poynting’s 
expression has the same form as (14a) of this paper. 

122R. W. Goranson, Thermodynamic Relations in Multi- 
component Systems (Carnegie Inst. Washington. Publ. No. 
408, 1930). For usage of activity see R. W. Goranson, J. 
Chem. Phys. 5, 107-112 (1937). 
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If the system is impervious to liquid and the 
load contact is a ‘‘perfect surface’ then only the 
external free face can be in contact with the 
solution. Actually the system is, in general, 
never impervious to the liquid and perfect con- 
tact surfaces non-existent. Penetration occurs 
along crystal boundaries, across crystal cleavages 
and cracks. The stressed surfaces carry the load 
non-uniformly and contain high points of 
localized stress. As these high points are dis- 
solved the stress distribution shifts to other high 
points and the process is repeated. 

Thus at the stressed face solubility is increased 
by a compressive load and lowered by a tensile 
load. At the free face solubility is increased by 
either a compression or tension but by a much 
smaller factor. It should be noted however that 
the integral expression for (9b) is 


w?—m)? pRT de 
E M 2% 


which may become large for large 7. p and E are 
here mean values for the range of integration. 
This relation is the plastic flow factor for tension. 

According to this interpretation it is found 
that, whereas under compressive load the solute 
concentration tends to increase at both stressed 
and free face, diffusion of solute away from these 
faces would leave the bulk of the solution super- 
saturated. Consequently deposition of solute 
should occur at unstressed places. If, however, 
the solution can supersaturate by an amount in 
excess of the relatively small increased solubility 
at the free faces, then crystallization will take 
place on these free surfaces which, under these 
conditions, act as nuclei for deposition. Our 
mechanism can thus be considered as a diffusion 
process along crystal grain boundaries wherein 
the solution acts as a transfer medium. This dif- 
fusion rate will approach a steady state when a 
dynamic equilibrium between rate of solution 
and rate of deposition has been established, and 
will be governed by the mobility of the solute in 
the solution, the path length, and the concen- 
tration head between the stressed and free faces. 

If, by changing the load, we change only the 
concentration head of the polycrystalline sub- 
stance then the change in the steady creep rate, 
é,, should be directly proportional to the change 


(10b) 
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in activity, or 


dln €; 
=B , (11) 
dr dr 


d\n ag 





where B is a constant. From (10a) and (11) we 
obtain 


a—mo=K |n (és/€0), (12) 


where K is a physical constant. This is the rela- 
tion obtained empirically'* for the steady creep 
rate of compressively loaded alabaster in contact 
with water. 

As compressive load is increased the solubility 
at the stressed surface eventually becomes 
larger than the amount by which the solution 
can supersaturate. Under these conditions the 
solution will drop the excess solute in any avail- 
able unstressed space. Consequently, the texture 
of the substance may thus become so loosened 
that it crumbles under the load. 

The remainder of the discussion parallels that 
for plastic flow in polycrystalline metals. As 
before, we divide our system into hypothetical 
isolated portions and derive expressions for the 
stressed face and free face. We assume also the 
coexistence of two phases—the bound atoms of 
the solid and the free migrating atoms with 
energies exceeding those necessary to break the 
cohesive bonds of the crystal lattice. At the 
melting point this will be equivalent to the heat 
of melting. 

The analogous relations to (9a) and (9b) con- 
necting compressive stress 7 and ‘melting’ 
temperature, 7;,, in degrees absolute are readily 


derived as 
(“ In =) 
dr SF 


din T» 7 
(——) Sa, (12b) 
dr FF pEh 


where fh is the heat of ‘‘melting’”’ and the other 
quantities are as before. 

The melting point at the stressed face is thus 
depressed for compression and raised for tension. 
This may be observed graphically from Fig. 3 


13 The experimental data were obtained by Griggs. To be 
offered for publication, with Goranson, in a future issue 0 
the Bulletin of the Geological Soctety of America. 
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where the thermodynamic potential ¢ is seen to 
increase with pressure and decrease with tension 
at the stressed face. The melting point at the 
free face is depressed for both compression and 
tension. Here (yo)rr=({)rr always increases. 
The ratio of lowering at the two faces is the 
same as that found for (9a) and (9b). 

The number of atoms with Maxwellian ener- 
gies in excess of that needed to break the lattice 
bonds, i.e., the number of free migrating atoms, 
increases with increase in temperature. At the 
melting point this fraction of the total number 
becomes unity. The factor 


6=T/Tn, (13) 


where JT denotes the constant experimental and 
T, the melting temperature, has been used in 
correlating temperature versus creep data for 
different metals. The behavior of a low melting 
point metal is thus considered as equivalent in 
behavior to a high melting point metal at a cor- 
respondingly higher temperature to give the 
same 6, other factors remaining unchanged. The 
argument should be carried out in terms of 
energy but, as compensation is here effected 
through the / term this expression is a sufficiently 
satisfactory approximation. 

On combining (12a) and (13) and integrating 
we have 


1 —1o= ph |n (0/6), (14) 


where p and fare here the mean values over the 
range of integration. 

For a sufficiently large load, i.e., for relatively 
large creep rates, we may in the same manner as 
before write 


4—mo=K' In (é,/é0), (15) 


where K’ is a physical constant and é, the steady 
state or minimum creep rate. This is the ex- 
pression that has been used for about thirty 
years to express the empirical relation between 
yield stress and minimum creep rate é, in metal 
mosaics. 

The expressions (12) and (15) were derived for 
relatively large creep rates and thus for com- 
pressive loads large enough to iron out the initial 
inequalities over a short period of time. Let us 
assume now that we are operating with either 
very small stresses or with a substance that does 


not “‘migrate’’ readily. For these conditions the 
number of migrating units will be relatively few 
and their paths much shorter. Thus, instead of 
being able to picture a statistical streaming 
action, we are slowed down to a hop-skip process 
and can no longer set up a steady mean flow. 
The crystallization process must be explicitly 
considered here and for these cases we have for 
mobility along the stressed surfaces. 


K (0/00) =Ke(—*0/oh (16a) 


and for crystallization along the same surfaces 


K (00/0) = Ke-(*— 70 /ph, (16b) 


The creep rate will be given by the difference 
between these two quantities or 


T To 
(€/€)) =K sinh ——. (17) 
ph 


The same reasoning and therefore a similar sinh 
expression also replaces (12) for very small creep 
rates. 

An expression of this type was obtained em- 
pirically'4 by combining analytically the empirical 
logarithmic relation for high creep rates (see Eq. 
(15)) with the linear relation between stress and 
creep rate observed for very low creep rates. 

So far we have considered’ only the relation 
between steady creep rate and stress. When 
compressive load is first applied the irregularities 
in the structure of the polycrystalline material 
set up an initial localization of stress at the 
raised points and thus, for large enough loads, 
in an initially high creep rate. This rate gradu- 
ally diminishes toward a steady state as the 
original inhomogeneities of texture become 
ironed out and the stress redistributes itself over 
larger and larger surface areas. The effective 
stress, for constant load, thus decreases with 
time. Finally the pore spaces in the texture 
become filled and the grains more or less reori- 
ented into their most stable crystallographic 
configuration for an axial load distribution. If 
now a series of such creep-time tests is made for 
varying compressive loads on initially identical 
specimens this last state should be reached at 


14 For example see H. Mussmann, Ann. d. Physik 31, 130 
(1938). 
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Fic. 4. Tensile creep tests on coarse-grained lead rods at 
room temperature. By suitable changes in coordinates this 
diagram can be made to express, to a fair approximation, 
the compressive creep tests of alabaster in water. 


about the same longitudinal contraction é€m 
provided none of the material has precipitated 
elsewhere. At this point the specimen has 
stabilized itself and from this point on should 
therefore behave more or less as a unit, i.e., 
approximate to the condition for a single crystal. 
Solution (or migration) then begins to work 
along crystallographic surfaces and will be most 
effective on surfaces across which the cohesive 
bonds are smallest, i.e., have the largest ¢, in 
the general 45° trend to the stress axis. If the 
load is large enough the specimen will then 
deform by a gliding action along 45° planes and 
the creep rate, measured by the rate of con- 
traction, will appear to accelerate until the 
specimen “‘fails.’’ This is graphically illustrated 
in Fig. 4.!® 

Plastic flow may also occur in single crystals 
but the effect in general is a gliding along 
cleavage, twinning, or parting planes from com- 
bined ‘‘melting’’ and snapping of bonds. Local- 
ized high stress regions may be set up in crystals, 
as was mentioned earlier, and thus an initially 
single crystal may, under a compressive load, 
finally become a mosaic of reoriented crystalline 
grains. 

The energy at the external crystal faces cannot 
be denoted by the energy of the interior because 
of loss of symmetry at these external faces; the 
difference is known as surface tension. For ionic 
crystals this difference will be greatest along 
edges and at corners and least at the centers of 


( on taken from J. McKeown, J. Inst. Metals 60, 207 
1 k 


GORANSON 


faces; this is the reason for the initial skeleton 
growth of such crystals. With a high symmetry 
type of structure this difference might become 
greatest at re-entrants. For our system solution 
or melting will be most rapid at the stressed 
surfaces of highest ¢ and growth by crystalliza- 
tion most rapid on the free faces of lowest ¢. 
These directions are, in general, indicated by the 
crystalline form and cleavage. For example, mica 
and related minerals will tend to grow with two 
crystal axes (cleavage planes) perpendicular to 
the axis of compressive load, asbestos and 
related minerals to elongate (one crystal axis) 
in the plane perpendicular to the stress. The best 
examples we have are the metamorphic rocks 
which have recrystallized according to this 
mechanism under high confining pressure (see 


below). 


ELAstTic AFTERWORKING 


There is a type of deformation which on load- 
ing is a contraction and on unloading a recovery 
by extension, according to the expression 


é=«,B/(t+7), 


where é denotes the strain rate, « the purely 
elastic portion of the strain, B a constant, ¢ the 
time and 7 a time constant. This type of re- 
covery, after unloading, is called “elastic after- 
working”’ and has been observed in glass fibers, 
in steels, and in rocks. 

This is the kind of deformation to be expected 
from materials which are relatively strong elas- 
tically. For these materials the strain is mainly 
an elastic one but, for a long continued applica- 
tion of a moderate load, a small amount of 
migration of lattice points will take place if the 
load is left on long enough. On release of load 
the material tends to recover elastically, i.e., 
instantaneously, but it cannot recover comi- 
pletely because under load lattice elements have 
migrated and then solidified in conformity with 
the equilibrium conditions existing while under 
load. These therefore formed bonds tying in this 
equilibrium state so that when load is released 
an opposing stress distribution is consequently 
initiated and the initial equilibrium conditions 
are attained only by a backward migration of 
these same elements in retracing their paths. 
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Work HARDENING 


The mechanism of work hardening should not 
require particular comment. There is a further 
consequence of this phenomenon: The reorien- 
tation and reorganization of a polycrystalline 
aggregate for greater elastic stability to a com- 
pressive load, for example, means that the ma- 
terial should then be weaker for a tensional load 
along the same axis. This is observable in push- 
pull load-extension diagrams of ‘‘overstrained”’ 
steels. 


EFFECT OF HYDROSTATIC PRESSURE 
ON STRENGTH 


If hydrostatic pressure is applied to a crystal 
lattice, r decreases by an amount ro—rp=€p 
and if the electrostatic force field is symmetric 
in space the decrease will be uniform in all direc- 
tions. If now an additional stress is applied by 
means of a superimposed unidirectional com- 
pressive load 7 then r along this axis will decrease 
to rx by the additional amount e,=7r,—7p. 

The lattice, assuming ‘normal’ behavior, 
expands laterally doing work against the con- 
fining pressure until a new equilibrium is estab- 
lished between the repulsive, attractive, and 
external forces. 

It might be anticipated that Poisson’s ratio 
should bear a functional relation to r. The 
evidence from seismic data indicates, however, 
that within the earth o remains approximately 
constant, at about 0.27, independent of pressure 
(depth). 

The lateral extension will increase as the 
longitudinal compressive load is increased and 
thus r will gradually move back through 7o to 
tm at which point the specimen becomes unstable 
and ruptures. The load for this condition will 
tepresent the compressive strength of the 
specimen. 

Rupture, here the “brittle” potential type of 
fracture, will occur then for the condition 


r=—9 dy 
f —d F—(t'm—T1o) 


F=—- dy F=<0 dF 
f —dF= 
F=—p dF o 


Where [(a)m— p] denotes the “elastic” com- 
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pressive strength of the specimen under confining 
pressure p. 
The slope of the r—F curve, 


dr vith 
dF jC;—iCgi-' 


flattens rapidly as r decreases and approaches 
zero as a limit. This initial change in slope may 
be observed graphically in the curves of com- 
pressibility plotted as a function of pressure. 
Hence the smaller the 7, and consequently the 
higher the hydrostatic confining pressure, the 
larger must be the unidirectional thrust for the 
same %?m—%, but this distance, as a matter of 
fact, is itself also increasing by virtue of increas- 
ing p. The compressive strength of the specimen 
should therefore accelerate rapidly with increas- 
ing confining pressure and become infinitely 
strong elastically as the confining pressure con- 
tinues to increase indefinitely. If, therefore, the 
compressive strength be plotted as ordinate and 
hydrostatic confining pressure as abscissa the 
curve should be found to rise gently at first and 
then to steepen rapidly, eventually becoming 
infinitely steep. 

The above conclusions from this theory, de- 
pending entirely on seismic data, have been held 
in abeyance until they could be further confirmed. 
Recently confirmation has been obtained from 
the experimental work of Griggs.1® Experimental 
work is also being conducted at the Geophysical 
Laboratory to verify the calculated strengths for 
steel. The initial pressure effect is complicated by 
the geometry of cracks but these cracks close up 
at a few thousand atmospheres pressure. 

The curve of compressive strength versus 
hydrostatic confining pressure for substances 
which have polymorphic pressure modifications 
should exhibit discontinuities at such trans- 
formation points. In fact if such transformations 
took place rapidly while the substance was under 
unidirectional load the structure should break 
down. These points represent therefore loci of 
instability, so that the above conclusions would 
be applicable only to the homogeneous regions 
above and below such loci. 


16D. T. Griggs, J. Geol. 44, 541-577 (1936). D. T. 
(os, and J. F. Bell, Bull. Geol. Soc. Am. 49, 1723-46 
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Thc energy levels rise very steeply and even- 
tually for some value >>? a state will be reached 
for which 


¢(T, 2, p= S(T, p), 


where these quantities are as before. 

If this condition becomes fulfilled failure will 
then take place by plastic deformation or ‘‘plastic 
shear.’’ In other words, a substance under axial 
compression will become not only stronger elas- 
tically but also more malleable with increasing 
hydrostatic confining pressure. This is shown 
very clearly from the experimental results of 
F. D. Adams and Coker.” Their results were not 
duplicated by Griggs but the reason for it should 
be fairly obvious. Adams and Coker’s confining 
pressure was obtained by means of a pliable steel 
wall, Griggs’ by means of a thin liquid. As soon 
as incipient melting takes place along, say, a 


17F, D, Adams and E. G. Coker, An Investigation into the 
Elastic Constants of Rocks (Carnegie Inst. Washington, Publ. 
No. 46, 1906). F. D. Adams and J. A. Bancroft, J. Geol. 
25, 597-658 (1917). 
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shear plane the substance is relatively free to 
shear in Griggs’ case because the only work 
involved is (a) that of breaking the remaining 
cohesive bonds across this plane and (0) that of 
exchanging regions with the pressure fluid, since 
the resistance to flow of the pressure fluid, for a 
finite rate, is a simple function of its viscosity. 
For Adams and Coker’s case we have, in addition 
to (a), the work of pushing aside the supporting 
steel wall which is here the total lateral support. 
Release of strain energy can take place rapidly 
in the former case but only at a slow rate in the 
latter case. 

Thus, although Griggs has expressed the belief 
that his results are more applicable to problems 
of geological deformation in depth and that 
rocks do not flow in the dry state, the conclusions 
arrived at herein do not support this contention 
as a general conclusion. These two apparently 
dissimilar experimental results are special cases 
of the same physical hypothesis described in 
this paper. 
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An empirical equation is given for the isothermal variation with density of the work content 
of pure hydrocarbons in the gaseous or liquid state. From this fundamental equation are de- 
rived (a) an equation of state, (b) an equation for the fugacity, and (c) an equation for the 
isothermal variation of the enthalpy. These equations summarize P-V-T properties of the 
gaseous or liquid phase, critical properties, vapor pressures, and latent heats of evaporation. 
A procedure is suggested for determining numerical values of the parameters in the equation. 
Such values are given for methane, ethane, propane, and u-butane. A comparison is made 
between observed properties of these hydrocarbons and those predicted by the equations. 


A. INTRODUCTION 


ECENT experimental studies of pure light 
hydrocarbons by Sage and Lacey and co- 
workers and by Beattie and co-workers are useful 
in developing an equation to represent the ther- 
modynamic properties of these substances. An 
equation for this purpose has several advantages: 





it permits interpolation of experimental data; it 
facilitates thermodynamic calculations involving 


‘integration and differentiation; it provides a 


concise summary of a large mass of data; and it 
provides a point of departure for the treatment 
of the thermodynamic properties of mixtures. 

Concurrent with this advance in experimental 
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knowledge there has taken place development of 
statistico-mechanical methods for deducing the 
thermodynamic properties of fluids in bulk from 
the interaction energies of aggregates of the 
constituent molecules. A method for the rigorous 
treatment of this problem, as applied to gases, 
was first given by Ursell! and has recently been 
extended and clarified by Mayer and co- 
workers.” A less rigorous treatment, leading to 
a single equation describing the thermodynamic 
properties of both the liquid and gaseous phases, 
has been given by Lennard-Jones and co- 
workers.* These theoretical treatments are par- 
ticularly valuable for suggesting the correct 
qualitative features of suitable equations.. They 
do not lead to simple expressions in terms of 
familiar functions, and in the present state of 
our knowledge of molecular fields, they do not 
lead to correct numerical values for the proper- 
ties of matter in bulk. It has, therefore, seemed 
worth while to attack the problem of representing 
the thermodynamic properties of hydrocarbons 
from a more empirical standpoint, and to use the 
theoretical results only to suggest the proper 
general characteristics of the formulation. 
These theoretical analyses demonstrate that 
the phase integral summarizing all thermo- 
dynamic properties of a fluid may be expressed 
as the product of two factors, one arising from 
the internal degrees of freedom of each molecule 
and the other arising from the interaction ener- 
gies between aggregates of two or more mole- 
cules. From the first factor may be derived 
equations for the specific heats, entropy, work 
content, and free energy of the fluid in the 
hypothetical ideal gas state. From the second 
factor may be derived equations for the difference 
between the thermodynamic properties of a real 
gas or liquid and the corresponding properties of 
the fluid in the hypothetical ideal gas state. The 
object of the present paper is to present an em- 
pirical equation for this second class of properties. 


Definition of terms 


Our equation expresses the isothermal varia- 
tion with density of the work content A of a mole 


. Ursell, Proc. Camb. Phil. Soc. 23, 685 (1927). 

* Mayer, J. Chem. Phys. 5, 67 (1937); Mayer and Acker- 
(1938) Ibid. 74; Mayer and Harrison, Ibid. 6, 87, 101 

* Lennard-Jones and Devonshire, Proc. Roy. Soc. Lon- 
don A163, 53 (1937). 





of gaseous or liquid hydrocarbons. The work 
content is also known as the Helmholtz free 
energy or the Gibbs y-function; it is related to 
the internal energy E, absolute temperature 7, 
and entropy S by the equation 


A=E-TS. 


We are concerned with the residual work content 
A, here defined as the difference between the 
work content of a mole of hydrocarbons at molal 
density d and absolute temperature T and its 
work content in the hypothetical ideal gas state 
at the same density and temperature. The 
residual work content A may be defined more 
rigorously by the equation: 


A=A-—RT In d—lim (A—RT Ind). (A1) 


The method of Ursell and Mayer may be used to 
prove that the limit occurring in Eq. (A1) exists 
and is finite. 

An equation of the form (A1) is a fundamental 
equation in the sense that from it a number of 
thermodynamic properties of the fluid may be 
computed, no one of which alone would permit 
the computation of all the others. The thermo- 
dynamic properties of a fluid which may be 
computed from the fundamental equation of this 
paper include: (1) P-V-T properties of the gas 
or liquid phase. (2) Critical properties. (3) The 
effect of pressure on the enthalpy. (4) The effect 
of pressure on the free energy; i.e., the fugacity. 
(5) Vapor pressures. For instance, an equation 
of state in which the pressure P is expressed as a 
function of the density and temperature may be 
obtained by the relation: 


P=RTd+4?(dA/dd)r. (A2) 
The fugacity f is given by: 
RT In f=A+RT In dRT+P/d—RT, (A3) 


and the isothermal variation of the enthalpy H 
by: 


. a 
H(T, d)—H(T,0)=— r( 





A/T 
) +P/d—RT. 
OT /a (A4) 
Conversely, the residual work content of the 
gas phase may be obtained by integration of an 
equation of state: 


@ P—RTd 
A= f altnonel* } (AS) 
0 d? 
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Relationship to other equations of state 


Of the many equations of state that have been 
proposed perhaps the best known are the van 
der Waals and the Beattie-Bridgeman equations. 
The van der Waals‘ equation, 


RTd 
1—bd 


P= (A6) 





—ad’, 


has the advantage of simplicity, and for this 
reason is often used in discussing the thermo- 
dynamic relationships of the gaseous and liquid 
states and the critical region between them. By 
substituting the van der Waals equation of state 
(A6) in (A5), Eq. (A7) for the residual work 
content of the gas phase is obtained: 


= —RT In (1—bd) —ad. (A7) 


When derived in this way, Eq. (A7) applies 
rigorously only to the gas phase. In the van der 
Waals theory of the continuity of the liquid and 
gaseous states the assumption is made that this 
equation is valid for the liquid phase, as well. 
From this standpoint, Eq. (A7) is the funda- 
mental equation of the van der Waals theory for 
both the gaseous and liquid phases. It may be 
used to calculate the equation of state and fuga- 
city of either phase, and from them, the vapor 
pressure. Vapor pressures calculated in this way 
from the van der Waals equation are in only 
qualitative agreement with observed vapor 
pressures. This need not invalidate the sig- 
nificant assumption of the van der Waals theory, 
that a single fundamental equation can be found 
for both the liquid and gaseous phases, as the 
discrepancy may be due to the recognized failure 
of the van der Waals equation to correspond 
accurately to observed P-V-T data. 

Among the many equations of state which 
have been proposed as improvements on the van 
der Waals equation, the Beattie-Bridgeman® 
equation, , 


P=RTd+(BoRT—Ao—Re/T?)d? 
+(—BobRT+Aca—RBoc/T?)d? 
+RBobcd*/T?, (A8) 


*Van der Waals, Die Continuitét des Gasférmigen und 
Flissigen Zustandes (Barth, Leipzig, 1899). 

5 Beattie and Bridgeman, Proc. Am. Acad. Arts Sci. 63, 
229 (1928). 
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has, perhaps, most satisfactorily combined 
mathematical simplicity with accurate represen- 
tation of P-V-T data over a wide temperature 
range. In its present form the Beattie-Bridgeman 
equation does not correspond accurately with the 
P-V-T properties of liquids or gases above their 
critical density, so that it cannot be used to cal- 
culate the residual work content above this 
density or the vapor pressure. 


The present equation 


The equation of state proposed in this article 
may be regarded as a modification of the Beattie- 
Bridgeman equation to enable it to represent 
more accurately the properties of fluids at high 
densities. The van der Waals assumption of the 
continuity of the liquid and gaseous states is 
retained, in that a common fundamental equa- 
tion and a common equation of state are used 
for both gaseous and liquid phases. The funda- 
mental equation proposed for expressing the 
dependence of the residual work content A of a 
mole of gaseous or liquid hydrocarbon on the 
molal density d and absolute temperature T is: 


A =(BoRT—Ao—Co/T?)d 
+ (ORT —a)d?/2+aad*/5 
cd*71—exp (—~yd?) exp (—yd?) 
=| dt —~é«iS 





} 9 


The equation of state corresponding to (A9) is: 
P=RTd+(BoRT—Ao—C)/T?)d? 
+(dRT —a)d'+aad® 
FS cd*(1+-yd?) exp (—~d?*) 
T? 





(A10) 


The isometrics of this equation of state have 
the same form as the isometrics of the Beattie- 
Bridgeman equation: 


P=Ty(d) — $(d@) -T(@)/T”. 


The parameters Bo, Ao, and Co of Eq. (A10) play 
the same role and have about the same numerical 
values as corresponding parameters in the 
Beattie-Bridgeman equation. The parameters 0, 
a, c, y, and @ are not directly related to any oc- 
curring in the Beattie-Bridgeman equation. 


(A11) 
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B. CONSIDERATIONS INFLUENCING ADOPTION OF 
PRESENT EQUATION 


The work of Ursell and Mayer may be used to 
prove that the residual work content A of a gas 
can be expressed by the series 


A=Q:(T)d+Q2(T)P+-:-, (B1) 


convergent in the region of stable existence of the 
gas phase. An equation for both the gaseous and 
liquid phases was sought which would reduce to 
(B1) at low densities, which would constitute a 
suitable compromise between simplicity of ap- 
plication and accurate correspondence with ob- 
served thermodynamic properties of gases and 
liquids, which would contain a minimum number 
of independent parameters to be fitted numeri- 
cally, and which would summarize as well as 
possible the following classes of properties: (1) 
P-V-T properties of the gas phase. (2) Critical 
properties. (3) Vapor pressures. It is believed 
that emphasis on these properties results in the 
most generally satisfactory equations for pure 
light hydrocarbons on which to base the predic- 
tion of liquid-vapor equilibria of their mixtures. 
This is the ultimate object of the project de- 
scribed in part in this article. 


TEMPERATURE, °C 


8 


150 


a 


200 225 


The functional form of the equation finally 
adopted as most satisfactorily filling these three 
requirements was determined primarily from an 
examination of the P-V-T properties of the gas 
phase; then, for each hydrocarbon, the best 
numerical values of the parameters in the equa- 
tion were estimated from a consideration of all 
three classes of properties. 


Equation of state 


In developing the equation of state, attention 
was first directed to the isometrics of (P — Rtd) /d? 
versus temperature; these are very nearly 
straight. For each hydrocarbon it was found that 
these isometrics could be represented by an 
equation of the form: 


(P—RTd)/d?=RTB(d)—A(d)—C(d)/T?. (B2) 


These isometrics are illustrated with the example 
of propane in Fig. 1. The plotted points are 
derived from the observed pressures of Beattie, 
Kay and Kaminsky;® the lines represent the 
equation of state for propane given in Section D. 

Equations, consistent with expression (B1), 
were then sought to express the dependence of 


6 Beattie, Kay and Kaminsky, J. Am. Chem. Soc. 59, 
1589 (1937). 
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Fic. 1. Isometrics of gaseous propane. 
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Fic. 2. Functions of density in equation of state. 


B(d), A(d), and C(d) on density. It was found 
that B(d) could be represented by the simple 
linear equation 


| 


B(d) = By+6d, (B3) 


but that plots of A(d) and C(d) against density 
showed marked curvature. A(d) was finally 
represented by an equation of the form: 


A(d) =A,»+ad(1—ad’*), (B4) 
and C(d) by an equation of the form: 
C(d) = Co—cd(1+ yd?) exp (— yd*). (B5) 


A simple polynomial would have sufficed to 
express C(d) for the gas phase, but in order to 
represent the fugacities of the liquid phase at 
high densities and low temperatures, the more 
complicated expression (B5) was required. 

The functions B(d), A(d), and C(d) for propane 
are illustrated in Fig. 2. The same qualitative 
features of these functions were found for each 
hydrocarbon studied. The properties of C(d) are 
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of especial interest. Since 
(0?P/dT?*) a= —6d?C(d)/T*, 


—C(d) has the sign of the curvature of the 
isometrics. It appears to be a general rule that 
the isometrics of the gas phase at densities less 
than the critical have negative curvature.’ Near 
the critical density the curvature of the iso- 
metrics becomes positive, as first noted by 
Sydney Young.’ The isometrics of the compressed 
liquid phase at still higher density have negative 
curvature, as shown by Bridgman.*’ These 
properties of the isometrics of real gases and 
liquids are paralleled by Eq. (B5) for C(d). As 
shown in Fig. 2, — C(d) is negative up to approx- 
imately the critical density. Between this 
density and about twice the critical density, it is 
positive. Above twice the critical density, in the 
region corresponding to the compressed liquid, 
C(d) is again negative. 

When expressions (B2), (B3), (B4), and (B5) 
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Fic. 3. Isotherms of propane. 


7 Young, Proc. Phys. Soc. London 13, 602 (1895). 
8 Bridgman, Proc. Am. Acad. Arts Sci. 66, 185 (1931). 
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are combined, the equation of state (B7) is 
obtained. This is the isometric form of 


P=RTd+d*) RT[Bo+bd]—[Ao-+ad—aad*] 
1 
~ppee nl oF exp (—7d?) |}, (B7) 


Eq. (A10), useful in computing pressures at 
constant density. This equation is integrable and 
differentiable with respect to density, volume, 
and temperature. Its eight parameters is the 
smallest number with which the P-V-T properties 
of the gas phase could be represented up to 
about twice the critical density. In Section E it 
is shown that the average error in calculated 
pressures up to this density is only 0.34 percent. 
Owing primarily to limitations in Eq. (B4) for 
A(d), the error at higher densities becomes much 
larger; better agreement at higher densities 
could not have been obtained without the use of 
several more parameters. 

To illustrate in terms of familiar properties 
the agreement of equation of state (B7) with 
observed data, isotherms of density versus 
pressure for propane are plotted in Fig. 3. This 
figure is typical of all hydrocarbons studied. The 
complex dependence of density on pressure, 
especially near the critical temperature, 96.8°C, 
is clearly shown. Up to 9 moles per liter, about 
twice the critical density, the equation of state 
accurately reproduces this dependence. Above 
this density, the equation predicts densities too 
low by from 2 to 3 percent, because of the 
limitations of the equation noted in the preceding 
paragraph. 


Critical properties 


Since the equation of state (B7) reproduces 
observed P-V-T properties of the gas phase in 
the critical region, it necessarily satisfies our 
second requirement, that it summarize observed 
critical properties. Except for methane, calcu- 
lated and observed critical temperatures and 
pressures agree within 0.25°C and 0.4 atmos., 
respectively. However, it should be noted that 
our equation describes what may be termed the 
“classical” picture of the critical region between 
liquid and vapor; at each temperature below the 
critical point it leads to only one pair of phases of 
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Fic. 4. Vapor pressure of propane. 


different densities which can be in equilibrium. 
This picture differs from the recent observations 
of Maass? and the predictions of the theory of 
Mayer,? which indicate that for several degrees 
below the critical temperature it is possible 
for an infinite number of phases of different 
density to coexist. 


Vapor pressures 


The fundamental equation (A9), summarizing 
P-V-T properties of both phases, critical proper- 
ties, and vapor pressures, was developed from 
the equation of state as follows. 


(1) An equation for the residual work content 
of the gas phase was obtained by substitution of 
equation of state (B7) in the thermodynamic 
relation (A5). 

(2) It was then assumed that the resulting 
Eq. (A9) was valid for the liquid phase also, as 
in the van der Waals treatment of the continuity 
between the gaseous and liquid states. 

This assumption may be tested by seeing how 
well vapor pressures computed from an equation 
of the form (A9) agree with observed data. 
Detailed comparison between observed and cal- 
culated vapor pressures is deferred to Section E. 
Here it may be stated that the average error in 
calculated vapor pressures above one atmosphere 
is only 1.1 percent. This justifies the use of Eq. 
(A9) as a fundamental equation for both the 
gaseous and liquid phases. 

Figure 4 illustrates vapor pressures of propane 
predicted by the equation and those observed by 
various authorities. Approximately the same 


9 Maass, Chem. Rev. 23, 17 (1938). 





340 BENEDICT, 


agreement is found for vapor pressures of 
methane, ethane, and n-butane. 

To calculate vapor pressures from the funda- 
mental equation (A9), Eq. (B8) for the fugacity 
f is first obtained from it by the thermodynamic 
relation (A3). 


RT In f=RT In dRT+2d(BoRT —Ao— C/T”) 
+3d?(bRT —a)/2+6aad*/5 
cd*. 1—exp (— yd?) 
= yd? 








—vyd? 
Pc tit a PO exp (- re") | (B8) 


A gaseous density d® and a liquid density d” are 
found by trial such that 


P(d°, T) =P(d*, T) (B9) 


and 
RT I|n f(d%, T) =RT In f(d4, T). 


P is the calculated vapor pressure at absolute 
temperature 7. , 


(B10) 


Isothermal variation of enthalpy 

From fundamental equation (A9) and ther- 
modynamic relation (A4) it is possible to derive 
an equation for the change in enthalpy H with 
density at constant temperature. We obtain: 
H(T, d)—H(T, 0) =(BoRT —2A9—4C,/T?)d 

+(2bRT —3a)d?/2+6aad*/5 

“|s —exp (—7d?) 


ime 
a” yd? 





Vv 


ome d? 
_ exp (—7 +a exp (— 14) (B11) 


Eq. (B11) may be used to compute the change 
in enthalpy of the gas upon isothermal expansion 
or to compute the latent heat of evaporation of 
the liquid. In the latter application, the densities 
of the saturated liquid and vapor which satisfy 
Eqs. (B9) and (B10) above should be used, 
rather than the observed saturation densities. 


Simplified equations for the gas phase 


It is an advantage of the present formulation 
that simpler equations can be used when the 
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TABLE I. Comparison of simplified and complete equations 
for propane. Density =2.565 moles /l. 
Temp. = 96.81°C. 











EQUATION SIMPLIFIED COMPLETE 
Pressure, atmos. 39.108 39.121 
Fugacity, atmos. 26.737 26.743 
Enthalpy, /-atmos./mole — 55.134 —55.147 








density does not exceed one-half the critical 
density. When fundamental equation (A9) is 
expanded in ascending powers of the density, we 
obtain: 


A = (BoRT —Ao— Co/T?)d 
+(bRT—a+c/T?)@/2 
+aad5/5—cy*d*/12T?+-+>. 


At low densities the terms in the fifth and higher 
powers of d may be neglected in comparison with 
the others. We then have for the fundamental 
equation: 


(B12) 


A=Q,d+Q.d’, 
where 
—_ anes oni 2 
Q:1=B)RT—Ao—C,/T?, (B13) 
Q2=(bRT—a+c/T?)/2; 
for the equation of state: 
P=RTd+(Q,d?+2Q2d*; B(14) 
for the fugacity: 
RT |n f=RT In dRT+2Qid+3Q2d?; (B15) 


and for the isothermal variation of the enthalpy: 


H(T, d)—H(T, 0) =(BoRT—2A0—4Co/T?)d 
+ (2bRT—3a+5c/T?)d?/2. (B16) 


The fundamental equation (B13) for the gas 
phase at low densities is in the form (B1) required 
by general statistico-mechanical considerations, 
except that a definite temperature dependence is 
here proposed for the coefficients Qi and Qs. It 
may be noted that Q; and Q are related to the 
coefficients in the “‘virial’”’ equation of state, 


P/RTd=1+61d+Bod?+:--, (B17) 

as follows: 
Bi1=Q,/RT (B18) 
and B2=2Q2/RT. (B19) 


Table I indicates the order of magnitude of 
the errors made by using the simplified equations. 
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Here, the pressure, fugacity, and isothermal 
change in enthalpy of propane at its critical 
temperature and one-half its critical density 
have been calculated with the simplied Eqs. 
(B14), (B15) and (B16) and with the complete 
Eqs. (B7), (B8) and (B11). At lower densities the 
errors resulting from using the simplified equa- 
tions are still smaller, so that for most purposes 
it is permissible to use them when the density is 
less than one-half the critical. On the other hand, 
the simplified equations may not be used for the 
gas phase at higher density or for the liquid 
phase, since the errors rapidly become larger as 
the density increases. 


C. DETERMINATION OF THE NUMERICAL VALUES 
OF THE PARAMETERS FROM 
EXPERIMENTAL DATA 


The foregoing general remarks on the con- 
siderations which led to the form of the present 
equation may be made more concrete by an 
outline of the procedure followed in determining 
from experimental data the most suitable 
numerical values of the parameters for a par- 
ticular hydrocarbon. This procedure will be 
found convenient in fitting the present equation 
to substances other than those treated in this 
article. 

(1) Interpolate experimental P-V-T data for 
the gas phase to even temperatures and den- 
sities. Pressures at these even temperatures and 
densities will be referred to as the observed 
pressures, Pops. 

(2) At each density fit an equation of the 
form: 


P=RTd+d*[RTB(d) ~A(d) —C(d)/T?] (C1) 


to the observed pressures by the method of least 
squares ; give equal weight to each pressure. 

(3) Plot B(d), determined in step (2), against 
density. Pass a straight line through the points; 
give greater weight to the points at the higher 
densities. The equation of this line is: 


B=By+bd. (C2) 


(4) Let the “cohesive pressure,” ops, be 
defined by the equation: 


Tors = RT (d+ Bod? + bd*) — Prous. (C3) 


At each density fit an equation of the form: 
a=d?[A'(d)+C'(d)/T?] (C4) 


to the cohesive pressures by the method of 
least squares ; give equal weight to each pressure. 

(5) Find by trial a value of y such that when 
C’(d) is plotted against d(1+ yd*) exp (— yd’) 
the points fall on an approximately straight line. 
Pass a straight line through these points; give 
greater weight to the points at the higher den- 
sities. The equation of this line is: 


C’(d) = Cy—cd(1+ yd?) exp (— yd”). (C5) 
(6) Let Acs be defined by the equation: 


a’A obs > Tobs 


2 


= [Coed +d") exp (—7d*)]. (C6) 


If B(d) and C’(d) have been properly determined, 
Abs will be practically independent of tempera- 
ture. ‘“‘Best’’ values of A(d), denoted by A’’(d), 
are selected as follows: At densities near the 
critical take the value of A», at the critical tem- 
perature; in this way the equation being fitted 
will be more nearly consistent with observed 
critical properties. At other densities use an 
average value of Aots. 

(7) Find by trial a value of a such that when 
A’'(d) is plotted against d(i—ad*) the points 
fall on an approximately straight line. Pass a 
straight line through these points; give greatest 
weight to the points in the neighborhood of the 
critical density. The equation of this line is: 


A" (d) =Ao+ad(1—ad*). (C7) 


(8) The values of Ao, a and a, given in Eq. 
(C7); Co, c and y, given in Eq. (C5); and Bo and 
b, given in Eq. (C2) constitute provisional values 
of these eight parameters. If these provisional 
values give a satisfactory representation of the 
observed pressures and critical properties, the 
values of b, a, c, y, and a may be accepted as 
final. 

(9) There remains the adjustment of Bo, Ao, 
and Cy, togive the best representation of observed 
vapor pressures. At each temperature, with pro- 
visional values of Bo, Ao, and Cy and final values 
of the other parameters, compute RT Inf at the 
observed vapor pressure for the gaseous and 
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TABLE II. Numerical values of parameters. Units: atmos- 
pheres, liters, g-mole, °K. Atomic weights: C= 12.000, 
H=1.0078. Universal gas constant: R=0.08207. 

T (°K) =t (°C)+273.13. 











METHANE ETHANE PROPANE n-BUTANE 
Bo 0.0426000 0.0627724 0.0973 130 0.124361 
Ao 1.85500 4.15556 6.87225 10.0847 
Co+10-6 0.0225700 0.179592 0.508256 0.992830 
b 0.00338004 | 0.0111220 0.0225000 0.0399983 
a 0.0494000 0.345160 0.947700 1.88231 
c-+10-6 0.00254500 0.0327670 0.129000 0.316400 
Y 0.0060000 0.0118000 0.0220000 0.0340000 
a 0.000 124359 | 0.000243389 | 0.000607175 | 0.00110132 




















liquid phases, as follows: Find by trial a liquid 
density d” and a vapor density d® such that 
when these densities are substituted in equation 
of state (B7), the computed pressures equal the 
observed vapor pressure. Compute the vapor 
fugacity f° and the liquid fugacity f” at d® and 
d", respectively, by means of Eq. (B8). 

By, Ao, and Cy have now to be adjusted so 
that f? and f” at the observed vapor pressure are 
made as nearly equal as possible. The effect on 
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TABLE III. Summarized comparison of observed properties of light hydrocarbons with those predicted by present equation. 


the calculated fugacities of adjusting Bo, Ao, and 
Cy may be estimated through its effect on Q,, 
defined as: 


Q,=BoRT —Ay—C)/T”. (C8) 
It may be proved that 
RT(d In f/0Q1) p=d, (C9) 


so that the effect on f of changing Bo, Ao, and 
Cy by ABo, AAo, and AC is given to the first 
order by: 





RTA |n f=dAQi 
=d(RTABy—AAy—AC)/T?). 





(C10) 





To make f% and f¥ equal at the observed vapor 
pressure, Q; must therefore be changed by an 
amount AQ,, given by: 














RT In f@/ ft 
d"—d* 











METHANE 







AVE. FOR 
4 HypRo- 
CARBONS 








n-BUTANE 





ETHANE PROPANE 





. Pressures of Gaseous Phase 


Source of Data 10, 11, 13 


Range of Density, moles/liter 2 to 18 
Range of Temperature, °C —70 to 200 
Ave. % Dev. of Pressures 0.34 
Observed Normal Density, g/1. 0.7168 
% Deviation —0.01 


. Critical Properties 


Source of Data = 


Deviation of Pressure, atmos. —1.5 —0.17 | —0.42 —0.05 —0.54 

Deviation of Temperature, °C —0.3 —0.25 | —0.16 —0.06 —0.19 

Deviation of Density, moles/liter +0.06 +0.06 0.00 +0.13 +0.06 
3. Properties of Saturated Liquid | 

Source of Data 12 14,17 18,19, 20 20, 22, 23 

Range of Temperature, °C — 139.50 to —85.13 | —50 to 25 —25to75 |4.44to121.11 

Range of Vapor Pressure, atmos. 4.51 to 41.81 5.45 to 41.37| 2.01 to 28.0 | 1.20 to 22.48 

Ave. % Dev. of Densities 2.7 2.7 2.3 1.2 22 

Ave. % Dev. of Vapor Pressures 2.4 0.4 1.0 0.7 ‘A 









13, 14, 15 6, 21 22,24 
0.5 to 10 1 to9 0.5 to 7 
25 to 275 96.81 to 275 150 to 300 
0.31 0.40 0.31 0.34 
1.3566 1.8325 (25°C) 2.70324 
+0.15 —0.005 —0.10 0.07 















15, 16 18 22 
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TABLE IV. Deviations of pressures calculated by equation of state from observed pressures for gaseous n-butane.” 























DENSITY, MOLES/LITER 0.5 1 2 4 5 6 7 8 
AVERAGE 
Temp. °C -PRESSURE, NORMAL ATMOSPHERES % DEV. 
150 Obs. 14.68 | 24.68 | 34.39 
% Dev. —0.07 | —0.04 | —0.06 0.06 
152.01 Obs. 37.46 | 37.76) 43.40] 73.27) 168.11 
% Dev. —0.11 | —0.82 | —0.51 | +0.33 | (—5.39) 0.44 
175 Obs. 15.91 | 27.50} 41.31] 47.89] 52.20] 58.30] 74.32 | 120.59 | 239.49 
% Dev. —0.06 | +0.04 | +0.56 | +0.69 | —0.33 | —1.60 | —0.57 | +0.24 | (—3.87) 0.51 
200 Obs. 17.11 | 30.24] 47.88} 58.99] 68.61} 81.96 | 109.23 | 172.87 | 315.71 
% Dev. —0.12 | +0.07 | +0.56 | +0.66 | —0.20 | —1.05 | —0.13 | +0.38 | (—3.31) 0.40 
225 Obs. 18.30 | 32.92] 54.29] 69.93} 85.20 | 106.23 | 145.02 | 225.65 
% Dev. —0.11 | +0.03 | +0.48 | +0.50 | —0.13 | —0.63 | +0.23 | +0.47 0.32 
250 Obs. 19.48 | 35.55 | 60.62} 80.74 | 101.84 | 130.80 | 181.11 | 278.51 
% Dev. —0.10 0.00 | +0.43 | +0.27 | —0.16 | —0.44 | +0.31 | +0.42 0.27 
275 Obs. 20.65 | 38.18] 66.82} 91.58 | 118.57 | 155.59 | 217.50 | 331.20 
% Dev. —0.10 | +0.05 | +0.31 | +0.17 | —0.19 | —0.37 | +0.30 | +0.25 0.22 
300 Obs. 21.81 | 40.76} 72.97 | 102.28 | 135.27 | 180.56 | 253.96 
% Dev. —0.05 | +0.05 | +0.22 | —0.01 | —0.32 | —0.36 | +0.17 0.17 
Average % Deviation 0.09 0.04 0.37 0.38 0.21 0.75 0.32 0.35 (4.19) 
Over-all Ave. % Dev. 0.31 






































These values of AQ; are to be fitted by an equa- 
tion of the form: 


Find by trial values of ABp, AAo, and ACy such 
that AQ; computed by (C11) is represented as 
well as possible, subject to the restrictions: 
(a) AQ; vanishes at the critical temperature. 
(b) AQ; is made as small as possible between the 
critical temperature and the highest temperature 
at which gaseous pressures are to be computed. 
Requirement (a) ensures that the calculated 
critical properties will not be affected by ad- 
justing Bo, Ao, and Cy: requirement (b) ensures 
that gaseous pressures will be changed as little 
as possible while fitting the vapor pressure data. 

To obtain the final value of Bo, Ao or Co, add 
ABo, AAo or ACy to the provisional value of the 
corresponding parameter. 

(10) Asa check on the process described in the 
preceding paragraphs, observed vapor pressures 
and observed P-V-T properties of the gas phase 
should be compared with those computed with 
final values of the eight parameters. 


D. NUMERICAL VALUES OF THE PARAMETERS 


In Table II are listed numerical values of the 
eight parameters Bo, Ao, Co, 0, a, c, y, and @ for 
the hydrocarbons methane, ethane, propane, and 
n-butane. These values were determined by the 
procedure outlined in Section C from the data 
listed in Table III of Section E, following. These 
data comprise: (1) P-V-T properties of the gas 
and liquid phases up to 1.8 times the critical 
density. (2) Critical properties. (3) Vapor pres- 
sures above one atmos. 

Atomic weights, the value of the universal gas 
constant R and the centigrade temperature of the 
ice point listed in Table II should be used in 
conjunction with the parameters given there, 
because they were employed in fitting the equa- 
tion to the experimental data. 

It is, of course, impossible to assign unique 
values accurate to six significant figures to each 
of these parameters. Nevertheless the full six 
figures given in Table II should be used in cal- 
culating the properties of liquids, and four or 
five should be used for gases. A change of one in 
the sixth place of one constant by itself will alter 
the calculated pressure of a liquid by about 0.01 
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TABLE V. Miscellaneous properties of n-butane. 


1. Properties of Gas at Low Pressure 
Obs. 


% Dev. 
Normal density, g/liter 2.70324 —0.10 
2. Critical Properties 
Obs. Obs.-Calc_ 
Pressure, atmos. 37.48 —0.05 
Temperature, °K 425.14 —0.06 
Density, moles/liter 3.88 +0.13 


3. Properties of Saturated Liquid 


Density, moles/1. Vapor Press., atmos. 


Temp., °C. Obs. % Dev. Obs. % Dev. 
4.44 10.26 +0.6 1.20 +41.2 
37.78 9.66% ++2.2 3.51% —1.0 
71.11 8.93% +1.8 8.23% —0.4 
100.00 8.07% +0.6 15.19% —0.5 
15.092 —1.1 
121.11 7.26% —0.8 22.48% +0.1 


atmos. However, it is possible to change all the 
constants simultaneously by about 5 percent 
without seriously prejudicing the calculated 
results. This point should be borne in mind in 
any attempted correlation of the constants with 
other properties of the hydrocarbons. 


E. COMPARISON OF OBSERVED AND CALCULATED 
PROPERTIES OF LIGHT HYDROCARBONS 


There is not space here for a detailed com- 
parison of observed properties of methane, 
ethane, propane, and n-butane with those pre- 
dicted by the present equation. Table III sum- 
marizes the results of such a comparison and 
indicates the source of the data used. Since the 
predictions of the equation correspond about 
equally well with experimental data for each 
hydrocarbon, the detailed comparison for n-bu- 
tane, given in Tables IV and V, may be con- 
sidered as typical of the results for methane, 
ethane, and propane as well. In these tables, 
“&% Dev.” is defined by 


‘ai Obs. Property — Calc. Property 





Obs. Property 


and ‘‘Ave. % Dev.” is the average of the absolute 
values of the individual deviations. 

Tables IV and V illustrate the agreement 
between observed data for n-butane and ‘values 
calculated by the equation of state. Similar 
tables for methane, ethane, and propane give 
about as good agreement between observed and 
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calculated properties as here illustrated for 
n-butane. A comparison between the columns in 
Table III helps to show this point. 

It was impossible to fit the points for x-butane 
(Table IV) at 8 moles/liter, about twice the 
critical density, and these were not used in deter- 
mining the parameters. Excluding these three 
points, the average deviation of the calculated 
pressures is 0.31 percent. Up to 4 moles/liter the 
average deviation is only 0.22 percent, whereas 
over approximately the same density range the 
average deviation from the Beattie-Bridgeman 
equation” is 0.36 percent. It may be said that the 
equation fits the data well for all four hydro- 
carbons up to about 1.8 times the critical density, 
but above this density calculated pressures are 
high. 

In Part 2 calculated critical properties are 
shown to be in good agreement with those ob- 
served by Beattie, Simard and Su.” Table III, 
Part 2, shows fair agreement for the other three 
hydrocarbons. Calculated critical pressures are 
somewhat high, though with the exception of 
methane the deviations are less than 0.5 atmos. 
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Calculated critical temperatures are from 0.06 to 
0.3°C high. 

In Part 3 of Table V calculated vapor pressures 
of n-butane are shown to agree with observed 
vapor pressures practically within experimental 
error. Calculated liquid densities differ from 
observed liquid densities with an average devi- 
ation of 1.2 percent. For the other hydrocarbons 
as well as for n-butane calculated liquid den- 
sities are from 0 to 3 percent low. Vapor pressures 
agree well down to pressures of the order of 2 
atmos.; at lower pressures calculated vapor 
pressures are low. 

The average deviations summarized in Table 
III indicate the error resulting when the present 
equations are used to compute thermodynamic 
properties of these hydrocarbons within the 
regions of density, pressure, and temperature in 
which the equations have been fitted. Caution 
should be exercised in using the present equa- 
tions at higher densities or lower temperatures. 
Extrapolation to higher temperatures or lower 
densities is less hazardous. 

The use of the equation for calculating thermo- 
dynamic properties other than those considered 
in evaluating the parameters is illustrated by 
the following examples: (1) the effect of pressure 
on the enthalpy of methane, and (2) the latent 
heat of evaporation of n-butane. 

Eucken and Berger®® have determined the 
effect of pressure on the enthalpy of methane 
calorimetrically below 20°C, and Budenholzer, 
Sage and Lacey®* have computed it from Joule- 
Thomson coefficients measured above 21.11°C. 
In Fig. 5 experimental values of H(0, T) 
—H(P,T) at —81.8, —50, and 20°C are com- 
pared with those calculated by Eqs. (B7) and 
(B11). The agreement at —81.8°C, the computed 


us and Berger, Zeits. f. ges. Kalte-ind. 41, 145 
% Budenholzer, Sage and Lacey, Ind. Eng. Chem. 31, 
369 (1939). 


critical temperature, and at 20°C, is probably 
within the error of the measurements. At — 50°C 
the agreement is less satisfactory. 

The calculated latent heat of evaporation of 
n-butane at 37.78°C is compared below with that 
measured calorimetrically by Sage, Webster and 
Lacey :** 


Calculated 202.6 
Observed 201.1 


The calculated latent heat is the difference in 
enthalpy between the liquid and vapor, each 
evaluated by Eq. (B11) at the calculated satura- 
tion density. 

The minor discrepancies between observed 
properties of pure hydrocarbons and those cal- 
culated by the present equation, although real, 
do not seriously impair the value of the present 
equation. It provides in a single expression a 
concise and reasonably accurate summary of 
several classes of thermodynamic data pertaining 
to hydrocarbons and covers a wider range of data 
than any equation proposed up to the present 
time. It may be applied both to the gas phase and 
to the saturated liquid. Perhaps the greatest 
value of the equation is that it provides a point 
of departure for the prediction of thermodynamic 
properties of mixtures. 

The next paper of this series will show how the 
present equation may be generalized to describe 
the thermodynamic properties of mixtures. Work 
is also planned on extending the equation to 
substances other than those treated in this first 


paper. 


liter-atmos. /mole 
liter-atmos. /mole 
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NYONE who has had much to do with 
special and general point positions in space 
groups of high symmetry has been delayed by 
the rather awkward notation in available tables, 
the best of which are the “International tables 
for the determination of crystal structures.’’! The 
coordinates of the different points of a single 
set, when written on the typewriter, take from 3 
characters as a minimum to 11 as a maximum; 
e.g., xyz; }—x, }—y, }—=. In addition it may be 
remarked that mistakes in copying fractions, 
printed or typed, are especially easy and that the 
derivation of point coordinates by application of 
symmetry operations is hindered by an unsym- 
metrical notation. This is written to point out 
that a consistent notation with exactly 6 charac- 
ters per point is easy to devise and remember. It 
is suggested that future tabulations take ad- 
vantage of this simplification. This suggestion 
was, in fact, made some years ago,’ but has not 
hitherto been published. 

The only fractions which appear with arbitrary 
parameters added to or subtracted from them are 
4%, \%, %, %, 24, 34, 36. Besides these there 
appear without additions the fractions 14, 3¢, 5¢, 
7g, and without fractional parts the multiplier 2 
affecting an arbitrary parameter. The problem is 
to express all these in exactly two typed charac- 
ters so that six consecutive characters will always 
be necessary and sufficient. 

Adopting a duodecimal basis and suppressing 
the denominator (10 on this basis, 12 on a 
decimal basis) the fractions to which parameters 
are added or subtracted become 2, 3, 4, 6, 8, 9, ¢ 
(using ¢ for “‘ten’’). The fractions that appear 
alone are similarly 16, 46, 76, 46. (Denominator 
100 suppressed.) The addition of two fractions 
can now be represented merely by writing them 
together, 3u for 3+, and as is customary in 


1 Leading title: Internationale Tabellen zur Bestimmung 
eV (Gebriider Borntraeger, Berlin, 
1935). 

2 At a meeting of the American Academy of Arts and 
Sciences, October 8, 1935. 
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some degree already, a negative quantity can be 
shown by overlining, 3a for 3—u. The multiplier 
2 is consistently suppressed, the symbol for the 
parameter being written twice, ## for 2%. 
Wherever a space would otherwise appear a zero 
is inserted. This completes the system of notation. 
It may be remarked that the letters chosen for 
representing parameters should be those which 
are not made unsightly or ambiguous by over- 
lining or underlining (the latter will be preferred 
by typists). The only conventional set of three 
consecutive letters which meets this criterion is 
the set u, v, w. If it is desired to use different 
letters for one, two, and three parameter sets the 
best choices are m; 7, 5; u, v, w. Another way of 
designating parameters without occupying more 
space than is available is to use superscript 
letters (or numbers 1, 2, 3) for added fractions, 
subscript letters (or numbers) for subtracted 
fractions. This is easier in typing but has the 
disadvantage of using a small type size in printing. 

As examples two point position sets are copied 
below in the present standard notation and 
translated into the suggested consistent notation. 


Space Group D344—H31c Point position 18(f) 


old (000; 3 30; 33 0)+ 
x £4; x, 2x, 1; 2%, Z, 4; 
#24; 2, 2%, 4; 2x, 2, 7. 

new (000000; 408000; 804000) + 


0u0730 ; Ouuu30 ; 7u0030 ;s 
070u90 ; 072%090 ; uu0u90. 


Space Group O1— F4,3 Point position 48(f ) 


old (000;033;203;330)+ 
4, x,4—x;4—x, 3,4; 2,4 —2,4;4, 2, 342; 
EEX, $125 E, FAH, $5 81H, GK; THX, 8%) 
x, 4+, 4; $1, g—X; g—X, §, 2; 2, 5—%, 3. 


new (000000; 006060; 600060; 606000) + 
160u3a ; 3%160u ; 0u3016; 160029u; 
9u1604 ; 027916 ; 460u3u; 3u460u; 
0u3u46 ; 4602924 ; 974602 ; 079746. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section should not exceed 600 words in length and must 


Transformations in Monolayers 


It appears from the letter to the Editor by Harkins and 
Boyd! that many parts of my article? have been overlooked. 

Points of transformation.—The thermodynamical defini- 
tions given on p. 932 are valid however indefinite the limits 
of the surface phase may be. 

Mesomorphous phase——This is analogous to a liquid 
crystal, not because it lies between the solid and liquid 
states, but for reasons and experimental evidence which 
have been given in a special paragraph, p. 945. 

Liquid phase——In other temperature ranges, the con- 
densation from the gaseous to the solid or mesomorphous 
phases is also a first-order transformation. This is not 
peculiar only to the liquid. As for the compressibility, there 
can be no comparison with the three-dimensional liquid, 
since in the film all the molecules are oriented with refer- 
ence to the same plane and the compressibility is measured 
only in that plane. The heat expansion of three-dimensional 
paraffin crystals is incomparably smaller in the direction 
of the chains than in the two other directions of the unit 
cell.3 

Isotropic film.—We are enabled to estimate the degree 
of anisotropy of liquid films by comparing the measured 
areas at different points of the isotherms to the area 
A=V?/, Actually it happens that, when T=T.., the areas 
at the point M are close to the calculated values. This is 
an experimental fact, an explanation of which I have 
tried to give, p. 941. . 

State of a film.—This is defined when three of the four 
variables A, x, T, pH, are fixed. Part of the misunderstand- 
ing of Harkins and Boyd comes from the fact that they 
compare molecular areas without defining the temperature 
and the pressure. Thus they compare to the numbers of 
column 4 of my Table I, which are relative to molecular 
areas of the liquid film under its vapor pressure and at 
temperatures T=T-c, values found for T7<7T..- (column 2, 
Table II of H. and B.) and even areas at the point 0 
(transition to the erystal). The value of 57A? for myristic 
acid is in agreement with Adam and Jessop’s classical 
isotherms (7 =34°), while the smaller area of 45A? may 
be found at 17° (see Fig. 17). 

General diagram.—My diagram has been incompletely 
and wrongly reproduced by Harkins and Boyd. Especially 
the isotherm relative to T>T.- is missing, while all the 
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reach the office of the Managing Editor not later than the 15th 
of the month preceding that of the issue in which the letter 
is to appear. No proof will be sent to the authors. The usual 
publication charge ($3.00 per page) will not be made and no 
reprints will be furnished free. 


discussion on the isotropy at the point M concerns pre- 
cisely that case. 

Areas in the solid and mesomorphous states —The state 
of the solidified film between 19.5 and 18.5? is metastable. 
This can be easily detected bv registration of the isotherm 
(see Figs. 7 and 14 and read legend); when using the ordi- 
nary method, there is no wonder if the film seems to col- 
lapse at 19.5A*. 

The areas of 18.5, 19.5, 20.5, 23.5A? per chain are average 
values. With the actual technique it is impossible to esti- 
mate the molecular area with a precision greater than 
+1 percent (errors on the area, volume of the solvent 
deposited, weighing of the substance). With the exception 
of a few studies on single crystals, x-ray analysis gives 
data with the same indetermination. So that even the 
values given by H. and B. in columns 3 and 4 agree in 
majority with the x-ray data. Besides, it should be kept 
in mind that palmitic acid at 20° is already in the ex- 
panded mesomorphous form (7;,=18°) and that the area 
of 23.5 for stearic acid is found at room temperature on 
substrates between distilled water and HCI N/500, while 
it rises to 24.5 when reaching NV/100. The important point 
is that, even when the area at the sublimation point is 
larger than 24A?, the value of 23.5 unquestionably appears 
at higher pressures as a second-order point (p. 946). 

Concerning alcohols, although their characteristic 
areas seem to agree with those for acids (see H. and B.’s 
own values), it would be more reasonable to compare them 
with those found by x-ray analysis of the same alcohols. 

Van der Waals forces do act by raising the melting and 
critical temperatures. But it seems doubtful whether, 
in the same state, the packing of the molecules increases 
with the chain length. This is a matter of discussion based 
on such feeble evidence as variations of less than 0.5 per- 
cent in the measured areas. For the moment, we are able 
to say that, within the precision of both methods, there is 
correspondence between the x-ray and monolayer meas- 
urements. 

D. G. DERVICHIAN 
Laboratoire de Chimie Physique, 
Faculté des Sciences, 


Paris, France, 
March 13, 1940. 


! Harkins and Boyd, J. Chem. Phys. 8, 129 (1940). 
2? Dervichian, J. Chem. Phys. 7, 931-48 (1939). 
3 Miiller, Proc. Roy. Soc. Al27, 417 (1930). 
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Mechanism of the Exchange Reaction Between Gaseous 
Bromine and Hydrogen Bromide 


Liberatore and Wiig! have presented conclusive evidence 
for a rapid homogeneous interchange reaction between Br2 
and HBr in the gaseous state at room temperatures. 
Their data show the reaction to be nearly but not quite 
at equilibrium in two minutes and to be homogeneous and 
nonphotochemical. As a mechanism, the authors suggest a 
chain involving the production of Br atoms by the ioniza- 
tion processes associated with the radioactive bromine 
radiation, the Br atoms reacting according to 


Br+HBr* = HBr-+Br* (1) 
or Br+BrBr* = BrBr+Br* (2) 


to establish uniform distribution of the radioactive isotopes. 

In this connection it may be of interest to consider ex- 
periments performed under conditions nearly identical with 
those of Liberatore and Wiig with the principal difference 
that the intensity of the radioactive radiations was nearly 
20,000-fold smaller. The actual number of radioactive 
atoms present during a typical experiment was 1.5105 
while Liberatore and Wiig give 3.5X10° as the value for 
their experiments. At these lower activities the reaction was 
within ten percent of completion in two minutes and was 
not appreciably affected either by painting the reaction 
vessel black or coating the inner surface with paraffin. 
For these reasons it seems certain that the suggested 
mechanism is not the principal one involved in the reac- 
tion. The alternative mechanism mentioned by the authors, 
a direct reaction between Br2 and HBr through a HBrs 
complex, appears to be the correct one. 

The data require that this reaction have an activation 
energy less than 11 kilocalories. Such a value does not 
appear to be entirely unreasonable. 

Additional experiments have shown that I, and HI 
undergo a similarly fast reaction both in the gaseous state 
and when dissolved in dry pentane. 


W. F. Lipsy 
Department of Chemistry, 
University of California, 
Berkeley, California, 
February 16, 1940. 


1L. C. Liberatore and E. O. Wiig, J. Chem. Phys. 8, 165 (1940). 





Thermal Diffusion Separation of Different Gases of the 
Same Molecular Weight 


Most of the work done by means of the recently de- 
veloped thermal diffusion method’ has been directed 
toward the separation of isotopes and isotopic compounds. 
In order to find experimentally some of the factors other 
than mass which influence separations, we have applied 
the thermal diffusion method to the separation of different 
gases of the same molecular weight (to the nearest whole 
mass unit). In such cases any separation which might occur 
would be due to the difference in structure of the molecules 
and to their different forces of attraction and repulsion. 

The apparatus used was a brass column nine feet long 
with an outer jacket cooled by tap water at 20°C and an 
inner tube heated by means of No. 22 Chromel wire to 
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TABLE I. 








COMPOSITION IN MOLE PERCENT 





AFTER DIFFUSION 
DIFFERENCE 
+0.4 


BEFORE 
DIFFUSION 
. COMBINATION BoTTOM 


- COz and C3Hs 





34.4% COz 
53.6 


. COz and N20 


69.1 
36.1% COz 
42.5 
54.5 
70.1 
26.8% CO 
38.0 
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. CO and Nz 
. CO and C2H, 


- Ne and CoH, 
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about 225°C. Mixtures of gases were placed into an inter- 
mediate chamber which was 0.7 cm thick and had an 
inner radius of 1.0 cm. The combinations of gases used 
were carbon dioxide and propane, nitrogen and ethylene, 
etc., the runs being made for various concentrations. 

The results of the work are summarized in Table I. 
The percentage compositions after thermal diffusion are 
the averages of at least four measurements. For each 
of the combinations listed except CO and Ng, the separa- 
tions were quite substantial. In the case of CO and N2, 
the molecules are so similar that no separation was found, 
the observed deviations being within the experimental 
error. The probable error in each case was about +0.4 
percent. 

Remarkably enough, carbon dioxide could be partially 
separated from nitrous oxide, the carbon dioxide tending 
to go to the top. No explanation for this is apparent, at 
least not when compared to the results for CO and Nz. In 
the other cases the larger of two molecules of the same 
molecular weight was found to go to the bottom. This 
result is probably to be related to the fact that the bulkier 
molecules have lower inverse power repulsive forces than 
do the compact molecules. 

The slight mass differences between carbon dioxide and 
propane or between ethylene and carbon monoxide or 
nitrogen are not sufficient to account for the observed 
separations. Even if the masses differed by a whole atomic 
unit, the separations on that account would not be as 
large as observed for a column only nine feet long. Ac- 
cordingly, the separation must be attributed to differences 
in structure. 

Experimental studies of other combinations are under 
way. Full details of the work will be published shortly. 
F. T. WALL 
C. E. Hottey, Jr. 


Noyes Chemical Laboratory, 
University of Illinois, 
Urbana, Illinois, 


February 29, 1940. 


(935). Clusius and G. Dichel, Naturwiss. 26, 546 (1938); 27, 148, 487 
2A, Bramley and A. K. Brewer, J. Chem. Phys. 7, 553, 972 (1939); 
Science 90, 165 (1939). 
3W. H. Furry, R. C. Jones and L. Onsager, Phys. Rev. 55, 1083 
(1939). (Theoretical paper.) 
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Windows for Studying the Optical Properties 
of Hydrogen Fluoride Gas 


Following the remarks of Freed, McMurry, and Rosen- 
baum! on the use of artificial corundum in optics, we should 
emphasize that the preliminary researches, interrupted in 
August by the present events, have led us to certain con- 
clusions relative to the materials which could constitute 
windows resisting the action of HF. These researches have 
as their aim, extension of our results regarding liquefied 
HCI to hydrogen fluoride.? 

In dilute solution, certain organic glasses, in particular 
Duroid, are inattackable by HF; those which are colorless 
may serve in the visible spectrum; they are deprived of 
double refraction and can be used in polarimetry, as Mlle 
F. Perrier has demonstrated. Duroid is yellow colored, but 
does not show in infra-red (under about 4 mm) important 
bands up to 1.2y; it allowed us to discover a wide and weak 
band about 0.93, distinct from that of the gas, in the aqueous 
solution of hydrogen fluoride (about 40 percent in weight). 

In contact with a concentrated solution of HF, gas or 
pure liquid, all organic substances which we have tried, 
have been attacked (probably by dehydration); the 
fluorite itself tarnishes more or less quickly, especially at a 
high temperature (perhaps by the formation of CaF», HF). 
We have taken instead artificial white corundum (white 
sapphire), of which we have been able to obtain, in France, 
some plates of 22 mm diameter, parallel to the axis. 
We have verified that this substance does not show bands 
of absorption between 1850A and 2.5y, and that it is not 
attacked by the gas HF, even at high temperatures (our 
experiments were made during an hour at 300°C). 

B. VoDAR 
R. FREYMANN 


Yreou Ta 
Laboratoires de Physique de la Sorbonne, 
Paris, France, 
January 20, 1940. 


1 Freed, McMurry and Rosenbaum, J. Chem. Phys. 7, 853 (1940). 
2 J. de phys. et rad. 9, 282 (1939). 





Exchange Reactions Between Gaseous Alkyl Bromides, 
Bromine and Hydrogen Bromide 


An investigation, with the aid of radioactive bromine, 
of the exchange reactions of ethyl bromide with bromine 
and with hydrogen bromide confirms our recent conclu- 
sion! that the exchange reaction between bromine and 
hydrogen bromide proceeds by an atom chain. That this 
last exchange is a bimolecular reaction? going through a 
HBrs complex and requiring an energy of activation less 
than 11 kcal. appears to be highly improbable. 

Gaseous Br2 containing radioactive bromine when 
mixed with ethyl bromide for various periods of time up 
to two hours shows no exchange at all. Similar results have 
been obtained with HBr and C;H;Br. The mixture of Br. 
and C,H;Br was also exposed to the full light of a 500-watt 
tungsten lamp for 9 hours without any detectable ex- 
change at room temperatures. When heated, exchange sets 
in. These results are all in agreement with an atom mecha- 
nism. Rapid exchange takes place when HBr and C.H;Br 
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are heated to 200-300°, but this probably occurs through 
the equilibrium C,H;Br—@C,H,+HBr. 

If the exchange reaction between HBr and Bro, which is 
very rapid at room temperatures, proceeds through a 
bimolecular reaction as Libby? suggests, then exchange 
should also occur between C,H;Br and Brz, since the 
strength of the C—Br bond (ca. 65 kcal.) is less than that 
of the H—Br bond (ca. 86 kcal.). Calculations by a semi- 
empirical method developed by Hirschfelder? of the energies 
of activation for these two reactions yield values of 37 
and 31 kcal., respectively. Similar calculations for the Br 
atom exchange reactions with C.H;Br and with HBr yield 
values of about 25 and less than 5 kcal., respectively. 
Thus, of these four reactions only the HBr+Br atom ex- 
change would be expected to occur at room temperatures, 
in agreement with our experiments. 

The atom chain reaction appears to be the only possible 
mechanism, though the Br atom concentrations in Libby's 
experiments might seem to be too small. In our HBr+Br. 
experiments the Br atom concentration resulting from 
radioactive processes was larger (30 : 1) than that due to 
thermal dissociation of bromine. The latter is much larger 
than the former in Libby’s experiments, however. The Br 
atom concentration at room temperatures may be estimated 
by extrapolation of the data of Perman and Atkinson‘ 
at 650-1050° or of Bodenstein and Cramer at 800-1200°. 

Calculations indicate that to account for Libby’s results 
under our experimental conditions in two minutes about 
225 collisions per cc of HBr with Br atoms would be re- 
quired, whereas from the data of Perman and Atkinson 
the thermal Br atom concentration is sufficient to yield 
only 10 collisions. Bodenstein and Cramer’s data give even 
fewer collisions. However, such an extrapolation might 
easily be too low by a factor of 100 or more. Thus Libby’s 
results may also be accounted for by an atom chain reac- 
tion initiated by bromine atoms produced thermally. 

More complete details of these investigations of the 
exchange reactions of Br. and HBr with C.H;Br and other 
alkyl bromides will appear shortly. 

L. C. LIBERATORE 
E. O. Wic 


University of Rochester, 
Rochester, New York, 
March 14, 1940. 


1 Liberatore and Wiig, J. Chem. ae. 8, 165 (1940). 

2 Libby, J. Chem. Phys. 8, 348 (1940). 

3 University of Wisconsin Symposium on the Kinetics of Homo- 
geneous Gas Reactions, June, 1939. 

*Lewis and Randall, Thermodynamics (McGraw-Hill Book Co., 
New York, 1923), p. 514. 

5 Bodenstein and Cramer, Zeits. f. Elektrochem. 22, 327 (1916). 





The Absorption Spectra of Liquid HF and of its Aqueous 
Solution in the Region \\8000—11,000 


The absorption spectra of liquid HF and of HF in 
aqueous solution have been photographed in the region 
4A8000-11,000 with a dispersion of about 60A per mm. 
In the case of the anhydrous HF a copper cell was used, 
with windows of MgO (artificial crystal) clamped directly 
onto the ends of the tube. The MgO is not attacked by 
anhydrous HF and appears to have greater mechanical 
strength than glass. 
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The absorption band of anhydrous HF in the liquid 
state is very similar to the band of liquid water at \9722 
as reported by Badger and Bauer,! though it is somewhat 
broader. The position of the band maximum determined at 
several temperatures is given in the following table. The 
large probable error indicated arises from the difficulty 
in estimating the maximum of a broad band in this region 
where the trend in plate sensitivity is large. 

TABLE I. The maximum of the band of liquid HF as a 

function of temperature. 

#(°C) Amax(A) 
19 9940 +75 
0 9990 

—30 10140 

—60 10300 

75 10380 

Salant and Kirkpatrick? have found that the third 
harmonic band of unassociated HF gas is located at \8790. 
From the data just given it is evident that the frequency 
decrease in going from the gas to the liquid is very large 
even at the boiling point, namely about 1320 cm™, This 
shift in the HF frequency, which must be due mainly to the 
formation of hydrogen bonds of the type F—H---F, is 
large in comparison with the shift in OH frequency on the 


THE EDITOR 


formation of the O—H---O linkage. The third harmonic 
OH band of methyl alcohol, for example, is found at 
49490 in the vapor* and at 10,070 in the liquid,! so that the 
frequency shift is about 600 cm7. 

The large difference in the magnitudes of the frequency 
shifts just mentioned is very interesting in view of the 
fact that the energies of the hydrogen bonds in associated 
HF and CH;0H are so nearly the same, being 6.7 and 
6.2 kcal. per mole,‘ respectively. 

The absorption spectrum of HF in aqueous solution has 
also been examined at HF concentrations of 25 and 50 
percent and at 0° and at 25°C. In no case does a cursory 
examination reveal any significant difference between the 
solution band and the 9722 band of water in either 
position or shape. 

AusTINn L. WAHRHAFTIG 


Gates and Crellin Laboratories of Chemistry, 
California Institute of Technology, 
Pasadena, California, 
March 11, 1940. 


1 Badger and Bauer, J. Chem. Phys. 5, 840 (1937). 

2 Salant and Kirkpatrick, Phys. Rev. 48, 945 (1935). 

3 Badger and Bauer, J. Chem. Phys. 4, 469 (1936). 

4L. Pauling, The Nature of the Chemical Bond (Cornell University 
Press, 1939), p. 313. 








